In this article, we explore the construction of Hamiltonians with long-range interactions and their corrections using the short-range behavior of the wave function. A key aspect of our investigation is the examination of the one-particle potential, kept constant in our previous work, and the effects of its optimization on the adiabatic connection. Our methodology involves the use of a parameter-dependent potential dependent on a single parameter to facilitate practical computations. We analyze the energy errors and densities in a two-electron system (harmonium) under various conditions, employing different confinement potentials and interaction parameters. The study reveals that while the mean-field potential improves the expectation value of the physical Hamiltonian, it does not necessarily improve the energy of the system within the bounds of chemical accuracy. We also delve into the impact of density variations in adiabatic connections, challenging the common assumption that a mean field improves results. Our findings indicate that as long as energy errors remain within chemical accuracy, the mean field does not significantly outperform a bare potential. This observation is attributed to the effectiveness of corrections based on the short-range behavior of the wave function, a universal characteristic that diminishes the distinction between using a mean field or not.

Continuer la lecture Partager

A general scheme for calculating ternary recombination rate constants of atomic species based on a hybrid quantum–classical nonadiabatic dynamics approach is presented and applied to the specific case of the ternary recombination of atomic ions of argon in cold argon plasmas. Rate constants are reported for both fine-structure states of the ion, and , T = 300 K, and for selected values of the reduced electric field. A thorough comparison with the literature data available for T = 300 K and a couple of close temperatures is performed with a favorable agreement achieved. It is shown that the excited ions may contribute to the formation of dimer ions, , as efficiently as the ground-state ions, , due to fast internal conversion of the electronic energy, which takes place in ternary collision complexes, .

Continuer la lecture Partager

We present the first theoretical line profile calculations of the ultraviolet spectral lines of carbon perturbed by helium using a semiclassical collision approach and high-quality ab initio potentials and electronic transition dipole moments. The temperature range is from 5000 to 8000 K. These results are important for astrophysical modelling of spectra in atmospheres of white dwarf stars showing atomic carbon in an helium atmosphere. Beyond the conventional symmetrical Lorentzian core at low He density, these lines exhibit a blue asymmetric behaviour. This blue asymmetry is a consequence of low maxima in the corresponding C–He potential energy difference curves at short internuclear distances. The collisional profiles are carefully examined and their perturber density dependence allow to understand the various line shapes of the observed carbon spectral lines in helium-rich white dwarf photosphere where the He perturber densities reach several 1021 cm−3.

Continuer la lecture Partager

Formation, distribution and behaviour of Complex Organic Molecules (COM's) in space is an important subject of research to the better understanding of the initial condition for the appearance of life on Earth. Furthermore, the study of high energy chemical processes in the interstellar medium (cosmic radiation's effect) and in solar system (solar wind's effect), is been of high interest. The aim of this work is to study astrophysical molecules trapped in interstellar ice systems under the effect of high energy radiation. These ices are characterised by being large systems, with large number of atoms. QM/MM hybrid method has become a very popular tool for molecular systems' simulations with a large number of atoms, appearing as a good compromise between accuracy and computational costs. We report the implementation of QM/MM hybrid method in the deMonNano software, using the Density Functional based Tight Binding (DFTB), an approximated DFT scheme, combined with Molecular Mechanic (MM) approach, namely Force Fields (FF) of class 1, such as OPLS-AA and AMBER-families of FFs. A complete implementation was performed using the QM/MM additive coupling scheme. In addition, the investigation of high energy chemical processes requires the explicit simulation of the electronic dynamics beyond the Born Oppenheimer approximation. As first step towards such dynamics, we will report the implementation of Real Time TD-DFTB in deMonNano, consisting in solving the Time-Dependent Schrödinger equation within the DFTB, where the electronic density matrix is propagated along time. We report a detailed introduction to new DFTB/MM and RT-TD-DFTB implementations as well as the complete study on glycine prebiotic molecule trapped in an interstellar ice. PAH interstellar systems will be also a matter of study.

Continuer la lecture Partager