Search

Fanpy is a free and open-source Python library for developing and testing multideterminant wavefunctions and related ab initio methods in electronic structure theory. The main use of Fanpy is to quickly prototype new methods by making it easier to convert the mathematical formulation of a new wavefunction ansatze to a working implementation. Fanpy is designed based on our recently introduced Flexible Ansatz for N-electron Configuration Interaction (FANCI) framework, where multideterminant wavefunctions are represented by their overlaps with Slater determinants of orthonormal spin-orbitals. In the simplest case, a new wavefunction ansatz can be implemented by simply writing a function for evaluating its overlap with an arbitrary Slater determinant. Fanpy is modular in both implementation and theory: the wavefunction model, the system's Hamiltonian, and the choice of objective function are all independent modules. This modular structure makes it easy for users to mix and match different methods and for developers to quickly explore new ideas. Fanpy is written purely in Python with standard dependencies, making it accessible for various operating systems. In addition, it adheres to principles of modern software development, including comprehensive documentation, extensive testing, quality assurance, and continuous integration and delivery protocols. This article is considered to be the official release notes for the Fanpy library.

Melamine is a precursor and building block for graphitic carbon nitride (g-CN) materials, a group of layered materials showing great promise for catalytic applications. The synthetic pathway to g-CN includes a polycondensation reaction of melamine by evaporation of ammonia. Melamine molecules in the crystal organize into wave-like planes with an interlayer distance of 3.3 Å similar to that of g-CN. Here we present an extensive investigation of the experimental electron density of melamine obtained from modelling of synchrotron radiation X-ray single-crystal diffraction data measured at 25 K with special focus on the molecular geometry and intermolecular interactions. Both intra- and interlayer structures are dominated by hydrogen bonding and π-interactions. Theoretical gas-phase optimizations of the experimental molecular geometry show that bond lengths and angles for atoms in the same chemical environment (C−N bonds in the ring, amine groups) differ significantly more for the experimental geometry than for the gas-phase-optimized geometries, indicating that intermolecular interactions in the crystal affects the molecular geometry. In the experimental crystal geometry, one amine group has significantly more sp 3-like character than the others, hinting at a possible formation mechanism of g-CN. Topological analysis and energy frameworks show that the nitrogen atom in this amine group participates in weak intralayer hydrogen bonding. We hypothesize that melamine condenses to g-CN within the layers and that the unique amine group plays a key role in the condensation process.

Eigenvectors of the reduced Bardeen-Cooper-Schrieffer (BCS) Hamiltonian, Richardson-Gaudin (RG) states, are used as a variational wavefunction ansatz for strongly correlated electronic systems. These states are geminal products whose coefficients are solutions of non-linear equations. Previous results showed an un-physical apparent avoided crossing in ground state dissociation curves for hydrogen chains. In this paper, it is shown that each seniority-zero state of the molecular Coulomb Hamiltonian corresponds directly to an RG state. However, the seniority-zero ground state does not correspond to the ground state of a reduced BCS Hamiltonian. The difficulty is in choosing the correct RG state. The systems studied showed a clear choice, and we expect that it should always be possible to reason physically which state to choose. Published under an exclusive license by AIP Publishing.

Trends in bond angle are identified in a systematic study of more than a thousand symmetric A(2)B triatomic molecules. We show that, in series where atoms A and B are each varied within a group, the following trends hold: (1) the A-B-A bond angle decreases for more polarizable central atoms B, and (2) the A-B-A angle increases for more polarizable outer atoms A. The physical underpinning is provided by the extended Debye polarizability model for the chemical bond angle, hence our present findings also serve as validation of this simple classical model. We use experimental bond angles from the literature and, where not available, we optimize molecular geometries with quantum chemical methods, with an open mind with regards to the stability of these molecules. We consider main group elements up to and including the sixth period of the periodic table.

