Born–Oppenheimer and Non-Born–Oppenheimer, Atomic and Molecular Calculations with Explicitly Correlated Gaussians

Efficient Exploitation of Numerical Quadrature with Distance-Dependent Integral Screening in Explicitly Correlated F12 Theory: Linear Scaling Evaluation of the Most Expensive RI-MP2-F12 Term

We present a linear scaling atomic orbital based algorithm for the computation of the most expensive exchange-type RI-MP2-F12 term by employing numerical quadrature in combination with CABS-RI to avoid six-center-three-electron integrals. Furthermore, a robust distance-dependent integral screening scheme, based on integral partition bounds [Thompson, T. H.; Ochsenfeld, C. J. Chem. Phys. 2019, 150, 044101], is used to drastically reduce the number of the required three-center-one-electron integrals substantially. The accuracy of our numerical quadrature/CABS-RI approach and the corresponding integral screening is thoroughly assessed for interaction and isomerization energies across a variety of numerical integration grids. Our method outperforms the standard density fitting/CABS-RI approach with errors below 1 μEh even for small grid sizes and moderate screening thresholds. The choice of the grid size and screening threshold allows us to tailor our ansatz to a desired accuracy and computational efficiency. We showcase the approach's effectiveness for the chemically relevant system valinomycin, employing a triple-ζ F12 basis set combination (C54H90N6O18, 5757 AO basis functions, 10,266 CABS basis functions, 735,783 grid points). In this context, our ansatz achieves higher accuracy combined with a 135× speedup compared to the classical density fitting based variant, requiring notably less computation time than the corresponding RI-MP2 calculation. Additionally, we demonstrate near-linear scaling through calculations on linear alkanes. We achieved an 817-fold acceleration for C80H162 and an extrapolated 28,765-fold acceleration for C200H402, resulting in a substantially reduced computational time for the latter─from 229 days to just 11.5 min. Our ansatz may also be adapted to the remaining MP2-F12 terms, which will be the subject of future work.

Gaussians for Electronic and Rovibrational Quantum Dynamics

The assumptions underpinning the adiabatic Born-Oppenheimer (BO) approximation are broken for molecules interacting with attosecond laser pulses, which generate complicated coupled electronic-nuclear wave packets that generally will have components of electronic and dissociation continua as well as bound-state contributions. The conceptually most straightforward way to overcome this challenge is to treat the electronic and nuclear degrees of freedom on equal quantum-mechanical footing by not invoking the BO approximation at all. Explicitly correlated Gaussian (ECG) basis functions have proved successful for non-BO calculations of stationary molecular states and energies, reproducing rovibrational absorption spectra with very high accuracy. In this Article, we present a proof-of-principle study of the ability of fully flexible ECGs (FFECGs) to capture the intricate electronic and rovibrational dynamics generated by short, high-intensity laser pulses. By fitting linear combinations of FFECGs to accurate wave function histories obtained on a large real-space grid for a regularized 2D model of the hydrogen atom and for the 2D Morse potential, we demonstrate that FFECGs provide a very compact description of laser-driven electronic and rovibrational dynamics.

Quantum Definition of Molecular Structure

Molecular structure, a key concept of chemistry, has remained elusive from the perspective of all-particle quantum mechanics, despite many efforts. Viewing molecular structure as a manifestation of strong statistical correlation between nuclear positions, we propose a practical method based on Markov chain Monte Carlo sampling and unsupervised machine learning. Application to the D3+ molecule unambiguously shows that it possesses an equilateral triangular structure. These results provide a major step forward in our understanding of the molecular structure from fundamental quantum principles.

Explicitly correlated Gaussians for high-precision variational calculations of S^{e}, P^{e}, and D^{e} states of quantum systems: An efficient algorithm

