1,721,037 research outputs found
A black-box, general purpose quadratic self-consistent field code with and without Cholesky decomposition of the two-electron integrals
Full text linkWe present the implementation of a quadratically convergent self-consistent field (QCSCF) algorithm based on an adaptive trust-radius optimisation scheme for restricted open-shell Hartree–Fock (ROHF), restricted Hartree–Fock (RHF), and unrestricted Hartree–Fock (UHF) references. The algorithm can exploit Cholesky decomposition (CD) of the two-electron integrals to allow calculations on larger systems. The most important feature of the QCSCF code lies in its black-box nature–probably the most important quality desired by a generic user. As shown for pilot applications, it does not require one to tune the self-consistent field (SCF) parameters (damping, Pulay's DIIS, and other similar techniques) in difficult-to-converge molecules. Also, it can be used to obtain a very tight convergence with extended basis sets–a situation often needed when computing high-order molecular properties–where the standard SCF algorithm starts to oscillate. Nevertheless, trouble may appear even with a QCSCF solver. In this respect, we discuss what can go wrong, focusing on the multiple UHF solutions of ortho-benzyne
A novel coupled-cluster singles and doubles implementation that combines the exploitation of point-group symmetry and Cholesky decomposition of the two-electron integrals
No full textA novel implementation of the coupled-cluster singles and doubles (CCSD) approach is presented that is specifically tailored for the treatment of large symmetric systems. It fully exploits Abelian point-group symmetry and the use of the Cholesky decomposition of the two-electron repulsion integrals. In accordance with modern CCSD algorithms, we propose two alternative strategies for the computation of the so-called particle–particle ladder term. The code is driven toward the optimal choice depending on the available hardware resources. As a large-scale application, we computed the frozen-core correlation energy of buckminsterfullerene (C60) with a polarized valence triple-zeta basis set (240 correlated electrons in 1740 orbitals)
Deuterium hyperfine splittings in the rotational spectrum of NH2D as revealed by Lamb-dip spectroscopy
Full text linkIn the context of radio-astronomical observations, laboratory experiments constitute a cornerstone in the interpretation of rich line surveys due to the concomitant presence of numerous emitting molecules. Here, we report the investigation of three different rotational transitions of mono-deuterated ammonia (NH2D), a species of astrophysical interest, for which the contribution of the deuterium nuclear spin to the rotational spectrum has been resolved for the first time in the millimeter- and submillimeter-wave domain. The effect of hyperfine interactions on the rotational spectrum has been unveiled by a combined theoretical and experimental approach. Quantum-chemical calculations based on coupled-cluster theory have been employed to evaluate the hyperfine parameters of nitrogen, hydrogen, and deuterium in NH2D. Subsequently, the Lamb-dip technique has been used to investigate the rotational spectrum of NH2D at high-resolution. In detail, three low-J transitions have been recorded at 86, 110, and 333 GHz with a frequency-modulation millimeter-/submillimeter-wave spectrometer. From the line profile analysis of the recorded spectra, the main terms responsible for the rotational hyperfine structure have been determined with good accuracy. Our work allows a comprehensive analysis of the rotational features of NH2D in radioastronomical spectra and a more accurate evaluation of its column density, especially in non-turbulent regions showing narrow linewidths
Rotational spectra of 1-chloro-2-fluoroethylene. I. Main isotopologues and deuterated species of the trans isomer
No full textGuided by theoretical predictions, the rotational spectra of the mono- and bideuterated species of trans-1-chloro-2-fluoroethylene, (CHCl)-Cl-35 = CDF, (CHCl)-Cl-37 = CDF, (CDCl)-Cl-35 = CHF, (CDCl)-Cl-37 = CHF, (CDCl)-Cl-35 = CDF, and (CDCl)-Cl-37 = CDF, have been recorded for the first time. Assignment of the Delta J=0 and Delta K-1=+1 bands with K-1=3,4,5,... (all isotopic species) as well as of several Delta J=+/- 1 and Delta K-1=+1 transitions (all isotopic species except (CHCl)-Cl-37 = CDF, (CDCl)-Cl-37 = CHF, and (CDCl)-Cl-37 = CDF) led to the accurate determination of the ground-state rotational constants, the quartic, and some sextic centrifugal distortion constants, as well as the nuclear quadrupole coupling constants for both Cl-35 and Cl-37 in good agreement with corresponding theoretical predictions based on high-level coupled-cluster calculations. Inconsistencies of the present spectroscopic parameters with respect to those reported earlier for the two main isotopologues, i.e., (CHCl)-Cl-35 = CHF and (CHCl)-Cl-37 = CHF, necessitated a reinvestigation of the rotational spectra for these two isotopic species. Supported by quantum chemical calculations, the previously recorded spectra are reassigned to a vibrationally excited state, while analysis of the Delta J=0 and Delta K-1=+1 as well as some Delta J=+/- 1 and Delta K-1=+1 transitions provided a revised set of spectroscopic parameters for the vibrational ground state of these two isotopic species. (c) 2006 American Institute of Physics
Cholesky Decomposition-Based Implementation of Relativistic Two-Component Coupled-Cluster Methods for Medium-Sized Molecules
No full textA Cholesky decomposition (CD)-based implementation of relativistic two-component coupled-cluster (CC) and equation-of-motion CC (EOM-CC) methods using an exact two-component Hamiltonian augmented with atomic-mean-field spin-orbit integrals (the X2CAMF scheme) is reported. The present CD-based implementation of X2CAMF-CC and EOM-CC methods employs atomic-orbital-based algorithms to avoid the construction of two-electron integrals and intermediates involving three and four virtual indices. Our CD-based implementation extends the applicability of X2CAMF-CC and EOM-CC methods to medium-sized molecules with the possibility to correlate around 1000 spinors. Benchmark calculations for uranium-containing small molecules were performed to assess the dependence of the CC results on the Cholesky threshold. A Cholesky threshold of 10-4 is shown to be sufficient to maintain chemical accuracy. Example calculations to illustrate the capability of the CD-based relativistic CC methods are reported for the bond-dissociation energy of the uranium hexafluoride molecule, UF6, with up to quadruple-ζ basis sets, and the lowest excitation energy in the solvated uranyl ion [UO22+(H2O)12]
Cholesky Decomposition in Spin-Free Dirac-Coulomb Coupled-Cluster Calculations
No full textWe present an implementation for the use of Cholesky decomposition (CD) of two-electron integrals within the spin-free Dirac-Coulomb (SFDC) scheme that enables to perform high-accuracy coupled-cluster (CC) calculations at costs almost comparable to those of their nonrelativistic counterparts. While for nonrelativistic CC calculations, atomic-orbital (AO)-based algorithms, due to their significantly reduced disk-space requirements, are the key to efficient large-scale computations, such algorithms are less advantageous in the SFDC case due to their increased computational cost in that case. Here, molecular-orbital (MO)-based algorithms exploiting the CD of the two-electron integrals allow us to reduce disk-space requirements and lead to computational cost in the CC step that is more or less the same as in the nonrelativistic case. The only remaining overhead in a CD-SFDC-CC calculation is due to the need to compute additional two-electron integrals, the somewhat higher cost of the Hartree-Fock calculation in the SFDC case, and additional cost in the transformation of the Cholesky vectors from the AO to the MO representation. However, these additional costs typically amount to less than 5-15% of the total wall time and are thus acceptable. We illustrate the efficiency of our CD scheme for SFDC-CC calculations on a series of illustrative calculations for the X(CO)4 molecules with X = Ni, Pd, Pt
Rotational spectra of 1-chloro-2-fluoroethylene. II. Equilibrium structures of the cis and trans isomer
No full textEquilibrium structures for the cis and trans isomer of 1-chloro-2-fluoroethylene are reported. The structures are obtained within a least-squares fit procedure using the available experimental ground-state rotational constants for various isotopic species of both forms. Vibrational effects were eliminated before the analysis using vibration-rotation interaction constants derived from computed quadratic and cubic force fields with the required quantum chemical calculations carried out using second-order Moller-Plesset perturbation as well as coupled-cluster (CC) theory. The semiexperimental or empirical equilibrium geometries obtained in this way agree well with the corresponding theoretical predictions obtained from CC calculations [at the CCSD(T) level] after extrapolation to the complete basis-set limit and inclusion of core-valence correlation corrections. The present results allow a detailed analysis of the geometrical differences between the two forms of 1-chloro-2-fluoroethylene. They are also compared to the structural data available for other halogenated ethylenes. (c) 2006 American Institute of Physics
Intramolecular structural parameters are key modulators of the gel-liquid transition in coarse grained simulations of DPPC and DOPC lipid bilayers
No full textThe capability of coarse-grained models based on the MARTINI mapping to reproduce the gel-liquid phase transition in saturated and unsaturated model lipids was investigated. We found that the model is able to reproduce a lower critical temperature for 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) with respect to 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC). Nonetheless, the appearance of a gel phase for DOPC is strictly dependent on the intramolecular parameters chosen to model its molecular structure. In particular, we show that the bending angle at the coarse-grained bead corresponding to the unsaturated carbon-carbon bond acts as an order parameter determining the temperature of the phase transition. Structural analysis of the molecular dynamics simulations runs evidences that in the gel phase, the packing of the lipophilic tails of DOPC assume a different conformation than in the liquid phase. In the latter phase, the DOPC geometry resembles that of the relaxed free molecule. DPPC:DOPC mixtures show a single phase transition temperature, indicating that the observation of a phase separation between the two lipids requires the simulation of systems with sizes much larger than the ones used here
Computation of NMR shieldings at the CASSCF level using gauge-including atomic orbitals and Cholesky decomposition
No full textWe present an implementation of coupled-perturbed complete active space self-consistent field (CP-CASSCF) theory for the computation of nuclear magnetic resonance chemical shifts using gauge-including atomic orbitals and Cholesky decomposed two-electron integrals. The CP-CASSCF equations are solved using a direct algorithm where the magnetic Hessian matrix-vector product is expressed in terms of one-index transformed quantities. Numerical tests on systems with up to about 1300 basis functions provide information regarding both the computational efficiency and limitations of our implementation
- …
