2024-10-30 | | Total: 12
A striking example of the need to accurately capture states of double-excitation character in molecules is seen in predicting photo-induced dynamics in small polyenes. Due to the coupling of electronic and nuclear motions, the dark 2$^1$Ag state, known to have double-excitation character, can be reached after an initial photo-excitation to the bright $^1$Bu state via crossings of their potential energy surfaces. However, the shapes of the surfaces are so poorly captured by most electronic structure methods, including adiabatic time-dependent density functional theory (TDDFT), complete active space self-consistent field, and equation of motion coupled cluster singles and doubles, that the crossing is missed or substantially mislocated. We demonstrate that the frequency-dependent kernel of dressed TDDFT beyond Tamm-Dancoff successfully captures the curve-crossing, providing an energy surface close to the highly accurate but expensive $\delta$-CR-EOMCC(2,3) benchmark reference. This, along with its accurate prediction of the excitation character of the state makes dressed TDDFT a practical and accurate route to electronic structure quantities needed in modeling ultrafast dynamics in molecules.
Efficient hybrid DFT simulations of solid state materials would be extremely beneficial for computational chemistry and materials science, but is presently bottlenecked by difficulties in computing Hartree-Fock (HF) exchange with plane wave orbital bases. We present a GPU-accelerated, Gaussian orbital based integral algorithm for systems with periodic boundary conditions, which takes advantage of Ewald summation to efficiently compute electrostatic interactions. We have implemented this approach into the TeraChem software package within the $\Gamma$ point approximation, enabling simulation of unit cells with hundreds or thousands of atoms at the HF or hybrid DFT level on a single GPU card. Our implementation readily parallelizes over multiple GPUs and paves the road to accurate simulation of the properties and dynamics of extended materials in both the ground and excited states.
Room temperature single photon sources (SPS) are crucial for developing the next generation quantum technologies. Quantum dots (QDs), recently, have been reported as promising materials as SPS at room temperature. By optimizing the single particle optical properties of a series of water-soluble, $CdTe_{x}Se_{1-x}$, here we provide an efficient SPS with increased single photon purity. The data revealed that second order photon correlation, $g^2(0)$ value decreases substantially from 0.21 in CdTe to 0.02 in $CdTe_{0.25}Se_{0.75}$ QDs. They also exhibited deterministic emissions with an increase in ON time exceeding 95% of the total time. This was accompanied by an increased photon count rate, substantially reduced blinking events, and extended single particle ON-time. The increased single photon emission in $CdTe_{x}Se_{1-x}$ is attributed to very fast electron trapping to dense trap states, which suppresses the multiexciton recombination.
Experiments show that light-matter strong coupling affects chemical properties, though the underlying mechanism remains unclear. We present an ab initio quantum electrodynamics coupled cluster method for collective light-matter strong coupling to investigate the interplay of local molecular details with the collective polaritonic excitation. The model is used to compute the potential energy surfaces of the argon dimer, which provides a general prototype for excimers, and analyze the ground and excited state vibrational levels. Potential energy surfaces are significantly modified for large light-matter couplings (single-molecule regime). On the other hand, in the collective regime (small light-matter coupling, large number of molecules), the ground state potential energy surface and the first vibrational levels of the excited state do not change significantly. However, collective strong coupling produces an abrupt transition in the vibrational landscape of the excimer, causing higher vibrational levels to behave similarly to ground state vibrations. Beyond a critical coupling strength, the excimer formation is thus inhibited for both the single-molecule and collective coupling regimes via different mechanisms. We propose this is a general feature of excimers under strong coupling, offering an additional aspect of polaritonic chemistry.
The accurate simulation of complex biochemical phenomena has historically been hampered by the computational requirements of high-fidelity molecular-modeling techniques. Quantum mechanical methods, such as ab initio wave-function (WF) theory, deliver the desired accuracy, but have impractical scaling for modelling biosystems with thousands of atoms. Combining molecular fragmentation with MP2 perturbation theory, this study presents an innovative approach that enables biomolecular-scale ab initio molecular dynamics (AIMD) simulations at WF theory level. Leveraging the resolution-of-the-identity approximation for Hartree-Fock and MP2 gradients, our approach eliminates computationally intensive four-center integrals and their gradients, while achieving near-peak performance on modern GPU architectures. The introduction of asynchronous time steps minimizes time step latency, overlapping computational phases and effectively mitigating load imbalances. Utilizing up to 9,400 nodes of Frontier and achieving 59% (1006.7 PFLOP/s) of its double-precision floating-point peak, our method enables us to break the million-electron and 1 EFLOP/s barriers for AIMD simulations with quantum accuracy.
During the dielectric breakdown process of thin solid-state nanopores, the application of high voltages may cause the formation of multi-nanopores on one chip, which number and sizes are important for their applications. Here, simulations were conducted to mimic the investigation of in situ nanopore detection with scanning ion conductance microscopy (SICM). Results show that SICM can provide accurate nanopore location and relative pore size. Detection resolution is influenced by the dimensions of the applied probe and separation between the probe and membranes, which can be enhanced under large voltages or a concentration gradient.
We investigate how growth of an oxide film will influence the deuterium recycling properties of Li on plasma facing surfaces. Lithium films on the walls or plasma-facing material surfaces of a fusion vacuum vessel improves plasma performance in part by removing residual impurity atoms from the plasma. Oxygen atoms, from residual water vapor or eroded oxide surfaces, are a prominent ubiquitous impurity that can create a layer of oxide on lithium films. The results of the work herein also indicate a possibility for determining the thickness of the oxide layer on such Li films by measuring the reflection probability of incident deuterium particles from the surface.
Optimizing the synthesis of zeolites and exploring novel frameworks offer pivotal opportunities and challenges in materials design. While inverse design proves highly effective for simpler crystals, its application to intricate structures like zeolites poses severe challenges. Here, we introduce an innovative inverse design workflow tailored to efficiently reproduce target zeolite frameworks in a binary coarse-grained model using enhanced sampling molecular dynamics simulations. This workflow integrates an evolutionary parameter optimization strategy with a variant of the seeding approach. Using this method, we successfully reproduce Z1 and SGT zeolites, and Type-I clathrates, find new optimal parameters for known phases, such as the SOD and CFI, and even discover novel frameworks, such as Z5. This is done within a simple coarse-grained model for a tetrahedra-forming component and a structure-directing agent. Our methodology not only enables the screening of synthesis protocols but also facilitates the discovery of hypothetical zeolites.
In this work, we perform a systematic computer simulation study of ice premelting, and explore the thickness and structure of quasi-liquid layers formed at the interface of ice with substrates of different hydrophilicity. Our study shows that interfacial premelting occurs on neutral substrates of whatever hydrophilicity, forming films of limited thickness for substrates with contact angles larger than ca. 50$^\circ$ but exhibiting complete interfacial premelting at smaller contact angles. Contrary to most experimental studies, we focus not only on the premelting behavior with temperature, but also with pressure, which is a matter of relevance in important situations such as ice friction. Our study is guided within a rigorous framework of surface thermodynamics, demonstrating that the premelting film structure is a function of a single thermodynamic variable. By this token we are able to relate properties measured along an isobar, with premelting films at arbitrary temperature and pressure. We also find that adhesion strength in atomically smooth surfaces is one to two orders of magnitude larger than those found in experiments, and conjecture that the reason is substrate roughness and the presence of organic adsorbents. Our theoretical framework also allows us to exploit our results on interfacial premelting in order to gain insight into heterogeneous ice nucleation. We show that apolar smooth substrates of whatever hydrophilicity are unlikely nucleators, and that too large hydrophilicity conspires also against ice nucleation. Furthermore, we exploit statistical-thermodynamic framework to shed light into the nature of the surface intermolecular forces promoting interfacial premelting, and provide a model to predict quasi-liquid layer thickness as a function of the substrate's hydrophilicity with great potential applications in fields ranging from earth sciences to aircraft engineering.
Constructing transferable descriptors for conformation representation of molecular and biological systems finds numerous applications in drug discovery, learning-based molecular dynamics, and protein mechanism analysis. Geometric graph neural networks (Geom-GNNs) with all-atom information have transformed atomistic simulations by serving as a general learnable geometric descriptors for downstream tasks including prediction of interatomic potential and molecular properties. However, common practices involve supervising Geom-GNNs on specific downstream tasks, which suffer from the lack of high-quality data and inaccurate labels leading to poor generalization and performance degradation on out-of-distribution (OOD) scenarios. In this work, we explored the possibility of using pre-trained Geom-GNNs as transferable and highly effective geometric descriptors for improved generalization. To explore their representation power, we studied the scaling behaviors of Geom-GNNs under self-supervised pre-training, supervised and unsupervised learning setups. We find that the expressive power of different architectures can differ on the pre-training task. Interestingly, Geom-GNNs do not follow the power-law scaling on the pre-training task, and universally lack predictable scaling behavior on the supervised tasks with quantum chemical labels important for screening and design of novel molecules. More importantly, we demonstrate how all-atom graph embedding can be organically combined with other neural architectures to enhance the expressive power. Meanwhile, the low-dimensional projection of the latent space shows excellent agreement with conventional geometrical descriptors.
The classical nucleation theory (CNT) and its modified versions provide a convenient framework for describing the nucleation process under the capillary approximation. However, these models often predict nucleation rates that depart significantly from simulation results, even for a simple Lennard-Jones fluid. This large discrepancy is likely due to the inaccurate estimation of the driving force for nucleation, which most traditional models estimate within the ideal solution approximation. In this study, we address this issue by directly calculating the driving force for nucleation using equations of state (EOS) and integrating this approach into the calculation of nucleation rates within the framework of CNT and its modified model. We apply this method to examine the condensation of a Lennard-Jones fluid and compare the resulting nucleation rates with molecular dynamics (MD) simulation data. Our results demonstrate that at relatively low supersaturation, where the capillary approximation is reasonable, our thermodynamic models exhibit excellent agreement with MD results, significantly outperforming traditional models. At moderate and high supersaturation, our approach continues to show a reasonable agreement with MD results. Furthermore, when comparing the results obtained by using different EOS, we find that more precise EOS generally yield better agreement with MD data.
Templated copolymerization, in which information stored in the sequence of a heteropolymer template is copied into another polymer product, is the mechanism behind all known methods of genetic information transfer. A key aspect of templated copolymerization is the eventual detachment of the product from the template. A second key feature of natural biochemical systems is that the template-binding free energies of both correctly-matched and incorrect monomers are heterogeneous. Previous work has considered the thermodynamic consequences of detachment, and the consequences of heterogeniety for polymerisation speed and accuracy, but the interplay of both separation and heterogeniety remains unexplored. In this work, we investigate a minimal model of templated copying that simultaneously incorporates both detachment from behind the leading edge of the growing copy and heterogeneous interactions. We first extend existing coarse-graining methods for models of polymerisation to allow for heterogeneous interactions. We then show that heterogeneous copying systems with explicit detachment do not exhibit the subdiffusive behaviour observed in the absence of detachment when near equilibrium. Next, we show that heterogeneity in correct monomer interactions tends to result in slower, less accurate copying, while heterogeneity in incorrect monomer interactions tends to result in faster, more accurate copying, due to an increased roughness in the free energy landscape of either correct or incorrect monomer pairs. Finally, we show that heterogeneity can improve on known thermodynamic efficiencies of homogeneous copying, but these increased thermodynamic efficiencies do not always translate to increased efficiencies of information transfer.