In this work we consider an efficient algorithm for variational calculations of quantum few-particle systems in S, P, and D states of the even parity using all-particle explicitly correlated Gaussian (ECG) basis sets. We primarily focus on the description of states where the dominant configuration contains either two particles in p states or a single particle in a d state (all other particles are in s states). The basis functions we consider are products of spherically symmetric ECGs and bipolar harmonics. We introduced a scheme for deriving expressions for matrix elements of the overlap, kinetic and potential energy, as well as their derivatives with respect to the nonlinear parameters of the Gaussians. This allowed us to improve the efficiency of numerical calculations of the matrix elements (which is the most critical part of any code that uses ECGs) by one to two orders of magnitude compared to previous implementations. We provide a complete set of formulas for all basic matrix elements using the formalism of matrix differential calculus and discuss some technical details relevant to their efficient implementation. Lastly, we report a few example calculations of the ^{2}D^{e} states of the Li atom, ^{4}P^{e} and ^{2}D^{e} states of the B atom, and ^{3}P^{e} states of the C atom.

Second-Order Self-Consistent Field Algorithms: From Classical to Quantum Nuclei

This work presents a general framework for deriving exact and approximate Newton self-consistent field (SCF) orbital optimization algorithms by leveraging concepts borrowed from differential geometry. Within this framework, we extend the augmented Roothaan-Hall (ARH) algorithm to unrestricted electronic and nuclear-electronic calculations. We demonstrate that ARH yields an excellent compromise between stability and computational cost for SCF problems that are hard to converge with conventional first-order optimization strategies. In the electronic case, we show that ARH overcomes the slow convergence of orbitals in strongly correlated molecules with the example of several iron-sulfur clusters. For nuclear-electronic calculations, ARH significantly enhances the convergence already for small molecules, as demonstrated for a series of protonated water clusters.

Quantum algorithms for quantum dynamics

Among the many computational challenges faced across different disciplines, quantum-mechanical systems pose some of the hardest ones and offer a natural playground for the growing field of quantum technologies. In this Perspective, we discuss quantum algorithmic solutions for quantum dynamics, reporting on the latest developments and offering a viewpoint on their potential and current limitations. We present some of the most promising areas of application and identify possible research directions for the coming years.

Achieving vibrational energies of diatomic systems with high quality by machine learning improved DFT method

When using ab initio methods to obtain high-quality quantum behavior of molecules, it often involves a lot of trial-and-error work in algorithm design and parameter selection, which requires enormous time and computational resource costs. In the study of vibrational energies of diatomic molecules, we found that starting from a low-precision DFT model and then correcting the errors using the high-dimensional function modeling capabilities of machine learning, one can considerably reduce the computational burden and improve the prediction accuracy. Data-driven machine learning is able to capture subtle physical information that is missing from DFT approaches. The results of ¹²C¹⁶O, ²⁴MgO and Na³⁵Cl show that, compared with CCSD(T)/cc-pV5Z calculation, this work improves the prediction accuracy by more than one order of magnitude, and reduces the computation cost by more than one order of magnitude.

Quantum Proton Effects from Density Matrix Renormalization Group Calculations

We recently introduced [J. Chem. Phys. 2020, 152, 204103] the nuclear-electronic all-particle density matrix renormalization group (NEAP-DMRG) method to solve the molecular Schrödinger equation, based on a stochastically optimized orbital basis, without invoking the Born-Oppenheimer approximation. In this work, we combine the DMRG method with the nuclear-electronic Hartree-Fock (NEHF-DMRG) approach, treating nuclei and electrons on the same footing. Inter- and intraspecies correlations are described within the DMRG method without truncating the excitation degree of the full configuration interaction wave function. We extend the concept of orbital entanglement and mutual information to nuclear-electronic wave functions and demonstrate that they are reliable metrics to detect strong correlation effects. We apply the NEHF-DMRG method to the HeHHe⁺ molecular ion, to obtain accurate proton densities, ground-state total energies, and vibrational transition frequencies by comparison with state-of-the-art data obtained with grid-based approaches and modern configuration interaction methods. For HCN, we improve on the accuracy of the latter approaches with respect to both the ground-state absolute energy and proton density, which is a major challenge for multireference nuclear-electronic state-of-the-art methods.

