No Arabic abstract
The method of increments (MoI) allows one to successfully calculate cohesive energies of bulk materials with high accuracy, but it encounters difficulties when calculating whole dissociation curves. The reason is that its standard formalism is based on a single Hartree-Fock (HF) configuration whose orbitals are localized and used for the many-body expansion. Therefore, in those situations where HF does not allow a size-consistent description of the dissociation, the MoI cannot yield proper results either. Herein we address the problem by employing a size-consistent multiconfigurational reference for the MoI formalism. This leads to a matrix equation where a coupling derived by the reference itself is employed. In principle, such approach allows one to evaluate approximate values for the ground as well as excited states energies. While the latter are accurate close to the avoided crossing only, the ground state results are very promising for the whole dissociation curve, as shown by the comparison with density matrix renormalization group (DMRG) benchmarks. We tested this two-state constant-coupling (TSCC)-MoI on beryllium rings of different sizes and studied the error introduced by the constant coupling.
We present an implementation of a fully self-consistent finite temperature second order Greens function perturbation theory (GF2) within the diagrammatic Monte Carlo framework. In contrast to the previous implementations of stochastic GF2 ({it J. Chem. Phys.},{bf 151}, 044144 (2019)), the current self-consistent stochastic GF2 does not introduce a systematic bias of the resulting electronic energies. Instead, the introduced implementation accounts for the stochastic errors appearing during the solution of the Dyson equation. We present an extensive discussion of the error handling necessary in a self-consistent procedure resulting in dressed Greens function lines. We test our method on a series of simple molecular examples.
We present an approximate scheme for analytical gradients and nonadiabatic couplings for calculating state-average density matrix renormalization group self-consistent-field wavefunction. Our formalism follows closely the state-average complete active space self-consistent-field (SA-CASSCF) emph{ansatz}, which employs a Lagrangian, and the corresponding Lagrange multipliers are obtained from a solution of the coupled-perturbed CASSCF (CP-CASSCF) equations. We introduce a definition of the matrix product state (MPS) Lagrange multipliers based on a single-site tensor in a mixed-canonical form of the MPS, such that a sweep procedure is avoided in the solution of the CP-CASSCF equations. We apply our implementation to the optimization of a conical intersection in 1,2-dioxetanone, where we are able to fully reproduce the SA-CASSCF result up to arbitrary accuracy.
We present the first implementation of a density matrix renormalization group algorithm embedded in an environment described by density functional theory. The frozen density embedding scheme is used with a freeze-and-thaw strategy for a self-consistent polarization of the orbital-optimized wavefunction and the environmental densities with respect to each other.
The electrical resistivity of the quasi-1D organic superconductor (TMTSF)2PF6 was recently measured at low temperature from the critical pressure needed to suppress the spin-density-wave state up to a pressure where superconductivity has almost disappeared. This data revealed a direct correlation between the onset of superconductivity at Tc and the strength of a non-Fermi-liquid linear term in the normal-state resistivity, going as r(T) = r0 + AT + BT2 at low temperature, so that A goes to 0 as Tc goes to 0. Here we show that the contribution of low-frequency antiferromagnetic fluctuations to the spin-lattice relaxation rate is also correlated with this non-Fermi-liquid term AT in the resistivity. These correlations suggest that anomalous scattering and pairing have a common origin, both rooted in the low-frequency antiferromagnetic fluctuations measured by NMR. A similar situation may also prevail in the recently-discovered iron-pnictide superconductors.
We present a matrix-product state (MPS)-based quadratically convergent density-matrix renormalization group self-consistent-field (DMRG-SCF) approach. Following a proposal by Werner and Knowles (JCP 82, 5053, (1985)), our DMRG-SCF algorithm is based on a direct minimization of an energy expression which is correct to second-order with respect to changes in the molecular orbital basis. We exploit a simultaneous optimization of the MPS wave function and molecular orbitals in order to achieve quadratic convergence. In contrast to previously reported (augmented Hessian) Newton-Raphson and super-configuration-interaction algorithms for DMRG-SCF, energy convergence beyond a quadratic scaling is possible in our ansatz. Discarding the set of redundant active-active orbital rotations, the DMRG-SCF energy converges typically within two to four cycles of the self-consistent procedure