Efficient and Robust Parameter Optimization of the Unitary Coupled-Cluster Ansatz.

The variational quantum eigensolver (VQE) framework has been instrumental in advancing near-term quantum algorithms. However, parameter optimization remains a significant bottleneck for VQE, requiring a large number of measurements for successful algorithm execution. In this paper, we propose sequential optimization with approximate parabola (SOAP) as an efficient and robust optimizer specifically designed for parameter optimization of the unitary coupled-cluster ansatz on quantum computers. SOAP leverages sequential optimization and approximates the energy landscape as quadratic functions, minimizing the number of energy evaluations required to optimize each parameter. To capture parameter correlations, SOAP incorporates the average direction from previous iterations into the optimization direction set. Numerical benchmark studies on molecular systems demonstrate that SOAP achieves significantly faster convergence and greater robustness to noise compared with traditional optimization methods. Furthermore, numerical simulations of up to 20 qubits reveal that SOAP scales well with the number of parameters in the ansatz. The exceptional performance of SOAP is further validated through experiments on a superconducting quantum computer using a 2-qubit model system.

Toward Chemical Accuracy with Shallow Quantum Circuits: A Clifford-Based Hamiltonian Engineering Approach.

Achieving chemical accuracy with shallow quantum circuits is a significant challenge in quantum computational chemistry, particularly for near-term quantum devices. In this work, we present a Clifford-based Hamiltonian engineering algorithm, namely CHEM, that addresses the trade-off between circuit depth and accuracy. Based on a variational quantum eigensolver and hardware-efficient ansatz, our method designs the Clifford-based Hamiltonian transformation that (1) ensures a set of initial circuit parameters corresponding to the Hartree-Fock energy can be generated, (2) effectively maximizes the initial energy gradient with respect to circuit parameters, (3) imposes negligible overhead for classical processing and does not require additional quantum resources, and (4) is compatible with any circuit topology. We demonstrate the efficacy of our approach using a quantum hardware emulator, achieving chemical accuracy for systems as large as 12 qubits with fewer than 30 two-qubit gates. Our Clifford-based Hamiltonian engineering approach offers a promising avenue for practical quantum computational chemistry on near-term quantum devices.

Time-Dependent Density Matrix Renormalization Group Method for Quantum Transport with Phonon Coupling in Molecular Junction.

Quantum transport in molecular junctions has attracted great attention. The charge motion in a molecular junction can cause geometric deformation, leading to strong electron phonon coupling, which was often overlooked. We have formulated a nearly exact method to assess the time-dependent current and occupation number in the molecular junction modeled by the electron-phonon coupled bridge state using the time-dependent density matrix renormalization group (TD-DMRG) method. The oscillation period and amplitude of the current are found to be dependent on the electron phonon coupling strength and energy level alignment with the electrodes. In an attempt to better understand these phenomena, we have devised a new approximation that explains the bistability phenomenon and the behavior of steady currents in the strong electron-phonon coupling regime. Comparisons have been made with the multilayer-multiconfiguration time-dependent Hartree (ML-MCTDH) method and the analytical result in the purely electronic limit. Furthermore, we explore the entropy of different orderings, extending to the electron phonon model problems. Regarding finite temperature, the thermal Bogoliubov transformation of both fermions and bosons is used and compared with imaginary time evolution results.

Rationale for white skin roll flap in unilateral cleft lip repair: A retrospective anthropometric measurement analysis.

BACKGROUND: Although the "white skin roll" of the lip is often considered a line, it is better defined as the subunit between the vermilion border and the upper lip horizontal groove. In many unilateral cleft lip repair techniques, this structure is approximated between both sides of the cleft without restoration. This study aimed to analyze the white skin roll height in patients with unilateral cleft lip. METHODS: This retrospective cohort study included 134 consecutive infants with unilateral cleft lip aged 3-6 months who underwent lip repair in a single institution between January 2019 and July 2021. White skin roll heights at the peak of the Cupid's bow on the non-cleft side (CPHIR), cleft medial element (CPHIL), and cleft lateral element (CPHIL') were measured, and differences in their averages were analyzed. RESULTS: The mean height was 1.70 ± 0.30 mm at CPHIR, 0.98 ± 0.33 mm at CPHIL, and 1.28 ± 0.32 mm at CPHIL.' The mean difference in height between CPHIR-CPHIL, CPHIR-CPHIL,' and CPHIL-CPHIL' groups was significant for each paired sample (p < 0.01). No difference was found between the complete and incomplete clefts or left and right clefts (p > 0.01). CONCLUSIONS: A significantly reduced mean height of the white skin roll was present more markedly on the cleft medial element than on the cleft lateral element. Therefore, we strongly support using a white skin roll flap on the cleft lateral element for unilateral cleft lip repair, embracing the concepts of subunits and lip contour lines.