Deformed explicitly correlated Gaussians

Deformed explicitly correlated Gaussian (DECG) basis functions are introduced, and their matrix elements are calculated. All matrix elements can be calculated analytically in a closed form, except the Coulomb one, which has to be approximated by a Gaussian expansion. The DECG basis functions can be used to solve problems with nonspherical potentials. One example of such potential is the dipole self-interaction term in the Pauli-Fierz Hamiltonian. Examples are presented showing the accuracy and necessity of deformed Gaussian basis functions to accurately solve light-matter coupled systems in cavity QED.

Orientational decoherence within molecules and emergence of the molecular shape

The question of classicality is addressed in relation with the shape of the nuclear skeleton of molecular systems. As the most natural environment, the electrons of the molecule are considered as continuously monitoring agents for the nuclei. For this picture, an elementary formalism of decoherence theory is developed and numerical results are presented for few-particle systems. The numerical examples suggest that the electron-nucleus Coulomb interaction is sufficient for inducing a blurred shape with strong quantum coherences in compounds of the lightest elements, H₂, D₂, T₂, and HeH⁺.

Nonadiabatic, Relativistic, and Leading-Order QED Corrections for Rovibrational Intervals of ^{4}He_{2}^{+} (X ^{2}Σ_{u}^{+})

The rovibrational intervals of the ^{4}He_{2}^{+} molecular ion in its X ^{2}Σ_{u}^{+} ground electronic state are computed by including the nonadiabatic, relativistic, and leading-order quantum-electrodynamics corrections. Good agreement of theory and experiment is observed for the rotational excitation series of the vibrational ground state and the fundamental vibration. The lowest-energy rotational interval is computed to be 70.937 69(10)  cm^{-1} in agreement with the most recently reported experimental value, 70.937 589(23)(60)_{sys}  cm^{-1} [L. Semeria et al., Phys. Rev. Lett. 124, 213001 (2020)PRLTAO0031-900710.1103/PhysRevLett.124.213001].

Hardware efficient quantum algorithms for vibrational structure calculations

We introduce a framework for the calculation of ground and excited state energies of bosonic systems suitable for near-term quantum devices and apply it to molecular vibrational anharmonic Hamiltonians. Our method supports generic reference modal bases and Hamiltonian representations, including the ones that are routinely used in classical vibrational structure calculations. We test different parametrizations of the vibrational wavefunction, which can be encoded in quantum hardware, based either on heuristic circuits or on the bosonic Unitary Coupled Cluster Ansatz. In particular, we define a novel compact heuristic circuit and demonstrate that it provides a good compromise in terms of circuit depth, optimization costs, and accuracy. We evaluate the requirements, number of qubits and circuit depth, for the calculation of vibrational energies on quantum hardware and compare them with state-of-the-art classical vibrational structure algorithms for molecules with up to seven atoms.

Development of nuclear basis sets for multicomponent quantum chemistry methods

The nuclear-electronic orbital (NEO) framework provides a practical approach for directly incorporating nuclear quantum effects and non-Born-Oppenheimer effects of specified nuclei, typically protons, into quantum chemistry calculations. Multicomponent wave function based methods, such as NEO coupled cluster singles and doubles, and multicomponent density functional theory (DFT), such as NEO-DFT, require the appropriate selection of electronic and nuclear basis sets. Although a wide array of electronic basis sets are available, systematically developed nuclear basis sets that balance accuracy and efficiency have been lacking. Herein, a series of nuclear basis sets are developed and shown to be accurate and efficient for describing both ground and excited state properties of multicomponent systems in which electrons and specified protons are treated quantum mechanically. Three series of Gaussian-type nuclear basis sets, denoted PB4, PB5, and PB6, are developed with varying levels of angular momentum. A machine-learning optimization procedure relying on the Gaussian process regression method is utilized to accelerate the optimization process. The basis sets are validated in terms of predictions of ground state energies, proton densities, proton affinities, and proton vibrational excitation energies, allowing the user to select the desired balance between accuracy and efficiency for the properties of interest. These nuclear basis sets will enhance the tractability of NEO methods for applications to a wide range of chemical systems.