We are grateful for constructive comments from K.-L. Cai, O. Denisov, F. de Soto, T. Frederico, J. Friedrich, C. Mezrag, V. Mokeev, W.-D. Nowak, C. Quintans, G. Salme and J. Segovia. Work supported by: National Natural Science Foundation of China (Grant nos. 12135007, 11805097); Helmholtz-Zentrum Dresden-Rossendorf High Potential Programme; Spanish Ministry of Science and Innovation (MICINN) (Grant nos. PID2020-113334GB, PID2019-107844GBC22); Generalitat Valenciana (Grant no. Prometeo/2019/087); Junta de Andalucia (Grant nos. P18-FR-5057, UHU-1264517, UHU EPIT2021); and STRONG-2020 "The strong interaction at the frontier of knowledge: fundamental research and applications" which received funding from the European Union's Horizon 2020 research and innovation programme (grant no. 824093)., Analyses of the pion valence-quark distribution function (DF), u(pi) (x; sigma), which explicitly incorporate the behaviour of the pion wave function prescribed by quantum chromodynamics (QCD), predict u(pi) (x similar or equal to 1; sigma) similar to (1 - x)(beta(sigma)), beta(sigma greater than or similar to m(p)) > 2, where mp is the proton mass. Nevertheless, more than forty years after the first experiment to collect data suitable for extracting the x similar or equal to 1 behaviour of up, the empirical status remains uncertain because some methods used to fit existing data return a result for up that violates this constraint. Such disagreement entails one of the following conclusions: the analysis concerned is incomplete; not all data being considered are a true expression of qualities intrinsic to the pion; or QCD, as it is currently understood, is not the theory of strong interactions. New, precise data are necessary before a final conclusion is possible. In developing these positions, we exploit a single proposition, viz. there is an effective charge which defines an evolution scheme for parton DFs that is all-orders exact. This proposition has numerous corollaries, which can be used to test the character of any DF, whether fitted or calculated., National Natural Science Foundation of China (NSFC) 12135007 11805097, Spanish Government European Commission PID2020-113334GB PID2019-107844GBC22, Generalitat Valenciana European Commission, General Electric Prometeo/2019/087, Junta de Andalucia P18-FR-5057 UHU-1264517 UHU EPIT2021, STRONG-2020 "The strong interaction at the frontier of knowledge: fundamental research and applications" - European Union's Horizon 2020 research and innovation programme 824093, Helmholtz-Zentrum Dresden-Rossendorf High Potential Programme

We extend exactly solvable model of a one-dimensional non-relativistic canonical semiconfined quantum harmonic oscillator with a mass that varies with position to the case where an external homogeneous field is applied. The problem is still exactly solvable and the analytic expression of the wave functions of the stationary states is expressed by means of generalised Laguerre polynomials, too. Unlike the case without any external field, when the energy spectrum completely overlaps with the energy spectrum of the standard non-relativistic canonical quantum harmonic oscillator, the energy spectrum is now still equidistant but depends on the semiconfinement parameter a. We also compute probabilities of the transitions for the model under the external field and discuss limit cases for the energy spectrum, wave functions and probabilities of transitions, when the semiconfinement parameter a goes to infinity.

We compute the differential yield for quark anti-quark dijet production in high-energy electron-proton and electron-nucleus collisions at small $x$ as a function of the relative momentum $\boldsymbol{P}_\perp$ and momentum imbalance $\boldsymbol{k}_\perp$ of the dijet system for different photon virtualities $Q^2$, and study the elliptic and quadrangular anisotropies in the relative angle between $\boldsymbol{P}_\perp$ and $\boldsymbol{k}_\perp$. We review and extend the analysis in [1], which compared the results of the Color Glass Condensate (CGC) with those obtained using the transverse momentum dependent (TMD) framework. In particular, we include in our comparison the improved TMD (ITMD) framework, which resums kinematic power corrections of the ratio $k_\perp$ over the hard scale $Q_\perp$. By comparing ITMD and CGC results we are able to isolate genuine higher saturation contributions in the ratio $Q_s/Q_\perp$ which are resummed only in the CGC. These saturation contributions are in addition to those in the Weizs\"ackerWilliams gluon TMD that appear in powers of $Q_s/k_\perp$. We provide numerical estimates of these contributions for inclusive dijet production at the future Electron-Ion Collider, and identify kinematic windows where they can become relevant in the measurement of dijet and dihadron azimuthal correlations. We argue that such measurements will allow the detailed experimental study of both kinematic power corrections and genuine gluon saturation effects., Comment: v3: typos corrected post-publication. v2: 52 pages, 13 figures, matches published version

