Abstract: While spectroscopic data on small hydrocarbons in interstellar media in combination with crossed molecular beam (CMB) experiments have provided a wealth of information on astrochemically relevant species, much of the underlying mechanistic pathways of their formation remain elusive. Therefore, in this work, the chemical reaction mechanisms of C(³P_J) + C₆H₆and C⁺(²P) + C₆H₆systems using the quantum mechanical molecular dynamics (QMMD) technique at the PBE0-D3(BJ) level of theory is investigated, mimicking a CMB experiment. Both the dynamics of the reactions as well as the electronic structure for the purpose of the reaction network are evaluated. The method is validated for the first reaction by comparison to the available experimental data. The reaction scheme for the C(³P_J) + C₆H₆system covers the literature data,<italic>e.g.</italic>the major products are the 1,2-didehydrocycloheptatrienyl radical (C₇H₅) and benzocyclopropenyl radical (C₆H₅–CH), and it reveals the existence of less common pathways for the first time. The chemistry of the C⁺(²P_J) + C₆H₆system is found to be much richer, and we have found that this is because of more exothermic reactions in this system in comparison to those in the C(³P_J) + C₆H₆system. Moreover, using the QMMD simulation, a number of reaction paths have been revealed that produce three distinct classes of reaction products with different ring sizes. All in all, at all the collision energies and orientations, the major product is the heptagon molecular ion for the ionic system. It is also revealed that the collision orientation has a dominant effect on the reaction products in both systems, while the collision energy mostly affects the charged system. These simulations both prove the applicability of this approach to simulate crossed molecular beams, and provide fundamental information on reactions relevant for the interstellar medium.

Abstract: Dielectric properties of nano-confined water are important in several areas of science, i.e., it is relevant in the dielectric double layer that exists in practically all heterogeneous fluid-based systems. Molecular dynamics simulations are used to predict the in-plane dielectric properties of confined water in planar channels of width ranging from sub-nanometer to bulk. Because of suppressed rotational degrees of freedom near the confining walls, the dipole of the water molecules tends to be aligned parallel to the walls, which results in a strongly enhanced in-plane dielectric constant (epsilon (parallel to)) reaching values of about 120 for channels with height 8 angstrom < h < 10 angstrom. With the increase in the width of the channel, we predict that epsilon (parallel to) decreases nonlinearly and reaches the bulk value for h > 70 angstrom. A stratified continuum model is proposed that reproduces the h > 10 angstrom dependence of epsilon (parallel to). For sub-nanometer height channels, abnormal behavior of epsilon (parallel to) is found with two orders of magnitude reduction of epsilon (parallel to) around h similar to 7.5 angstrom, which is attributed to the formation of a particular ice phase that exhibits long-time (similar to mu s) stable ferroelectricity. This is of particular importance for the understanding of the influence of confined water on the functioning of biological systems.

Abstract: Research progress on single layer group III monochalcogenides has been increasing rapidly owing to their interesting physics. Herein, we investigate the dynamically stable single layer forms of XBi (X = Ge, Si or Sn) using density functional theory calculations. Phonon band dispersion calculations and ab initio molecular dynamics simulations reveal the dynamical and thermal stability of the considered monolayers. Raman spectra calculations indicate the existence of 5 Raman active phonon modes, 3 of which are prominent and can be observed in possible Raman measurements. The electronic band structures of the XBi single layers were investigated with and without the effects of spin-orbit coupling (SOC). Our results show that XBi single layers show semiconducting properties with narrow band gap values without SOC. However, only single layer SiBi is an indirect band gap semiconductor, while GeBi and SnBi exhibit metallic behaviors when adding spin-orbit coupling effects. In addition, the calculated linear elastic parameters indicate the soft nature of the predicted monolayers. Moreover, our predictions for the thermoelectric properties of single layer XBi reveal that SiBi is a good thermoelectric material with increasing temperature. Overall, it is proposed that single layer XBi structures can be alternative, stable 2D single layers with varying electronic and thermoelectric properties.