Computer program ATOM-MOL-nonBO for performing calculations of ground and excited states of atoms and molecules without assuming the Born-Oppenheimer approximation using all-particle complex explicitly correlated Gaussian functions

In this work, we describe a computer program called ATOM-MOL-nonBO for performing bound state calculations of small atoms and molecules without assuming the Born-Oppenheimer approximation. All particles forming the systems, electrons and nuclei, are treated on equal footing. The wave functions of the bound states are expanded in terms of all-particle one-center complex explicitly correlated Gaussian functions multiplied by Cartesian angular factors. As these Gaussian functions are eigenfunctions of the operator representing the square of the total angular momentum of the system, the problem separates and calculations of states corresponding to different values of the total rotational quantum number can be solved independently from each other. Due to thorough variational optimization of the Gaussian exponential parameters, the method allows us to generate very accurate wave functions. The optimization is aided by analytically calculated energy gradient determined with respect to the parameters. Three examples of calculations performed for diatomic and triatomic molecules are shown as an illustration of calculations that can be performed with this program. Finally, we discuss the limitations, applicability range, and bottlenecks of the program.

Nuclear-electronic all-particle density matrix renormalization group

We introduce the Nuclear-Electronic All-Particle Density Matrix Renormalization Group (NEAP-DMRG) method for solving the time-independent Schrödinger equation simultaneously for electrons and other quantum species. In contrast to the already existing multicomponent approaches, in this work, we construct from the outset a multi-reference trial wave function with stochastically optimized non-orthogonal Gaussian orbitals. By iterative refining of the Gaussians' positions and widths, we obtain a compact multi-reference expansion for the multicomponent wave function. We extend the DMRG algorithm to multicomponent wave functions to take into account inter- and intra-species correlation effects. The efficient parameterization of the total wave function as a matrix product state allows NEAP-DMRG to accurately approximate the full configuration interaction energies of molecular systems with more than three nuclei and 12 particles in total, which is currently a major challenge for other multicomponent approaches. We present the NEAP-DMRG results for two few-body systems, i.e., H₂ and H₃ ⁺, and one larger system, namely, BH₃.

Simulation of a hydrogen atom in a laser field using the time-dependent variational principle

The time-dependent variational principle is used to optimize the linear and nonlinear parameters of Gaussian basis functions to solve the time-dependent Schrödinger equation in one and three dimensions for a one-body soft Coulomb potential in a laser field. The accuracy is tested comparing the solution to finite difference grid calculations using several examples. The approach is not limited to one particle systems and the example presented for two electrons demonstrates the potential to tackle larger systems using correlated basis functions.

Molecular Vibrational Frequencies with Multiple Quantum Protons within the Nuclear-Electronic Orbital Framework

The nuclear-electronic orbital (NEO) approach treats all electrons and specified nuclei, typically protons, on the same quantum mechanical level. Proton vibrational excitations can be calculated using multicomponent time-dependent density functional theory (NEO-TDDFT) for fixed classical nuclei. Recently the NEO-DFT(V) approach was developed to enable the calculation of molecular vibrational frequencies for modes composed of both classical and quantum nuclei. This approach uses input from NEO-TDDFT to construct an extended NEO Hessian that depends on the expectation values of the quantum protons as well as the classical nuclear coordinates. Herein strategies are devised for extending these approaches to molecules with multiple quantum protons in a self-contained, effective, and computationally practical manner. The NEO-TDDFT method is shown to describe vibrational excitations corresponding to collective nuclear motions, such as linear combinations of proton vibrational excitations. The NEO-DFT(V) approach is shown to incorporate the most significant anharmonic effects in the molecular vibrations, particularly for the hydrogen stretching modes. These theoretical strategies pave the way for a wide range of multicomponent quantum chemistry applications.