TenCirChem: An Efficient Quantum Computational Chemistry Package for the NISQ Era.

TenCirChem is an open-source Python library for simulating variational quantum algorithms for quantum computational chemistry. TenCirChem shows high-performance in the simulation of unitary coupled-cluster circuits, using compact representations of quantum states and excitation operators. Additionally, TenCirChem supports noisy circuit simulation and provides algorithms for variational quantum dynamics. TenCirChem's capabilities are demonstrated through various examples, such as the calculation of the potential energy curve of H2O with a 6-31G(d) basis set using a 34-qubit quantum circuit, the examination of the impact of quantum gate errors on the variational energy of the H2 molecule, and the exploration of the Marcus inverted region for charge transfer rate based on variational quantum dynamics. Furthermore, TenCirChem is capable of running real quantum hardware experiments, making it a versatile tool for both simulation and experimentation in the field of quantum computational chemistry.

Feasibility and Toxicity of Interval-Compressed Chemotherapy in Asian Children and Young Adults with Sarcoma.

Twelve Asian patients with sarcoma received interval-compressed (ic-) chemotherapy scheduled every 14 days with a regimen of vincristine (2 mg/m2), doxorubicin (75 mg/m2), and cyclophosphamide (1200-2200 mg/m2) (VDC) alternating with a regimen of ifosfamide (9000 mg/m2) and etoposide (500 mg/m2) (IE), with filgrastim (5-10 mcg/kg/day) between cycles. Carboplatin (800 mg/m2) was added for CIC-rearranged sarcoma. The patients were treated with 129 cycles of ic-VDC/IE with a median interval of 19 days (interquartile range [IQR], 15-24 days. Median nadirs (IQR) were neutrophil count, 134 (30-396) × 106/L at day 11 (10-12), recovery by day 15 (14-17) and platelet count, 35 (23-83) × 109/L at day 11 (10-13), recovery by day 17 (14-21). Fever and bacteremia were observed in 36% and 8% of cycles, respectively. The diagnoses were Ewing sarcoma (6), rhabdomyosarcoma (3), myoepithelial carcinoma (1), malignant peripheral nerve sheath tumor (1), and CIC-DUX4 Sarcoma (1). Seven of the nine patients with measurable tumors responded (one CR and six PR). Interval-compressed chemotherapy is feasible in the treatment of Asian children and young adults with sarcomas.

Toward practical quantum embedding simulation of realistic chemical systems on near-term quantum computers.

Quantum computing has recently exhibited great potential in predicting chemical properties for various applications in drug discovery, material design, and catalyst optimization. Progress has been made in simulating small molecules, such as LiH and hydrogen chains of up to 12 qubits, by using quantum algorithms such as variational quantum eigensolver (VQE). Yet, originating from the limitations of the size and the fidelity of near-term quantum hardware, the accurate simulation of large realistic molecules remains a challenge. Here, integrating an adaptive energy sorting strategy and a classical computational method-the density matrix embedding theory, which respectively reduces the circuit depth and the problem size, we present a means to circumvent the limitations and demonstrate the potential of near-term quantum computers toward solving real chemical problems. We numerically test the method for the hydrogenation reaction of C6H8 and the equilibrium geometry of the C18 molecule, using basis sets up to cc-pVDZ (at most 144 qubits). The simulation results show accuracies comparable to those of advanced quantum chemistry methods such as coupled-cluster or even full configuration interaction, while the number of qubits required is reduced by an order of magnitude (from 144 qubits to 16 qubits for the C18 molecule) compared to conventional VQE. Our work implies the possibility of solving industrial chemical problems on near-term quantum devices.

Computational Method for Evaluating the Thermoelectric Power Factor for Organic Materials Modeled by the Holstein Model: A Time-Dependent Density Matrix Renormalization Group Formalism.

Organic/polymeric materials are of emerging importance for thermoelectric conversion. The soft nature of these materials implies strong electron-phonon coupling, often leading to carrier localization. This poses great challenges for the conventional Boltzmann transport description based on relaxation time approximation and band structure calculations. In this work, combining the Kubo formula with the finite-temperature time-dependent density matrix renormalization group (FT-TD-DMRG) in the grand canonical ensemble, we developed a nearly exact algorithm to calculate the thermoelectric power factor PF = α2 σ, where α is the Seebeck coefficient and σ is the electrical conductivity, and apply the algorithm to Holstein Hamiltonian with electron-phonon coupling to model organic materials. Our algorithm can provide a unified description covering the weak coupling limit described by the bandlike Boltzmann transport to the strong coupling hopping limit.

Hybrid Quantum-Classical Boson Sampling Algorithm for Molecular Vibrationally Resolved Electronic Spectroscopy with Duschinsky Rotation and Anharmonicity.

Using a photonic quantum computer for boson sampling has demonstrated a tremendous advantage over classical supercomputers. It is highly desirable to develop boson sampling algorithms for realistic scientific problems. In this work, we propose a hybrid quantum-classical sampling (HQCS) algorithm to calculate the optical spectrum for complex molecules considering Duschinsky rotation effects and anharmonicity. The classical sum-over-states method for this problem has a computational complexity that exponentially increases with system size. The HQCS algorithm creates an intermediate harmonic potential energy surface (PES) to bridge the initial and final PESs. The magnitude and sign of the overlap between the initial and the intermediate state are estimated by boson sampling and classical algorithms, respectively. The overlap between the intermediate and the final state is efficiently evaluated by classical algorithms. The feasibility of HQCS is demonstrated in calculations of the emission spectrum of a Morse model as well as the pyridine molecule.

On the fly swapping algorithm for ordering of degrees of freedom in density matrix renormalization group.

Density matrix renormalization group (DMRG) and its time-dependent variants have found widespread applications in quantum chemistry, includingab initioelectronic structure of complex bio-molecules, spectroscopy for molecular aggregates, and charge transport in bulk organic semiconductors. The underlying wavefunction ansatz for DMRG, matrix product state (MPS), requires mapping degrees of freedom (DOF) into a one-dimensional topology. DOF ordering becomes a crucial factor for DMRG accuracy. In this work, we propose swapping neighboring DOFs during the DMRG sweeps for DOF ordering, which we term 'on the fly swapping' (OFS) algorithm. We show that OFS is universal for both static and time-dependent DMRG with minimum computational overhead. Examples are given for one dimensional antiferromagnetic Heisenberg model,ab initioelectronic structure of N2molecule, and the S1/S2internal conversion dynamics of pyrazine molecule. It is found that OFS can indeed improve accuracy by finding better DOF ordering in all cases.

A general charge transport picture for organic semiconductors with nonlocal electron-phonon couplings.

The nonlocal electron-phonon couplings in organic semiconductors responsible for the fluctuation of intermolecular transfer integrals has been the center of interest recently. Several irreconcilable scenarios coexist for the description of the nonlocal electron-phonon coupling, such as phonon-assisted transport, transient localization, and band-like transport. Through a nearly exact numerical study for the carrier mobility of the Holstein-Peierls model using the matrix product states approach, we locate the phonon-assisted transport, transient localization and band-like regimes as a function of the transfer integral (V) and the nonlocal electron-phonon couplings (ΔV), and their distinct transport behaviors are analyzed by carrier mobility, mean free path, optical conductivity and one-particle spectral function. We also identify an "intermediate regime" where none of the established pictures applies, and the generally perceived hopping regime is found to be at a very limited end in the proposed regime paradigm.

A general automatic method for optimal construction of matrix product operators using bipartite graph theory.

Constructing matrix product operators (MPOs) is at the core of the modern density matrix renormalization group (DMRG) and its time dependent formulation. For the DMRG to be conveniently used in different problems described by different Hamiltonians, in this work, we propose a new generic algorithm to construct the MPO of an arbitrary operator with a sum-of-products form based on the bipartite graph theory. We show that the method has the following advantages: (i) it is automatic in that only the definition of the operator is required; (ii) it is symbolic thus free of any numerical error; (iii) the complementary operator technique can be fully employed so that the resulting MPO is globally optimal for any given order of degrees of freedom; and (iv) the symmetry of the system could be fully employed to reduce the dimension of MPO. To demonstrate the effectiveness of the new algorithm, the MPOs of Hamiltonians ranging from the prototypical spin-boson model and the Holstein model to the more complicated ab initio electronic Hamiltonian and the anharmonic vibrational Hamiltonian with the sextic force field are constructed. It is found that for the former three cases, our automatic algorithm can reproduce exactly the same MPOs as the optimally hand-crafted ones already known in the literature.

Applying Marcus theory to describe the carrier transports in organic semiconductors: Limitations and beyond.

Marcus theory has been successfully applied to molecular design for organic semiconductors with the aid of quantum chemistry calculations for the molecular parameters: the intermolecular electronic coupling V and the intramolecular charge reorganization energy λ. The assumption behind this is the localized nature of the electronic state for representing the charge carriers, being holes or electrons. As far as the quantitative description of carrier mobility is concerned, the direct application of Marcus semiclassical theory usually led to underestimation of the experimental data. A number of effects going beyond such a semiclassical description will be introduced here, including the quantum nuclear effect, dynamic disorder, and delocalization effects. The recently developed quantum dynamics simulation at the time-dependent density matrix renormalization group theory is briefly discussed. The latter was shown to be a quickly emerging efficient quantum dynamics method for the complex system.

Finite-Temperature TD-DMRG for the Carrier Mobility of Organic Semiconductors.

A large number of nonadiabatic dynamical studies have been applied to reveal the nature of carrier transport in organic semiconductors with different approximations. We present here a "nearly exact" graphical-process-unit-based finite-temperature time-dependent density matrix renormalization group (TD-DMRG) method to evaluate the carrier mobility in organic semiconductors, as described by the electron-phonon model, in particular, in rubrene crystal, one of the prototypical organic semiconductors, with parameters derived from first-principles. We find that (i) TD-DMRG is a general and robust method that can bridge the gap between hopping and band pictures, covering a wide range of electronic coupling strengths and (ii) with realistic parameters, TD-DMRG is able to account for the experimentally observed "band-like" transport behavior (∂µ/∂T < 0) in rubrene. We further study the long-standing puzzle of the isotope effect for charge transport and unambiguously demonstrate that the negative isotope effect (∂µ/∂m < 0 where m is the atomic mass) should be universal.

Finite Temperature Dynamical Density Matrix Renormalization Group for Spectroscopy in Frequency Domain.

We present a novel methodology through casting the dynamical density matrix renormalization group (DDMRG) into the matrix product state (MPS) formulation to calculate the spectroscopy at finite temperature for molecular aggregates. The frequency domain algorithm can avoid the time evolution accumulation of error and is naturally suitable for parallelization, in addition to facile graphic processing unit (GPU) acceleration. The high accuracy is demonstrated by simulating the optical spectra of vibronic model systems ranging from an exactly solvable dimer model to a more complex real-world perylene bisimide (PBI) J-aggregate. The relationship between the 0-0 emission strength and the exciton thermal coherent length is discussed for linearly stacked aggregates. The computing performance largely boosted by GPU demonstrates that DDMRG emerges as a promising method to study dynamical properties for complex systems.

Numerical assessment for accuracy and GPU acceleration of TD-DMRG time evolution schemes.

The time dependent density matrix renormalization group (TD-DMRG) has become one of the cutting edge methods of quantum dynamics for complex systems. In this paper, we comparatively study the accuracy of three time evolution schemes in the TD-DMRG, the global propagation and compression method with the Runge-Kutta algorithm (P&C-RK), the time dependent variational principle based methods with the matrix unfolding algorithm (TDVP-MU), and with the projector-splitting algorithm (TDVP-PS), by performing benchmarks on the exciton dynamics of the Fenna-Matthews-Olson complex. We show that TDVP-MU and TDVP-PS yield the same result when the time step size is converged and they are more accurate than P&C-RK4, while TDVP-PS tolerates a larger time step size than TDVP-MU. We further adopt the graphical processing units to accelerate the heavy tensor contractions in the TD-DMRG, and it is able to speed up the TDVP-MU and TDVP-PS schemes by up to 73 times.

Understanding Carrier Transport in Organic Semiconductors: Computation of Charge Mobility Considering Quantum Nuclear Tunneling and Delocalization Effects.

Despite great attention to the charge transport in organic semiconductors (OSC) over the last decades, the underlying mechanism is still controversial. After our theoretical position in 2009, the quantum nuclear tunneling effect has been proven by more and more experiments to play an essential role for charge transport in organic and polymeric materials. On the other hand, back in the 1970s, it was proposed that the nature of charge transport could be analyzed by the isotope effect, which, however, has not been confirmed either experimentally or theoretically. In this Perspective, we review the understanding of microscopic mechanisms on charge transport by using different transport mechanisms from hopping to band transport. Particularly, we point out that the isotope effect, which is absent in the semiclassical Marcus theory, should be negative for the localized charge transport with quantum nuclear tunneling. We conclude that the quantum nuclear tunneling effect dominates the charge transport in OSCs.

Risk factors for complications following immediate tissue expander based breast reconstruction in Taiwanese population.

BACKGROUND/PURPOSE: Breast cancer patients in Asia show considerable disparities from Caucasian patients, such as younger age of onset and lower rates of smoking, obesity, and diabetes. Findings of prior studies regarding risk factors associated with complications in tissue expander may not hold for Asian populations, since most of these studies involved Caucasian patients. In this study, we surveyed risk factors in the Taiwanese population, providing additional evidence about the important differences and discuss the implications for clinical practice. METHODS: Patients who underwent immediate, two-stage, tissue expander breast reconstruction from December 2008 to August 2014 in the National Taiwan University Hospital, Taipei, Taiwan were included. Follow-up observations of all patients were conducted until December 2014. Complications occurring during the tissue expander stage were evaluated. Multivariate regression modeling was used to identify risk factors for complications. RESULTS: A total of 246 consecutive, immediate, smooth round tissue expander placements were performed for breast reconstruction. The most common complication was skin necrosis (4.9%), followed by wound dehiscence (4.1%). In the multivariate model, body mass index (BMI) ≥ 24 kg/m2 was the only risk factor that reached statistical significance (odds ratio: 2.41, 95% confidence interval: 1.17-4.96). CONCLUSION: We provided evidence that racial disparities have an impact on the risk factors for complications associated with tissue expander breast reconstruction. BMI≥24 kg/m2 was the only risk factor significantly associated with complications. Clinically, BMI≥24 kg/m2, rather than the standard definition of obesity (BMI > 30 kg/m2), may be a more suitable cutoff point for risk in patients of Asian ethnicity.

Treatment and long-term follow-up of oral cancer postoperative sialorrhea with dermal sling operation.

INTRODUCTION: Reconstruction of a full-thickness defect that includes oral commissure presents a considerable challenge to maxillofacial and plastic surgeons. The goals of reconstruction are both functional and cosmetic. Sialorrhea, or drooling, is a major problem after flap reconstruction and influences the quality of life of the patient. In this article, we report on our experience performing a dermal sling operation to treat postoperative sialorrhea in patients with oral cancer. MATERIALS AND METHODS: Preoperative and postoperative levels of sialorrhea were evaluated based on the Drooling Severity and Frequency Scale. Dermal sling operations were performed on 27 patients from January 2000 to December 2013. In these patients, 12 cases were reviewed and followed up over 1 year. RESULTS: Of the 12 patients, 11 were men and one was a woman, with the mean age of 58 years (range, 40-79 years). There were no operative complications. The mean preoperative score was 4.75 (range, 3-7), and the mean postoperative score was 3.83 (range, 2-5). This change was significant (P=0.005), with valuation with the Wilcoxon signed rank test. The mean time of follow-up was 3.5 years (range, 1.1-7.7 years). CONCLUSIONS: The dermal sling operation is an acceptable treatment for postoperative sialorrhea in patients with commissure-involved oral cancer.

Regional identification, partition, and integral phase unwrapping method for processing moiré interferometry images.

We present a new method of regional identification, partition, and integral (RIPI) phase unwrapping for processing images, especially those with low quality, obtained from moiré interferometry experiments. By introducing the principle of preorder traversal of a general tree in data structures and then by applying the idea of a regional integral, the proposed method makes regional partition and phase evaluation much easier and more accurate, and it also overcomes the common faults that can occur when conventional approaches, such as line defects, are used. Examples are given to demonstrate the advantage and applicability of the proposed RIPI method when processing experimental images. It is shown that the proposed method works well for global phase distribution, and, at the same time, local mutational information is preserved and limited to its vicinity without affecting other parts.