Source author record

Fabien Brieuc

Fabien Brieuc appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.chem-ph cond-mat.stat-mech

Catalog footprint

What is connected

4works

2topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Infrared spectra at coupled cluster accuracy from neural network representations

Infrared spectroscopy is key to elucidate molecular structures, monitor reactions and observe conformational changes, while providing information on both structural and dynamical properties. This makes the accurate prediction of infrared spectra based on first-principle theories a highly desirable pursuit. Molecular dynamics simulations have proven to be a particularly powerful approach for this task, albeit requiring the computation of energies, forces and dipole moments for a large number of molecular configurations as a function of time. This explains why highly accurate first principles methods, such as coupled cluster theory, have so far been inapplicable for the prediction of fully anharmonic vibrational spectra of large systems at finite temperatures. Here, we push cutting-edge machine learning techniques forward by using neural network representations of energies, forces and in particular dipoles to predict such infrared spectra fully at "gold standard" coupled cluster accuracy as demonstrated for protonated water clusters as large as the protonated water hexamer, in its extended Zundel configuration. Furthermore, we show that this methodology can be used beyond the scope of the data considered during the development of the neural network models, allowing for the computation of finite-temperature infrared spectra of large systems inaccessible to explicit coupled cluster calculations. This substantially expands the hitherto existing limits of accuracy, speed and system size for theoretical spectroscopy and opens up a multitude of avenues for the prediction of vibrational spectra and the understanding of complex intra- and intermolecular couplings.

preprint2022arXiv

Neural Network Interaction Potentials for para-Hydrogen with Flexible Molecules

The study of molecular impurities in $para$-hydrogen ($p$H$\rm_2$) clusters is key to push forward our understanding of intra- and intermolecular interactions including their impact on the superfluid response of this bosonic quantum solvent. This includes tagging with one or very few $p$H$\rm_2$, the microsolvation regime, and matrix isolation. However, the fundamental coupling between the bosonic $p$H$\rm_2$ environment and the (ro-)vibrational motion of molecular impurities remains poorly understood. Quantum simulations can in provide the necessary atomistic insight, but very accurate descriptions of the involved interactions are required. Here, we present a data-driven approach for the generation of $impurity\cdots p$H$\rm_2$ interaction potentials based on machine learning techniques which retain the full flexibility of the impurity. We employ the well-established adiabatic hindered rotor (AHR) averaging technique to include the impact of the nuclear spin statistics on the symmetry-allowed rotational quantum numbers of $p$H$\rm_2$. Embedding this averaging procedure within the high-dimensional neural network potential (NNP) framework enables the generation of highly-accurate AHR-averaged NNPs at coupled cluster accuracy, namely CCSD(T$^*$)-F12a/aVTZcp in an automated manner. We apply this methodology to the water and protonated water molecules, as representative cases for quasi-rigid and highly-flexible molecules respectively, and obtain AHR-averaged NNPs that reliably describe the H$\rm _2$O$\cdots p$H$\rm_2$ and H$\rm _3$O$^+\cdots p$H$\rm_2$ interactions. Using path integral simulations we show for the hydronium cation that umbrella-like tunneling inversion has a strong impact on the first and second $p$H$\rm_2$ microsolvation shells. The data-driven nature of our protocol opens the door to the study of bosonic $p$H$\rm_2$ quantum solvation for a wide range of embedded impurities.

preprint2021arXiv

Transferability of machine learning potentials: Protonated water neural network potential applied to the protonated water hexamer

A previously published neural network potential for the description of protonated water clusters up to the protonated water tetramer, H$^+$(H$_2$O)$_4$, at essentially converged coupled cluster accuracy (J. Chem. Theory Comput. 16, 88 (2020)) is applied to the protonated water hexamer, H$^+$(H$_2$O)$_6$ -- a system that the neural network has never seen before. Although being in the extrapolation regime, it is shown that the potential not only allows for quantum simulations from ultra-low temperatures $\sim$ 1 K up to 300 K, but that it is able to describe the new system very accurately compared to explicit coupled cluster calculations. This transferability of the model is rationalized by the similarity of the atomic environments encountered for the larger cluster compared to the environments in the training set of the model. Compared to the interpolation regime the quality of the model is reduced by roughly one order of magnitude, but most of the difference to the coupled cluster reference comes from global shifts of the potential energy surface, while local energy fluctuations are well recovered. These results suggest that the application of neural network potentials in extrapolation regimes can provide useful results and might be more general than usually thought.

preprint2016arXiv

Zero-Point Energy Leakage in Quantum Thermal Bath Molecular Dynamics Simulations

The quantum thermal bath (QTB) has been presented as analternative to path-integral-based methods to introduce nuclear quantumeffects in molecular dynamics simulations. The method has proved to beefficient, yielding accurate results for various systems. However, the QTBmethod is prone to zero-point energy leakage (ZPEL) in highly anharmonicsystems. This is a well-known problem in methods based on classicaltrajectories where part of the energy of the high-frequency modes istransferred to the low-frequency modes leading to a wrong energydistribution. In some cases, the ZPEL can have dramatic consequences onthe properties of the system. Thus, we investigate the ZPEL by testing theQTB method on selected systems with increasing complexity in order to studythe conditions and the parameters that influence the leakage. We also analyze the consequences of the ZPEL on the structuraland vibrational properties of the system. We find that the leakage is particularly dependent on the damping coefficient and thatincreasing its value can reduce and, in some cases, completely remove the ZPEL. When using sufficiently high values for thedamping coefficient, the expected energy distribution among the vibrational modes is ensured. In this case, the QTB methodgives very encouraging results. In particular, the structural properties are well-reproduced. The dynamical properties should beregarded with caution although valuable information can still be extracted from the vibrational spectrum, even for large values ofthe damping term.