H3 + as a five-body problem described with explicitly correlated Gaussian basis sets

Various explicitly correlated Gaussian (ECG) basis sets are considered for the solution of the molecular Schrödinger equation with particular attention to the simplest polyatomic system, H₃ ⁺. Shortcomings and advantages are discussed for plain ECGs, ECGs with the global vector representation, floating ECGs and their numerical projection, and ECGs with complex parameters. The discussion is accompanied with particle density plots to visualize the observations. In order to be able to use large complex ECG basis sets in molecular calculations, a numerically stable algorithm is developed, the efficiency of which is demonstrated for the lowest rotationally and vibrationally excited states of H₂ and H₃ ⁺.

The any particle molecular orbital/molecular mechanics approach

A computational scheme is proposed to broaden the range of applications of multicomponent methodologies for the study of local properties of big molecular systems existing in the gas phase and in solvated environments. This scheme extends the any particle molecular orbital (APMO) approach in the quantum mechanics/molecular mechanics (QM/MM) framework. As a first assessment of the performance of the proposed approach, we estimate the proton affinities (PAs) of seventy amines in the gas phase and the proton binding energies (PBEs) in the gas phase and in an explicitly solvated environment of the sixty-one protons present in the chignolin protein. These calculations are performed with the QM/MM versions of the APMO second-order proton propagator (APMO-PP2) and the APMO extended Koopmans' theorem (APMO-KT) approaches. Calculated PAs and PBEs show significant reductions in the computational effort with a reduced loss in accuracy. These results suggest that the APMO/MM scheme might be used as a low-cost multi-component alternative for studies of local properties in big molecular systems. Graphical Abstract QMMM regions and CPU times for the APMO/MM approach.

Quantum Chemistry in the Age of Quantum Computing

Practical challenges in simulating quantum systems on classical computers have been widely recognized in the quantum physics and quantum chemistry communities over the past century. Although many approximation methods have been introduced, the complexity of quantum mechanics remains hard to appease. The advent of quantum computation brings new pathways to navigate this challenging and complex landscape. By manipulating quantum states of matter and taking advantage of their unique features such as superposition and entanglement, quantum computers promise to efficiently deliver accurate results for many important problems in quantum chemistry, such as the electronic structure of molecules. In the past two decades, significant advances have been made in developing algorithms and physical hardware for quantum computing, heralding a revolution in simulation of quantum systems. This Review provides an overview of the algorithms and results that are relevant for quantum chemistry. The intended audience is both quantum chemists who seek to learn more about quantum computing and quantum computing researchers who would like to explore applications in quantum chemistry.

Multicomponent coupled cluster singles and doubles and Brueckner doubles methods: Proton densities and energies

The nuclear-electronic orbital (NEO) framework enables computationally practical coupled cluster calculations of multicomponent molecular systems, in which all electrons and specified nuclei, typically protons, are treated quantum mechanically. In addition to energies, computing accurate proton densities is essential for the calculation of reliable molecular properties, including vibrationally averaged geometries and vibrational frequencies. Herein, the Lagrangian formalism for the multicomponent coupled cluster with single and double excitations (NEO-CCSD) method is derived and implemented. The multicomponent coupled cluster with double excitations method using optimized Brueckner orbitals, denoted as NEO-BCCD, is also developed. Both of these methods are used to compute the proton densities for two molecular systems. The results illustrate that orbital relaxation effects, which can be included either indirectly with the NEO-CCSD method or directly with the NEO-BCCD method, are critical for computing even qualitatively accurate proton densities. Both methods are also able to provide accurate proton affinities and vibrationally averaged optimized geometries. This Lagrangian formalism will enable the calculation of other properties such as analytical nuclear gradients and Hessians with NEO coupled cluster methods. Moreover, the accuracy of these methods may be improved systematically by the inclusion of higher-order excitations. Thus, this work provides the foundation for a wide range of future methodological developments and applications within the NEO framework.