Abstract: Two-dimensional materials composed of transition metal carbides and nitrides (MXenes) are poised to revolutionize energy conversion and storage. In this work, we used density functional theory (DFT) to investigate the adsorption of Mg and Na adatoms on five M2CS2 monolayers (where M = Mo, Nb, Ti, V, and Zr) for battery applications. We assessed the stability of the adatom (i.e. Na and Mg)-monolayer systems by calculating adsorption and formation energies, as well as voltages as a function of surface coverage. For instance, we found that Mo2CS2 cannot support a full layer of Na nor even a single Mg atom. Na and Mg exhibit the strongest binding on Zr2CS2, followed by Ti2CS2, Nb2CS2 and V2CS2. Using the nudged elastic band method (NEB), we computed promising diffusion barriers for both dilute and nearly full ion surface coverage cases. In the dilute ion adsorption case, a single Mg and Na atom on Ti2CS2 experience similar to 0.47 eV and similar to 0.10 eV diffusion barriers between the lowest energy sites, respectively. For a nearly full surface coverage, a Na ion moving on Ti2CS2 experiences a similar to 0.33 eV energy barrier, implying a concentration-dependent diffusion barrier. Our molecular dynamics results indicate that the three (one) layers (layer) of the Mg (Na) ion on both surfaces of Ti2CS2 remain stable at T = 300 K. While, according to voltage calculations, Zr2CS2 can store Na up to three atomic layers, our MD simulations predict that the outermost layers detach from the Zr2CS2 monolayer due to the weak interaction between Na ions and the monolayer. This suggests that MD simulations are essential to confirm the stability of an ion-electrode system – an insight that is mostly absent in previous studies.

Abstract: Inspired by the recent successful growth of Ti2C and Ti3C2 monolayers, here, we investigate the structural, electronic, and mechanical properties of functionalized Ti2C and Ti3C2 monolayers by means of density functional theory calculations. The results reveal that monolayers of Ti2C and Ti3C2 are dynamically stable metals. Phonon band dispersion calculations demonstrate that two-surface functionalization of Ti2C and Ti(3)C(2)via chalcogenides (S, Se, and Te), halides (F, Cl, Br, and I), and oxygen atoms results in dynamically stable novel functionalized monolayer materials. Electronic band dispersions and density of states calculations reveal that all functionalized monolayer structures preserve the metallic nature of both Ti2C and Ti3C2 except Ti2C-O-2, which possesses the behavior of an indirect semiconductor via full-surface oxygen passivation. In addition, it is shown that although halide passivated Ti3C2 structures are still metallic, there exist multiple Dirac-like cones around the Fermi energy level, which indicates that semi-metallic behavior can be obtained upon external effects by tuning the energy of the Dirac cones. In addition, the computed linear-elastic parameters prove that functionalization is a powerful tool in tuning the mechanical properties of stiff monolayers of bare Ti2C and Ti3C2. Our study discloses that the electronic and structural properties of Ti2C and Ti3C2 MXene monolayers are suitable for surface modification, which is highly desirable for material property engineering and device integration.

Abstract: Motivated by the recent successful synthesis of highly crystalline ultrathin BiTeCl and BiTeBr layered sheets [Debarati Hajra et al., ACS Nano, 2020, 14, 15626], herein for the first time, we carry out a comprehensive study on the structural and electronic properties of BiTeCl and BiTeBr Janus monolayers using density functional theory (DFT) calculations. Different structural and electronic parameters including the lattice constant, bond lengths, layer thickness in the z-direction, different interatomic angles, work function, charge density difference, cohesive energy and Rashba coefficients are determined to acquire a deep understanding of these monolayers. The calculations show good stability of the studied single layers. BiTeCl and BiTeBr monolayers are semiconductors with electronic bandgaps of 0.83 and 0.80 eV, respectively. The results also show that the semiconductor-metal transformation can be induced by increasing the number of layers. In addition, the engineering of the electronic structure is also studied by applying an electric field, and mechanical uniaxial and biaxial strain. The results show a significant change of the bandgaps and that an indirect-direct band-gap transition can be induced. This study highlights the positive prospect for the application of BiTeCl and BiTeBr layered sheets in novel electronic and energy conversion systems.

