No Arabic abstract
We present a generally-applicable computational framework for the efficient and accurate characterization of molecular structural patterns and acid properties in explicit solvent using H$_2$O$_2$ and CH$_3$SO$_3$H in phenol as an example. In order to address the challenges posed by the complexity of the problem, we resort to a set of data-driven methods and enhanced sampling algorithms. The synergistic application of these techniques makes the first-principle estimation of the chemical properties feasible without renouncing to the use of explicit solvation, involving extensive statistical sampling. Ensembles of neural network potentials are trained on a set of configurations carefully selected out of preliminary simulations performed at a low-cost density-functional tight-binding (DFTB) level. The energy and forces of these configurations are then recomputed at the hybrid density functional theory (DFT) level and used to train the neural networks. The stability of the NN model is enhanced by using DFTB energetics as a baseline, but the efficiency of the direct NN (i.e., baseline-free) is exploited via a multiple-time step integrator. The neural network potentials are combined with enhanced sampling techniques, such as replica exchange and metadynamics, and used to characterize the relevant protonated species and dominant non-covalent interactions in the mixture, also considering nuclear quantum effects.
I use first principles calculations to investigate the thermal conductivity of $beta$-In$_2$O$_3$ and compare the results with that of $alpha$-Al$_2$O$_3$, $beta$-Ga$_2$O$_3$, and KTaO$_3$. The calculated thermal conductivity of $beta$-In$_2$O$_3$ agrees well with the experimental data obtain recently, which found that the low-temperature thermal conductivity in this material can reach values above 1000 W/mK. I find that the calculated thermal conductivity of $beta$-Ga$_2$O$_3$ is larger than that of $beta$-In$_2$O$_3$ at all temperatures, which implies that $beta$-Ga$_2$O$_3$ should also exhibit high values of thermal conductivity at low temperatures. The thermal conductivity of KTaO$_3$ calculated ignoring the temperature-dependent phonon softening of low-frequency modes give high-temperature values similar that of $beta$-Ga$_2$O$_3$. However, the calculated thermal conductivity of KTaO$_3$ does not increase as steeply as that of the binary compounds at low temperatures, which results in KTaO$_3$ having the lowest low-temperature thermal conductivity despite having acoustic phonon velocities larger than that of $beta$-Ga$_2$O$_3$ and $beta$-In$_2$O$_3$. I attribute this to the fact that the acoustic phonon velocities at low frequencies in KTaO$_3$ is less uniformly distributed because its acoustic phonon branches are more dispersive compared to the binary oxides, which causes enhanced momentum loss even during the normal phonon-phonon scattering processes. I also calculate thermal diffusivity using the theoretically obtained thermal conductivity and heat capacity and find that all four materials exhibit the expected $T^{-1}$ behavior at high temperatures. Additionally, the calculated ratio of the average phonon scattering time to Planckian time is larger than the lower bound of 1 that has been observed empirically in numerous other materials.
OH 231.8+4.2, a bipolar outflow around a Mira-type variable star, displays a unique molecular richness amongst circumstellar envelopes (CSEs) around O-rich AGB and post-AGB stars. We report line observations of the HCO+ and H13CO+ molecular ions and the first detection of SO+, N2H+, and (tentatively) H3O+ in this source. SO+ and H3O+ have not been detected before in CSEs around evolved stars. These data have been obtained as part of a full mm-wave and far-IR spectral line survey carried out with the IRAM 30 m radio telescope and with Herschel/HIFI. Except for H3O+, all the molecular ions detected in this work display emission lines with broad profiles (FWHM 50-90 km/s), which indicates that these ions are abundant in the fast bipolar outflow of OH 231.8. The narrow profile (FWHM 14 km/s) and high critical densities (>1e6cm-3 ) of the H3O+ transitions observed are consistent with this ion arising from denser, inner (and presumably warmer) layers of the fossil remnant of the slow AGB CSE at the core of the nebula. From rotational diagram analysis, we deduce excitation temperatures of Tex 10-20 K for all ions except for H3O+, which is most consistent with Tex 100 K. Although uncertain, the higher excitation temperature suspected for H3O+ is similar to that recently found for H2O and a few other molecules, which selectively trace a previously unidentified, warm nebular component.The column densities of the molecular ions reported here are in the range Ntot [1-8]x1e13 cm-2, leading to beam-averaged fractional abundances relative to H2 of X(HCO+) 1e-8, X(H13CO+) 2e-9, X(SO+) 4e-9, X(N2H+) 2e-9, and X(H3O+) 7e-9 cm-2. We have performed chemical kinetics models to investigate the formation of these ions in OH 231.8 as the result of standard gas phase reactions initiated by cosmic-ray and UV-photon ionization. (abridged).
We calculated reaction rate constants including atom tunneling of the reaction of dihydrogen with the hydroxy radical down to a temperature of 50 K. Instanton theory and canonical variational theory with microcanonical optimized multidimensional tunneling (CVT/$mu$OMT) were applied using a fitted potential energy surface [J. Chem. Phys. 138, 154301 (2013)]. All possible protium/deuterium isotopologues were considered. Atom tunneling increases at about 250 K (200 K for deuterium transfer). Even at 50 K the rate constants of all isotopologues remain in the interval $ 4 cdot 10^{-20}$ to $4 cdot 10^{-17}$ cm$^3$ s$^{-1}$ , demonstrating that even deuterat
In a previous study of the L1157 B1 shocked cavity, a comparison between NH$_3$(1$_0$-$0_0$) and H$_2$O(1$_{rm 10}$--1$_{rm 01}$) transitions showed a striking difference in the profiles, with H$_2$O emitting at definitely higher velocities. This behaviour was explained as a result of the high-temperature gas-phase chemistry occurring in the postshock gas in the B1 cavity of this outflow. If the differences in behaviour between ammonia and water are indeed a consequence of the high gas temperatures reached during the passage of a shock, then one should find such differences to be ubiquitous among chemically rich outflows. In order to determine whether the difference in profiles observed between NH$_3$ and H$_2$O is unique to L1157 or a common characteristic of chemically rich outflows, we have performed Herschel-HIFI observations of the NH$_3$(1$_0$-0$_0$) line at 572.5 GHz in a sample of 8 bright low-mass outflow spots already observed in the H$_2$O(1$_{rm 10}$--1$_{rm 01}$) line within the WISH KP. We detected the ammonia emission at high-velocities at most of the outflows positions. In all cases, the water emission reaches higher velocities than NH$_3$, proving that this behaviour is not exclusive of the L1157-B1 position. Comparisons with a gas-grain chemical and shock model confirms, for this larger sample, that the behaviour of ammonia is determined principally by the temperature of the gas.
We report the crystal growth and structural and magnetic properties of quasi two-dimensional $S=1/2$ quantum magnet Cu[C$_6$H$_2$(COO)$_4$][H$_3$N-(CH$_2$)$_2$-NH$_3$]$cdot$3H$_2$O. It is found to crystallize in a monoclinic structure with space group $C2/m$. The CuO$_4$ plaquettes are connected into a two-dimensional framework in the $ab$-plane through the anions of [C$_6$H$_2$(COO)$_4$]$^{4-}$ (pyromellitic acid). The [H$_3$N-(CH$_2$)$_2$-NH$_3$]$^{2+}$$cdot$3H$_2$O groups are located between the layers and provide a weak interlayer connection via hydrogen (H...O) bonds. The temperature dependent magnetic susceptibility is well described by $S=1/2$ frustrated square lattice ($J_1-J_2$) model with nearest-neighbor interaction $J_1/k_{rm B} simeq 5.35$ K and next-nearest-neighbor interaction $J_2/k_{rm B} simeq -0.01$ K. Even, our analysis using frustrated rectangular lattice ($J_{1a,b}-J_2$) model confirms almost isotropic nearest-neighbour interactions ($J_{rm 1a}/k_{rm B} simeq 5.31$ K and $J_{rm 1b}/k_{rm B} simeq 5.38$ K) in the $ab$-plane and $J_2/k_{rm B}simeq-0.24$ K. Further, the isothermal magnetization at $T=1.9$ K is also well described by a non-frustrated square lattice model with $J_1/k_{rm B} simeq 5.2$ K. Based on the $J_2/J_1$ ratio, the compound can be placed in the N{e}el antiferromagnetic state of the $J_1 - J_2$ phase diagram. No signature of magnetic long-range-order was detected down to 2 K.