Effective non-adiabatic Hamiltonians for the quantum nuclear motion over coupled electronic states

The quantum mechanical motion of the atomic nuclei is considered over a single- or a multidimensional subspace of electronic states which is separated by a gap from the rest of the electronic spectrum over the relevant range of nuclear configurations. The electron-nucleus Hamiltonian is block-diagonalized up to O(εⁿ⁺¹) through a unitary transformation of the electronic subspace, and the corresponding nth-order effective Hamiltonian is derived for the quantum nuclear motion. Explicit but general formulas are given for the second- and the third-order corrections. As a special case, the second-order Hamiltonian corresponding to an isolated electronic state is recovered which contains the coordinate-dependent mass-correction terms in the nuclear kinetic energy operator. For a multidimensional, explicitly coupled electronic band, the second-order Hamiltonian contains the usual Born-Oppenheimer terms and nonadiabatic corrections, but generalized mass-correction terms appear as well. These, earlier neglected terms, perturbatively account for the outlying (discrete and continuous) electronic states not included in the explicitly coupled electronic subspace.

Molecular Vibrational Frequencies within the Nuclear-Electronic Orbital Framework

A significant challenge for multicomponent quantum chemistry methods is the calculation of vibrational frequencies for comparison to experiment. The nuclear-electronic orbital (NEO) approach treats specified nuclei, typically key protons, quantum mechanically. The Born-Oppenheimer separation between the quantum and classical nuclei prevents the direct calculation of vibrational frequencies corresponding to modes composed of both types of nuclei. Herein an effective strategy for calculating the vibrational frequencies of the entire molecule within the NEO framework is devised and implemented. This strategy requires diagonalization of an extended NEO Hessian that depends on the expectation values of the quantum nuclei as well as the coordinates of the classical nuclei and is constructed with input from multicomponent time-dependent density functional theory (NEO-TDDFT). Application of this NEO-DFT(V) approach to molecular systems illustrates that it accurately incorporates the most significant anharmonic effects. This general theoretical formulation opens up a broad spectrum of new directions for multicomponent quantum chemistry.

Side chain torsion dictates planarity and ionizability of green fluorescent protein's chromophore leading to spectral perturbations

Spectral characteristics of fluorescent proteins (FPs) are well studied, and through protein engineering, several FP variants constituting entire visible spectrum have been created. One of the most common mechanisms attributed to spectral shifts in FP is excited state proton transfer (ESPT), hydroxyl moiety protonation and deprotonation, along with chromophore cis-trans isomerism. The most widely studied FPs are those derived from avGFP (Aequorea victoria GFP) and Dsred (Discosoma coral). Apart from the above mechanism, certain interacting residues are said to play a vital role in altering the proton transfer pathway leading to numerous spectral variants. Similarly, the hydrogen-bonded networks solely cannot dictate the energy landscape of FPs. Non-bonded interactions also can create secondary harmonic shifts by dipole-dipole inductions. Side chain contacts tend to alter the topological and torsional geometry, thereby disturbing the chromophore's planarity. Side chain torsional variations have almost been unaccounted for their distortions in FPs. We hypothesize the torsional landscape and altered residual interactions as prominent factors for the spectral shifts. Through our 200 ns molecular dynamics investigation, we prospect that van der Waals packing in Dsred is more compact than that of avGFP, thus creating a low solvent occupiable environment and reduced solvent interactions having higher red spectral shift. The torsional changes of wild avGFP, S65T avGFP and Dsred have been studied to comprehend the inter-residual contact distance and the geometrical descriptors. Communicated by Ramaswamy H. Sarma.

Explicitly correlated Gaussian functions with shifted-center and projection techniques in pre-Born-Oppenheimer calculations