Abstract: In the present work, the identification of two hexagonal phases of germanium oxides (namely GeO2 and GeO) through the vibrational and optical properties is reported using density functional theory calculations. While structural optimizations show that single-layer GeO2 and GeO crystallize in 1T and buckled phases, phonon band dispersions reveal the dynamical stability of each structure. First-order off-resonant Raman spectral predictions demonstrate that each free-standing single-layer possesses characteristic peaks that are representative for the identification of the germanium oxide phase. On the other hand, electronic band dispersion analysis shows the insulating and large-gap semiconducting nature of single-layer GeO2 and GeO, respectively. Moreover, optical absorption, reflectance, and transmittance spectra obtained by means of G(0)W(0)-BSE calculations reveal the existence of tightly bound excitons in each phase, displaying strong optical absorption. Furthermore, the excitonic gaps are found to be at deep UV and visible portions of the spectrum, for GeO2 and GeO crystals, with energies of 6.24 and 3.10 eV, respectively. In addition, at the prominent excitonic resonances, single-layers display high reflectivity with a zero transmittance, which is another indication of the strong light-matter interaction inside the crystal medium.

Abstract: We calculate bubble nucleation rates in a Lennard-Jones fluid through explicit molecular dynamics simulations. Our approach-based on a recent free energy method (dubbed reweighted Jarzynski sampling), transition state theory, and a simple recrossing correction-allows us to probe a fairly wide range of rates in several superheated and cavitation regimes in a consistent manner. Rate predictions from this approach bridge disparate independent literature studies on the same model system. As such, we find that rate predictions based on classical nucleation theory, direct brute force molecular dynamics simulations, and seeding are consistent with our approach and one another. Published rates derived from forward flux sampling simulations are, however, found to be outliers. This study serves two purposes: First, we validate the reliability of common modeling techniques and extrapolation approaches on a paradigmatic problem in materials science and chemical physics. Second, we further test our highly generic recipe for rate calculations, and establish its applicability to nucleation processes.

Abstract: The added value of supervised Machine Learning (ML) methods to determine the Absolute Configuration (AC) of compounds from their Vibrational Circular Dichroism (VCD) spectra was explored. Among all ML methods considered, Random Forest (RF) and Feedforward Neural Network (FNN) yield the best performance for identification of the AC. At its best, FNN allows near-perfect AC determination, with accuracy of prediction up to 0.995, while RF combines good predictive accuracy (up to 0.940) with the ability to identify the spectral areas important for the identification of the AC. No loss in performance of either model is observed as long as the spectral sampling interval used does not exceed the spectral bandwidth. Increasing the sampling interval proves to be the best method to lower the dimensionality of the input data, thereby decreasing the computational cost associated with the training of the models.

Abstract: In the present work, a method was utilized to control the in‐situ partial hydrolysis of a diphosphonate ester in presence of a titania precursor and in function of acid content and its impact on the hybrid nanoparticles was assessed. The hydrolysis degree of organodiphosphonate ester linkers during the formation of hybrid organic‐inorganic metal oxide nanoparticles, are relatively underexplored . Quantitative solution NMR spectroscopy revealed that during the synthesis of TiO2 nanoparticles, an increase in acid concentration introduces a higher degree of partial hydrolysis of the TEPD linker into diverse acid/ester derivatives of TEPD. Increasing the HCl/Ti ratio from 1 to 3, resulted in an increase in degree of partial hydrolysis of the TEPD linker in solution from 4% to 18.8% under the here applied conditions. As a result of the difference in partial hydrolysis, the linker‐TiO2 bonding was altered. Upon subsequent drying of the colloidal TiO2 solution, different textures, at nanoscale and macroscopic scale, were obtained dependent on the HCl/Ti ratio and thus the degree of hydrolysis of TEPD. Understanding such linker‐TiO2 nanoparticle surface dynamics is crucial for making hybrid organic‐inorganic materials (i.e. (porous) metal phosphonates) employed in applications such as electronic/photonic devices, separation technology and heterogeneous catalysts.

Abstract: The interface heat transfer of two layers induced by van der Waals (vdW) contacts is theoretically investigated, based on first-principles calculations at low temperatures. The results suggest that out-of-plane acoustic phonons with low frequencies dominate the interface thermal transport due to the vdW interaction. The interface thermal conductivity is proportional to the cubic of temperature at very low temperatures, but becomes linearly proportional to temperature as temperature increases. We show that manipulating the strain alters vdW coupling, leading to increased interfacial thermal conductivity at the interface. Our findings provide valuable insights into the interface heat transport in vdW heterostructures and support further design and optimization of electronic and optoelectronic nanodevices based on vdW contacts. The heat transfer induced by van der Waals contacts is dominated by ZA phonons. The interface thermal conductivity is proportional to the cubic of temperature, but becomes linearly proportional to temperature as temperature increases.