We study coherent diffractive photon and vector meson production in electron-proton and electron-nucleus collisions within the Color Glass Condensate effective field theory. We show that electron-photon and electron-vector meson azimuthal angle correlations are sensitive to non-trivial spatial correlations in the gluon distribution of the target, and perform explicit calculations using spatially dependent McLerran-Venugopalan initial color charge configurations coupled to the numerical solution of small $x$ JIMWLK evolution equations. We compute the cross-section differentially in $Q^2$ and $|t|$ and find sizeable anisotropies in the electron-photon and electron-$\mathrm{J}/��$ azimuthal correlations ($v_{1,2} \approx 2 - 10 \%$) in electron-proton collisions for the kinematics of the future Electron-Ion Collider. In electron-gold collisions these modulations are found to be significantly smaller ($v_{1,2}, 32 figures, 18 figures

The covalent nature of the low-barrier N−H−N hydrogen bonds in the negative thermal expansion material H 3[Co(CN) 6] has been established by using a combination of X-ray and neutron diffraction electron density analysis and theoretical calculations. This finding explains why negative thermal expansion can occur in a material not commonly considered to be built from rigid linkers. The pertinent hydrogen atom is located symmetrically between two nitrogen atoms in a double-well potential with hydrogen above the barrier for proton transfer, thus forming a low-barrier hydrogen bond. Hydrogen is covalently bonded to the two nitrogen atoms, which is the first experimentally confirmed covalent hydrogen bond in a network structure. Source function calculations established that the present N−H−N hydrogen bond follows the trends observed for negatively charge-assisted hydrogen bonds and low-barrier hydrogen bonds previously established for O−H−O hydrogen bonds. The bonding between the cobalt and cyanide ligands was found to be a typical donor–acceptor bond involving a high-field ligand and a transition metal in a low-spin configuration.

We develop an extension of the Gutzwiller Approximation (GA) formalism that includes the effects of Coulomb interactions of arbitrary range (including density density, exchange, pair hopping and Coulomb assisted hopping terms). This formalism reduces to the ordinary GA formalism for the multi-band Hubbard models in the presence of only local interactions. This is accomplished by combining the 1/z expansion -- where z is the coordination number and only the leading order terms contribute in the limit of infinite dimensions -- with a P^2-I expansion, where P is the Gutzwiller projector. The method is conveniently formulated in terms of a Gutzwiller Lagrange function. We apply our theory to the extended single band Hubbard model. Similarly to the usual Brinkman-Rice mechanism we find a Mott transition. A valence skipping transition is observed, where the occupation of the empty and doubly occupied states for the Gutzwiller wavefunction is enhanced with respect to the uncorrelated Slater determinant wavefunction., 13+23 pages (1+4 figures) comments welcome; v2: typos corrected, submitted to PRB; v3: Added Section II which presents our results in a simple intuitive manner

In previous work on the performance of multiconfigurational second-order perturbation theory (CASPT2) in describing spin state energetics in first-row transition metal systems [ Pierloot et al. J. Chem. Theory Comput. 2017 , 13 , 537 - 553 ], we showed that standard CASPT2 works well for valence correlation but does not describe the metal semicore (3s3p) correlation effects accurately. This failure is partially responsible for the well-known bias toward high-spin states of CASPT2. In this paper, we expand our previous work and show that this bias could be partly removed with a combined CASPT2/CC approach: using high-quality CASPT2 with extensive correlation-consistent basis sets for valence correlation and low-cost CCSD(T) calculations with minimal basis sets for the metal semicore (3s3p) correlation effects. We demonstrate that this approach is efficient by studying the spin state energetics of a series of iron complexes modeling important intermediates in oxidative catalytic processes in chemistry and biochemistry. On the basis of a comparison with bare CCSD(T) results from this and previous work, the average error of the CASPT2/CC approach is estimated at around 2 kcal mol-1 in favor of high spin states. ispartof: JOURNAL OF CHEMICAL THEORY AND COMPUTATION vol:14 issue:5 pages:2446-2455 ispartof: location:United States status: published

For the design of efficient singlet fission chromophores, knowledge of the factors that govern the singlet fission rate is important. This rate is approximately proportional to the electronic coupling between the lowest (diabatic) spin singlet state that is populated following photoexcitation state and a so-called (TT)-T-1 state. The latter state is characterised by two triplets, each localised on one of two neighbouring molecules, which are coupled into a singlet. Here, we show the applicability of a nonorthogonal configuration interaction approach for the calculation of this electronic coupling. The advantages of this rigorous approach are that (1) the coupling can be calculated directly, (2) it includes important correlation and orbital relaxation effects, and (3) it has a clear chemical interpretation in terms of molecular states. This approach is applied to calculate the electronic coupling for a biradicaloid molecule, viz, the bis(inner salt) of 2,5-dihydroxy-1,4-dimethyl-pyrazinium. The biradicaloid molecule is, based on the energetic criteria, a promising candidate for singlet fission. We show that the electronic coupling between the molecules is also sufficiently large for singlet fission, rendering molecules based on this chemical moiety interesting singlet fission chromophores. (C) 2017 The Author(s). Published by Elsevier B.V.

In this study, single crystals of urea ninhydrin monohydrate (UNMH) have been grown by slow evaporation method. The grown crystals were characterized by FT-IR, FT-Raman and UV-Vis-NIR spectroscopies. The Kurtz and Perry powder method was employed to confirm the near-zero SHG efficiency of the as-grown centrosymmetric UNMH crystal. The third order nonlinearity of the crystal has been studied by the open aperture Z-scan method. The nonlinear absorption coefficient is calculated and the potentiality of UNMH in optical limiting applications is identified. The molecular geometry and the origin of optical non-linearity at the molecular level have been investigated by the density functional theory. The normal coordinate analysis was carried out to assign the molecular vibrational modes. Vibrational spectral studies confirms the presence of weak O-H ... O and moderate O-H ... O type hydrogen bonds in the molecule as well as O-H ... O, N-H ... O and blue-shifted C-H ... O type H-bonds in the crystal. The intramolecular charge transfer interactions and the electronic absorption mechanisms have been discussed. The static and the dynamic values of hyperpolarizabilities for UNMH were estimated theoretically by DFT methods. (C) 2017 Elsevier B.V. All rights reserved.

We report on accurate and efficient calculations of vibrationally resolved emission spectra for oligothiophenes from anharmonic vibrational configuration interaction wave-function calculations in reduced vibrational spaces. These reduced spaces are chosen based on the independent mode displaced harmonic oscillator model. Good agreement with experiment is obtained for all-trans oligothiophenes with two to five rings also when employing only a few active modes. Vibrational modes incorporating inter-ring carbon-carbon stretches and a ring breathing mode are found to be the main players in the vibrational progression for the emission from the first excited electronic state for all investigated oligothiophene derivatives. The presented framework is here illustrated for oligothiophenes, but we have made no underlying system-dependent assumptions and believe it to become a valuable tool for the rational design of fluorescence biomarkers.

Building on the success of Quantum Monte Carlo techniques such as diffusion Monte Carlo, alternative stochastic approaches to solve electronic structure problems have emerged over the past decade. The full configuration interaction quantum Monte Carlo (FCIQMC) method allows one to systematically approach the exact solution of such problems, for cases where very high accuracy is desired. The introduction of FCIQMC has subsequently led to the development of coupled cluster Monte Carlo (CCMC) and density matrix quantum Monte Carlo (DMQMC), allowing stochastic sampling of the coupled cluster wave function and the exact thermal density matrix, respectively. In this Article, we describe the HANDE-QMC code, an open source implementation of FCIQMC, CCMC and DMQMC, including initiator and semistochastic adaptations. We describe our code and demonstrate its use on three example systems; a molecule (nitric oxide), a model solid (the uniform electron gas), and a real solid (diamond). An illustrative tutorial is also included. Public domain – authored by a U.S. government employee

A methodology is proposed to investigate electron transfer reactions between redox-active biomolecular systems (e.g. a protein) and inorganic surfaces. The whole system is modelled at the atomistic level using classical molecular dynamics - making an extensive sampling of the system's configurations possible - and the energies associated with the redox-active complex reduction are calculated using a hybrid quantum/classical approach along the molecular dynamics trajectory. The non-adiabaticity is introduced a posteriori using a Monte Carlo approach based on the Landau-Zener theory extended to treat a metal surface. This approach thus allows us to investigate the role of the energy fluctuations, determined by the dynamical evolution of the system, as well as the role of non-adiabaticity in affecting the kinetic rate of the electron transfer reaction. Most notably, it allows us to investigate the two contributions separately, hence achieving a detailed picture of the mechanisms that determine the rate. The analysis of the system configurations also allows us to relate the estimated electronic coupling to the structural and dynamic properties of the system. As a test case, the methodology is here applied to study the electron transfer reaction between cytochrome c and a gold surface. The results obtained explain the different electron transfer rates experimentally measured for two different concentrations of proteins on the electrode surface.

International audience; The fast and reliable determination of wave functions and electron densities of macromolecules has been one of the goals of theoretical chemistry for a long time, and in this context, several linear scaling techniques have been successfully devised over the years. Different approaches have been adopted to tackle this problem, and one of them exploits the fact that, according to the traditional chemical perception, molecules can be seen as constituted of recurring units (e.g., functional groups) with well-defined chemical features. This has led to the development of methods in which the global wave functions or electron densities of macromolecules are obtained by simply transferring density matrices or fuzzy electron densities associated with molecular fragments. In this context, we propose an alternative strategy that aims at quickly reconstructing wave functions and electron densities of proteins through the transfer of extremely localized molecular orbitals (ELMOs), which are orbitals strictly localized on small molecular units and, for this reason, easily transferable from molecule to molecule. To accomplish this task we have constructed original libraries of ELMOs that cover all the possible elementary fragments of the 20 natural amino acids in all their possible protonation states and forms. Our preliminary test calculations have shown that, compared to more traditional methods of quantum chemistry, the transfers from the novel ELMO databanks allow to obtain wave function and electron densities of large polypeptides and proteins at a significantly reduced computational cost. Furthermore, notwithstanding expected discrepancies, the obtained electron distributions and electrostatic potentials are in very good agreement with those obtained at Hartree–Fock and density functional theory (DFT) levels. Therefore, the results encourage to use the new libraries as alternatives to the popular pseudoatom-databases of crystallography in the refinement of crystallographic structures of macromolecules. In particular, in this context, we have already envisaged the coupling of the ELMO databanks with the promising Hirshfeld atom refinement technique to extend the applicability of the latter to very large systems.

We demonstrate that charge migration can be ‘engineered’ in arbitrary molecular systems if a single localised orbital – that diabatically follows nuclear displacements – is ionised. Specifically, we describe the use of natural bonding orbitals in Complete Active Space Configuration Interaction (CASCI) calculations to form cationic states with localised charge, providing consistently well-defined initial conditions across a zero point energy vibrational ensemble of molecular geometries. In Ehrenfest dynamics simulations following localised ionisation of -electrons in model polyenes (hexatriene and decapentaene) and -electrons in glycine, oscillatory charge migration can be observed for several femtoseconds before dephasing. Including nuclear motion leads to slower dephasing compared to fixed-geometry electron-only dynamics results. For future work, we discuss the possibility of designing laser pulses that would lead to charge migration that is experimentally observable, based on the proposed diabatic orbital approach. KEYWORDS: Ehrenfest method, coupled electron-nuclear dynamics, charge migration, localised orbital

Accurate characterization of the high-lying excited state manifolds of organic molecules is of fundamental importance for the interpretation of the rich response detected in time-resolved nonlinear electronic spectroscopies. Here, we have characterized the singlet excited state manifold of benzophenone (BP), a versatile organic photoinitiator and a well-known DNA photosensitizer. Benchmarks of various multiconfigurational/multireference (RASSCF/PT2) and time-dependent density functional theory (TD-DFT) approaches allowed assignments of experimental linear absorption signals of BP in the ultraviolet (UV) region, with unprecedented characterization of ground state absorptions in the far UV. Experimental transient absorption spectra obtained by UV-vis pump-probe spectroscopy at very short time delays are shown to be directly comparable to theoretical estimates of excited state absorptions (from the low-lying nOπ∗ and π π∗ singlet states) in the Franck-Condon region. Multireference computations provided reliable interpretation of the PP spectra, with TD-DFT results yielding a fair agreement as long as electronic transitions featuring double excitations contributions are not involved. These results lay the groundwork for further computational studies and interpretation of experimental nonlinear electronic spectra of benzophenone in more complex systems, such as BP/DNA adducts.

In methods like geminal-based approaches or coupled cluster that are solved using the projected Schr\"odinger equation, direct computation of the 2-electron reduced density matrix (2-RDM) is impractical and one falls back to a 2-RDM based on response theory. However, the 2-RDMs from response theory are not $N$-representable. That is, the response 2-RDM does not correspond to an actual physical $N$-electron wave function. We present a new algorithm for making these non-$N$-representable 2-RDMs approximately $N$-representable, i.e. it has the right symmetry and normalization and it fulfills the $P$-, $Q$- and $G$-conditions. Next to an algorithm which can be applied to any 2-RDM, we have also developed a 2-RDM optimization procedure specifically for seniority-zero 2-RDMs. We aim to find the 2-RDM with the right properties that is the closest (in the sense of the Frobenius norm) to the non-N-representable 2-RDM by minimizing the square norm of the difference between the initial 2-RDM and the targeted 2-RDM under the constraint that the trace is normalized and the 2-RDM, $Q$- and $G$-matrices are positive semidefinite, i.e. their eigenvalues are non-negative. Our method is suitable for fixing non-N-respresentable 2-RDMs which are close to being N-representable. Through the N-representability optimization algorithm we add a small correction to the initial 2-RDM such that it fulfills the most important N-representability conditions., Comment: 13 pages, 8 figures

A simple algorithm for the Kohn–Sham inversion problem is presented. The method is found to converge toward a nearby $$\nu$$ -representable Kohn–Sham density irrespective of the fact whether the initial target density has been $$\nu$$ -representable or not. For the proposed procedure, the target density can be of general nature. The algorithm can handle Hartree–Fock and post-Hartree–Fock, spin-unpolarized and polarized states equally well. Additionally, experimental densities and even general gedanken densities can be treated. The algorithm is easy to implement and does not require an additional procedure to adjust eigenvalues.

The high resolution near edge X-ray absorption fine structure spectrum of nitrogen displays the vibrational structure of the core-excited states. This makes nitrogen well suited for assessing the accuracy of different electronic structure methods for core excitations. We report high resolution experimental measurements performed at the SOLEIL synchrotron facility. These are compared with theoretical spectra calculated using coupled cluster theory and algebraic diagrammatic construction theory. The coupled cluster singles and doubles with perturbative triples model known as CC3 is shown to accurately reproduce the experimental excitation energies as well as the spacing of the vibrational transitions. The computational results are also shown to be systematically improved within the coupled cluster hierarchy, with the coupled cluster singles, doubles, triples, and quadruples method faithfully reproducing the experimental vibrational structure. This is the authors’ accepted and refereed manuscript to the article. This article may be downloaded for personal use only. Any other use requires prior permission of the author and AIP Publishing. The following article appeared in Journal of Chemical Physics and may be found at https://doi.org/10.1063/1.5011148

We report an implementation of non-adiabatic coupling (NAC) forces within the equation-of-motion coupled-cluster with single and double excitations (EOM-CCSD) framework via the summed-state approach. Using illustrative examples, we compare NAC forces computed with EOM-CCSD and multi-reference (MR) wave functions (for selected cases, we also consider configuration interaction singles). In addition to the magnitude of the NAC vectors, we analyze their direction, which is important for the calculations of the rate of non-adiabatic transitions. Our benchmark set comprises three doublet radical-cations (hexatriene, cyclohexadiene, and uracil), neutral uracil, and sodium-doped ammonia clusters. When the characters of the states agree among different methods, we observe good agreement between the respective NAC vectors, both in the Franck-Condon region and away. In the cases of large discrepancies between the methods, the disagreement can be attributed to the difference in the states' character, which, in some cases, is very sensitive to electron correlation, both within single-reference and multi-reference frameworks. The numeric results confirm that the accuracy of NAC vectors depends critically on the quality of the underlying wave functions. Within their domain of applicability, EOM-CC methods provide a viable alternative to MR approaches.

We present a new variational tree tensor network state (TTNS) ansatz, the three-legged tree tensor network state (T3NS). Physical tensors are interspersed with branching tensors. Physical tensors have one physical index and at most two virtual indices, as in the matrix product state (MPS) ansatz of the density matrix renormalization group (DMRG). Branching tensors have no physical index, but up to three virtual indices. In this way, advantages of DMRG, in particular a low computational cost and a simple implementation of symmetries, are combined with advantages of TTNS, namely incorporating more entanglement. Our code is capable of simulating quantum chemical Hamiltonians, and we present several proof-of-principle calculations on LiF, N$_2$ and the bis($\mu$-oxo) and $\mu - \eta^2 : \eta^2$ peroxo isomers of $[\mathrm{Cu}_2\mathrm{O}_2]^{2+}$., Comment: 14 pages, 8 figures

With an increase of computational capabilities, ab initio molecular dynamics becomes the natural choice for exploring the nuclear dynamics of solids. As based on classical mechanics, the validity of this approach is, in-principle, limited to the high-T regime, whilst low-temperature simulations require inclusion of quantum effects. The methods commonly used to account for nuclear quantum effects are based on the path-integral formalism, which become, however, particularly time consuming when high accuracy methods are used for calculating forces. Recently, new efficient alternative approaches to account for quantum nature of nuclei have been proposed, using so-called quantum thermostats. In this work, we examine the simulations performed with the quantum colored-noise thermostat introduced by Ceriotti [Phys. Rev. Lett., 103: 030603, 2009]. We present the tests of portable implementation of the quantum thermostat in the ABIN program, which has been extended to periodic systems through the interface to CASTEP, a leading spectroscopy-oriented plane-wave density functional theory code. The range of applicability of quantum-thermostatted molecular dynamics simulations for the interpretation of neutron scattering data was examined and compared to classical molecular dynamics and lattice-dynamics simulations, using solid formic acid case as a test bed. We find that the approach is particularly useful for the modeling of low-temperature inelastic neutron scattering spectra as well as provides some theoretical estimate for the low-limit of the mean kinetic energy. While finding the quantum-thermostat to seriously affect the dynamic properties of the title system, we illustrate to which extent the unperturbed response can be successfully recovered.

The role of a strong organobase, guanidine, in sulfuric acid-driven new-particle formation is studied using state-of-the-art quantum chemical methods and molecular cluster formation simulations. Cluster formation mechanisms at the molecular level are resolved, and theoretical results on cluster stability are confirmed with mass spectrometer measurements. New-particle formation from guanidine and sulfuric acid molecules occurs without thermodynamic barriers under studied conditions, and clusters are growing close to a 1:1 composition of acid and base. Evaporation rates of the most stable clusters are extremely low, which can be explained by the proton transfers and symmetrical cluster structures. We compare the ability of guanidine and dimethylamine to enhance sulfuric acid-driven particle formation and show that more than 2000-fold concentration of dimethylamine is needed to yield as efficient particle formation as in the case of guanidine. At similar conditions, guanidine yields 8 orders of magnitude higher particle formation rates compared to dimethylamine. Highly basic compounds such as guanidine may explain experimentally observed particle formation events at low precursor vapor concentrations, whereas less basic and more abundant bases such as ammonia and amines are likely to explain measurements at high concentrations.

Many population analysis methods are based on the precept that molecules should be built from fragments (typically atoms) that maximally resemble the isolated fragment. The resulting molecular building blocks are intuitive (because they maximally resemble well-understood systems) and transferable (because if two molecular fragments both resemble an isolated fragment, they necessarily resemble each other). Information theory is one way to measure the deviation between molecular fragments and their isolated counterparts, and it is a way that lends itself to interpretation. For example, one can analyze the relative importance of electron transfer and polarization of the fragments. We present key features, advantages, and disadvantages of the information-theoretic approach. We also codify existing information-theoretic partitioning methods in a way, that clarifies the enormous freedom one has within the information-theoretic ansatz.

We give a brief description of the variational and diffusion quantum Monte Carlo methods and their application to the study of self-interstitial defects in silicon. The diffusion quantum Monte Carlo calculations give formation energies for the most stable defects of about 4.9 eV, which is considerably larger than the values obtained in density functional theory methods. The quantum Monte Carlo results indicate a value for the formation+migration energy of the self-interstitial contribution to self-diffusion of about 5 eV, which is consistent with the experimental data.

We present the inner products of eigenstates in integrable Richardson-Gaudin models from two different perspectives and derive two classes of Gaudin-like determinant expressions for such inner products. The requirement that one of the states is on-shell arises naturally by demanding that a state has a dual representation. By implicitly combining these different representations, inner products can be recast as domain wall boundary partition functions. The structure of all involved matrices in terms of Cauchy matrices is made explicit and used to show how one of the classes returns the Slavnov determinant formula. This framework provides a further connection between two different approaches for integrable models, one in which everything is expressed in terms of rapidities satisfying Bethe equations, and one in which everything is expressed in terms of the eigenvalues of conserved charges, satisfying quadratic equations., 21+16 pages, minor revisions compared to the previous version

We present the extension of variational Monte Carlo (VMC) to the calculation of electronic excitation energies and oscillator strengths using time-dependent linear-response theory. By exploiting the analogy existing between the linear method for wave-function optimisation and the generalised eigenvalue equation of linear-response theory, we formulate the equations of linear-response VMC (LR-VMC). This LR-VMC approach involves the first-and second-order derivatives of the wave function with respect to the parameters. We perform first tests of the LR-VMC method within the Tamm-Dancoff approximation using single-determinant Jastrow-Slater wave functions with different Slater basis sets on some singlet and triplet excitations of the beryllium atom. Comparison with reference experimental data and with configuration-interaction-singles (CIS) results shows that LR-VMC generally outperforms CIS for excitation energies and is thus a promising approach for calculating electronic excited-state properties of atoms and molecules., Advances in Quantum Chemistry, 2017, Novel Electronic Structure Theory: General Innovations and Strongly Correlated Systems

In analogy to the previously reported tetranuclear magnesium hydride cluster with a bridged dianionic bis-beta-diketiminate ligand, a related zinc hydride cluster has been prepared. The crystal structures of these magnesium and zinc hydride complexes are similar: the metal atoms are situated at the corners of a tetrahedron in which the vertices are bridged either by dianionic bis-b-diketiminate ligands or hydride ions. Both structures are retained in solution and show examples of H-center dot center dot center dot H- NMR coupling (Mg: 8.5 Hz; Zn: 16.0 Hz). The zinc hydride cluster [NN-(ZnH)(2)](2) thermally decomposes at 90 degrees C and releases 1.8 equivalents of H-2. In contrast to magnesium hydride clusters, there is no apparent relationship between cluster size and thermal decomposition temperature for the zinc hydrides. DFT calculations reproduced the structure of the zinc hydride cluster reasonably well and charge density analysis showed no bond paths between the hydride ions. This contrasts with calculations on the analogous magnesium hydride cluster in which a counter-intuitive H-center dot center dot center dot H- bond path was observed. Forcing a reduced H-center dot center dot center dot H- distance in the zinc hydride cluster, however, gave rise to a H-center dot center dot center dot H- bond path. Such weak interactions could play a role in H-2 desorption. The presumed molecular product after H-2 release, a Zn(I) cluster, could not be characterized experimentally but DFT calculations predicted a cluster with two localized Zn-Zn bonds.

Unimolecular decompositions of neutral (NH2CHO) and protonated (NH3CHO+) formamide, an active precursor of biomolecules in prebiotic chemistry, are investigated in the ground (S-0) and first triplet (T-1) and singlet (S-1) excited states. Different decomposition channels including the homolytic bond dissociations, dehydration, decarbonylation, dehydrogenation, etc., are explored using coupled-cluster theory (CCSD(T)/CBS method) for both S-0 and T-1 states and RASPT2(18,15)/6-31 G(d,p) computations for the S-1 state. On S-1 and T-1 energy surfaces, formamide preferentially follows C-N homolytic bond cleavages forming NH2 + HCO radical pairs. Formation of HCN and HNC from dehydration of neutral and protonated formamide via formimic acid and aminohydroxymethylene isomers has higher energy barriers. A strong stabilization upon triplet excitation of the two latter isomers significantly facilitates the interconversions between isomers, and thus considerably reduces the energy barriers for dehydration pathways. The most probable pathways for HCN and HNC generation are found to be dehydration of formamide in the T-1 state. Dehydration pathways from the neutral S-1 and protonated T-1 forms lead to stable complexes of HCN and HNC with water but are associated with large energy barriers. Overall, in the lower-lying excited states of either neutral or protonated formamide, dehydration is not competitive with homolytic C-N bond cleavages, which finally lead to formation of CO.

The electric dipole polarizability quantifies the low-energy behaviour of the dipole strength and is related to critical observables such as the radii of the proton and neutron distributions. Its computation is challenging because most of the dipole strength lies in the scattering continuum. In this paper we combine integral transforms with the coupled-cluster method and compute the dipole polarizability using bound-state techniques. Employing different interactions from chiral effective field theory, we confirm the strong correlation between the dipole polarizability and the charge radius, and study its dependence on three-nucleon forces. We find good agreement with data for the 4He, 40Ca, and 16O nuclei, and predict the dipole polarizability for the rare nucleus 22O., 12 pages, 11 figures

A first quantized free photon is a complex massless vector field $A=(A^\mu)$ whose field strength satisfies Maxwell's equations in vacuum. We construct the Hilbert space $\mathscr{H}$ of the photon by endowing the vector space of the fields $A$ in the temporal-Coulomb gauge with a positive-definite and relativistically invariant inner product. We give an explicit expression for this inner product, identify the Hamiltonian for the photon with the generator of time translations in $\mathscr{H}$, determine the operators representing the momentum and the helicity of the photon, and introduce a chirality operator whose eigenfunctions correspond to fields having a definite sign of energy. We also construct a position operator for the photon whose components commute with each other and with the chirality and helicity operators. This allows for the construction of the localized states of the photon with a definite sign of energy and helicity. We derive an explicit formula for the latter and compute the corresponding electric and magnetic fields. These turn out to diverge not just at the point where the photon is localized but on a plane containing this point. We identify the axis normal to this plane with an associated symmetry axis, and show that each choice of this axis specifies a particular position operator, a corresponding position basis, and a position representation of the quantum mechanics of photon. In particular, we examine the position wave functions determined by such a position basis, elucidate their relationship with the Riemann-Silberstein and Landau-Peierls wave functions, and give an explicit formula for the probability density of the spatial localization of the photon., Comment: 35 pages, 2 figures, slightly expanded published version

We present results of the first ab initio lattice QCD calculation of the normalization constants and first moments of the leading twist distribution amplitudes of the full baryon octet, corresponding to the small transverse distance limit of the associated S-wave light-cone wave functions. The P-wave (higher twist) normalization constants are evaluated as well. The calculation is done using N-f = 2 + 1 flavors of dynamical (clover) fermions on lattices of different volumes and pion masses down to 222 MeV. Significant SU(3) flavor symmetry violation effects in the shape of the distribution amplitudes are observed., SCOAP3

Semilocal density functional theory is the most used computational method for electronic structure calculations in theoretical solid-state physics and quantum chemistry of large systems, providing good accuracy with a very attractive computational cost. Nevertheless, because of the non-locality of the exchange-correlation hole outside a metal surface, it was always considered inappropriate to describe the correct surface asymptotics. Here, we derive, within the semilocal density functional theory formalism, an exact condition for the image-like surface asymptotics of both the exchange-correlation energy per particle and potential. We show that this condition can be easily incorporated into a practical computational tool, at the simple meta-generalized-gradient approximation level of theory. Using this tool, we also show that the Airy-gas model exhibits asymptotic properties that are closely related to the ones at metal surfaces. This result highlights the relevance of the linear effective potential model to the metal surface asymptotics., Comment: 6 pages, 4 figures

[Image: see text] The flexibility of three MIL-47(V)-type materials (MIL-47, COMOC-2, and COMOC-3) has been explored by constructing the pressure versus volume and free energy versus volume profiles at various temperatures ranging from 100 to 400 K. This is done with first-principles-based force fields using the recently proposed QuickFF parametrization protocol. Specific terms were added for the materials at hand to describe the asymmetry of the one-dimensional vanadium-oxide chain and to account for the flexibility of the organic linkers. The force fields are used in a series of molecular dynamics simulations at fixed volumes but varying unit cell shapes. The three materials show a distinct pressure–volume behavior, which underlines the ability to tune the mechanical properties by varying the linkers toward different applications such as nanosprings, dampers, and shock absorbers.

Density matrix embedding theory (DMET) provides a theoretical framework to treat finite fragments in the presence of a surrounding molecular or bulk environment, even when there is significant correlation or entanglement between the two. In this work, we give a practically oriented and explicit description of the numerical and theoretical formulation of DMET. We also describe in detail how to perform self-consistent DMET optimizations. We explore different embedding strategies with and without a self-consistency condition in hydrogen rings, beryllium rings, and a sample S$_{\text{N}}$2 reaction. The source code for the calculations in this work can be obtained from \url{https://github.com/sebwouters/qc-dmet}., Comment: 41 pages, 10 figures