Numerical projection methods are elaborated for the calculation of eigenstates of the non-relativistic many-particle Coulomb Hamiltonian with selected rotational and parity quantum numbers employing shifted explicitly correlated Gaussian functions, which are, in general, not eigenfunctions of the total angular momentum and parity operators. The increased computational cost of numerically projecting the basis functions onto the irreducible representations of the three dimensional rotation-inversion group is the price to pay for the increased flexibility of the basis functions. This increased flexibility allowed us to achieve a substantial improvement for the variational upper bound to the Pauli-allowed ground-state energy of the H 3 + = { p + , p + , p + , e - , e - } molecular ion treated as an explicit five-particle system. We compare our pre-Born-Oppenheimer result obtained for this molecular ion with rotational-vibrational calculations carried out on a potential energy surface.

Solving the Schrödinger equation of atoms and molecules: Chemical-formula theory, free-complement chemical-formula theory, and intermediate variational theory

Chemistry is governed by the principle of quantum mechanics as expressed by the Schrödinger equation (SE) and Dirac equation (DE). The exact general theory for solving these fundamental equations is therefore a key for formulating accurately predictive theory in chemical science. The free-complement (FC) theory for solving the SE of atoms and molecules proposed by one of the authors is such a general theory. On the other hand, the working theory most widely used in chemistry is the chemical formula that refers to the molecular structural formula and chemical reaction formula, collectively. There, the central concepts are the local atomic concept, transferability, and from-atoms-to-molecule concept. Since the chemical formula is the most successful working theory in chemistry ever existed, we formulate our FC theory to have the structure reflecting the chemical formula. Our basic postulate is that as far as the SE is the principle of chemistry, its solutions for chemistry should have the structure that can be related to the chemical formulas. So, in this paper, we first formulate a theory that designs the wave function to reflect the structure of the chemical formula. We call this theory chemical formula theory (CFT). In the CFT, we place the valence ground and excited states of each atom at each position of the chemical formula of the molecule and let them interact using their free valences to form the ground and excited states of the molecule. The principle there is the variational principle so that the ground and excited states obtained satisfy the orthogonality and Hamiltonian-orthogonality relations. Then, we formulate the exact FC theory starting from the initial functions produced by the CFT. This FC theory is referred to as free-complement chemical-formula theory (FC-CFT), which is expected to describe efficiently the solution of the SE by the above reason. The FC-CFT wave function is modified from that of CFT. Since this modification is done by the exact SE, its analysis may give some insights to chemists that assist their chemistry. Thus, this theory would be not only exact but also conceptually useful. Furthermore, the intermediate theory between CFT and FC-CFT would also be useful. There, we use only integratable functions and apply the variational principle so that we refer to this theory as FC-CFT-variational (FC-CFT-V). It is an advanced theory of CFT. Since the variational method is straightforward and powerful, we can do extensive chemical studies in a reasonable accuracy. After finishing such studies, if we still need an exact level of solutions, we add the remaining functions of the FC-CFT and perform the exact calculations. Furthermore, when we deal with large and even giant molecules, the inter-exchange (iExg) theory for the antisymmetry rule introduced previously leads to a large simplification. There, the inter-exchanges between distant electron pairs fade away so that only Coulombic interactions survive. Further in giant systems, even an electrostatic description becomes possible. Then, the FC-CFT for exactly solving the SE would behave essentially to order N for large and giant molecular systems, though the pre-factor should be very large and must be minimized.

T-DM1 and brain metastases: Clinical outcome in HER2-positive metastatic breast cancer

BACKGROUND

We reported the results of an Italian large retrospective analysis that evaluated the effectiveness and safety of T-DM1 in 'field-practice' breast cancer patients. We performed a sub-analysis to investigate the clinical activity of T-DM1 in patients with brain metastases (BMs).

METHODS

The records of 87 adult women with HER2-positive breast cancer and BMs treated with T-DM1 were reviewed. Their clinical outcomes were compared with those of 216 patients without central nervous system (CNS) involvement.

RESULTS

Response to T-DM1 treatment in BMs was available for 53 patients in the BM group (60.9%): two patients reported a complete response (3.8%), 11 patients obtained partial response (20.7%; overall response rate: 24.5%), 16 patients had a stable disease (30.1%). Regarding extracranial disease, a total of 77 and 191 patients were evaluable for response in BM group and non-BM group, respectively. The overall response rate was 35.1% in the BM group and 38.3% in the non-BM group; disease control rate was 53.3% and 66.6%, respectively. At a median follow-up of 16 months (range: 1-55), median cumulative progression-free survival (PFS) was 7 months (95% CI: 5.4-8.6) in the BM group and 8 months (95% CI: 5.7-10.3) in the non-BM group. In the second-line setting, PFS was 5 (95% CI: 3.1-6.9) versus 11 (95% CI: 7.1-14.9) months (p = 0.01). Overall survival was 14 months (95% CI: 12.2-15.8) in the BM group and 32 months (95% CI: 24.4-39.6) in the non-BM group (p < 0.0001).

CONCLUSIONS

T-DM1 is active in breast cancer patients with BMs.

Generalized elimination of the global translation from explicitly correlated Gaussian functions

This paper presents the multi-channel generalization of the center-of-mass kinetic energy elimination approach [B. Simmen et al., Mol. Phys. 111, 2086 (2013)] when the Schrödinger equation is solved variationally with explicitly correlated Gaussian functions. The approach has immediate relevance in many-particle systems which are handled without the Born-Oppenheimer approximation and can be employed also for Dirac-type Hamiltonians. The practical realization and numerical properties of solving the Schrödinger equation in laboratory-frame Cartesian coordinates are demonstrated for the ground rovibronic state of the H₂⁺={p⁺,p⁺,e^-} ion and the H₂ = {p⁺, p⁺, e^-, e^-} molecule.

Toward a muon-specific electronic structure theory: effective electronic Hartree-Fock equations for muonic molecules

An effective set of Hartree-Fock (HF) equations are derived for electrons of muonic systems, i.e., molecules containing a positively charged muon, conceiving the muon as a quantum oscillator, which are completely equivalent to the usual two-component HF equations used to derive stationary states of the muonic molecules. In these effective equations, a non-Coulombic potential is added to the orthodox coulomb and exchange potential energy terms, which describes the interaction of the muon and the electrons effectively and is optimized during the self-consistent field cycles. While in the two-component HF equations a muon is treated as a quantum particle, in the effective HF equations it is absorbed into the effective potential and practically transformed into an effective potential field experienced by electrons. The explicit form of the effective potential depends on the nature of muon's vibrations and is derivable from the basis set used to expand the muonic spatial orbital. The resulting effective Hartree-Fock equations are implemented computationally and used successfully, as a proof of concept, in a series of muonic molecules containing all atoms from the second and third rows of the Periodic Table. To solve the algebraic version of the equations muon-specific Gaussian basis sets are designed for both muon and surrounding electrons and it is demonstrated that the optimized exponents are quite distinct from those derived for the hydrogen isotopes. The developed effective HF theory is quite general and in principle can be used for any muonic system while it is the starting point for a general effective electronic structure theory that incorporates various types of quantum correlations into the muonic systems beyond the HF equations.

On the physical interpretation of the nuclear molecular orbital energy

Recently, several groups have extended and implemented molecular orbital (MO) schemes to simultaneously obtain wave functions for electrons and selected nuclei. Many of these schemes employ an extended Hartree-Fock approach as a first step to find approximate electron-nuclear wave functions and energies. Numerous studies conducted with these extended MO methodologies have explored various effects of quantum nuclei on physical and chemical properties. However, to the best of our knowledge no physical interpretation has been assigned to the nuclear molecular orbital energy (NMOE) resulting after solving extended Hartree-Fock equations. This study confirms that the NMOE is directly related to the molecular electrostatic potential at the position of the nucleus.