Due to its protective properties of biological samples at low temperatures and under desiccation, dimethyl sulfoxide (DMSO) in aqueous solutions has been studied widely by many experimental approaches and molecular dynamics (MD) simulations. In the case of the latter, AMBER is among the most commonly used force fields for simulations of biomolecular systems; however, the parameters for DMSO published by Fox and Kollman in 1998 have only been tested for pure liquid DMSO. We have conducted an MD simulation study of DMSO in a water mixture and computed several structural and dynamical properties such as of the mean density, self-diffusion coefficient, hydrogen bonding and DMSO and water ordering. The AMBER force field of DMSO is seen to reproduce well most of the experimental properties of DMSO in water, with the mixture displaying strong and specific water ordering, as observed in experiments and multiple other MD simulations with other non-polarizable force fields.
One of the central assumptions when a particle moves through a window in microporous materials is that interaction of the diffusing particle with the silicon (Si) and aluminum (Al) atoms of the framework can be neglected, as the presence of bulkier oxygen in the host structure is thought to hinder close proximity of the diffusing particle to Si and Al. We examine this assumption, exploring the diffusion path and cross-checking the bottleneck associated with the diffusion process. Our study reveals that short-range interactions between the diffusing species and Si/Al of the host have a significant effect on the diffusion process. Guest–host interaction energy increases significantly if interaction between Si and Al atoms with the diffusing species is considered. The self-diffusion coefficient (D) decreases significantly in the linear regime, whereas in the anomalous regime, surprisingly, D increases. The increase in D is due to a decrease in the activation energy in the anomalous regime, whereas in the linear regime, activation energy increases, thus D decreases.
Graphical abstracta Interaction energies (Ea) for different LJ potential for guest–guest interactions (σgg) along the diffusion path; b correspondingdiffusivity values
The present study reports the geometries, electronic structures, growth behavior, and stabilities of neutral and ionized copper-doped germanium clusters containing 1–20 Ge atoms within the framework of linear combination of atomic orbitals density functional theory (DFT) under the spin-polarized generalized gradient approximation. It was found that Cu-capped Gen (or Cu-substituted Gen+1) and Cu-encapsulated Gen clusters mostly occur in the ground state at a particular cluster size (n). In order to explain the relative stabilities of the ground-state clusters, parameters such as the average binding energy per atom (BE), the embedding energy (EE), and the fragmentation energy (FE) of the clusters were calculated, and the resulting values are discussed. To explain the chemical stabilities of the clusters, parameters such as the energy gap between the highest occupied and the lowest unoccupied molecular orbitals (the HOMO–LUMO gap), the ionization energy (IP), the electron affinity (EA), the chemical potential (μ), the chemical hardness (η), and the polarizability were calculated, and the resulting values are also discussed. Natural atomic orbital (NAO) and natural bond orbital (NBO) analyses were also used to determine the electron-counting rule that should be applied to the most stable Ge10Cu cluster. Finally, the relevance of the calculated results to the design of Ge-based superatoms is discussed.
Figure Contributions of the valance orbitals of the Ge and Cu atom(s) to the HOMO of the ground-state icosahedral Ge10Cu cluster obtained from NBO analysis. The numbers below the clusters represent the occupancies of the HOMO orbitals
A perfectly planar Al13+ cluster (CI) and a quasi-planar Al13+ cluster (CII) have been found for the first time. Both clusters have a triangular core surrounded by a set of ten Al atoms in the form of a ring. These cationic clusters have substantial aromatic character. The planar CI cluster has local antiaromatic patches within global aromatic sea. It is doubly aromatic having both σ and π aromatic character. The quasi-planar CII cluster is also aromatic but it has more σ-delocalization.
In this work, we address the effects of molecular doping on the electronic properties of fluorinated and chlorinated silicon nanowires (SiNWs), in comparison with those corresponding to hydrogen-passivated SiNWs. Adsorption of n-type dopant molecules on hydrogenated and halogenated SiNWs and their chemisorption energies, formation energies, and electronic band gap are studied by using density functional theory calculations. The results show that there are considerable charge transfers and strong covalent interactions between the dopant molecules and the SiNWs. Moreover, the results show that the energy band gap of SiNWs changes due to chemical surface doping and it can be further tuned by surface passivation. We conclude that a molecular based ex-situ doping, where molecules are adsorbed on the surface of the SiNW, can be an alternative path to conventional doping.
Hydrogen molecule adsorption on frameworks consisting of alkaline earth metal atoms (Be, Mg, or Ca) in LTL zeolite was investigated via density functional theory. A 24T zeolite cluster model was used in this study. HOMO and LUMO energy, chemical potential, chemical hardness, electronegativity, adsorption energy, and adsorption enthalpy values were calculated. The Mg-LTL and Ca-LTL clusters were found to have much lower chemical potentials and adsorption energies than those of the Be-LTL cluster. Additionally, the calculations indicated that the Mg-LTL and Ca-LTL clusters are softer (considering their lower chemical hardness values) and more chemically reactive than the Be-LTL cluster. The calculated hydrogen adsorption enthalpies were ?14.7 and ?9.4 kJ/mol for the Mg-LTL and Ca-LTL clusters, respectively, which are significantly larger than the enthalpy of liquefaction for the hydrogen molecule. These results imply that the Mg-LTL and Ca-LTL zeolite structures are promising cryoadsorbents for hydrogen storage.
Graphical abstract Hydrogen adsorption was theoretically investigated on Be-, Ca- and Mg-LTL clusters. Ca- and Mg-LTL zeolites are potential cryoadsorbent materials for hydrogen storage.
At high temperature, silicon oxycarbide (SiCO) exhibits excellent mechanical properties and thermal stability. The incorporation of boron in SiCO results in improved performance in creep temperatures. In this work, large-scale molecular dynamics calculations were applied to obtain amorphous SiCO structures containing boron. Phase separation of C–C, B–C and Si–O was achieved for three compositions, and silicon-centered mixed-bond tetrahedrons were reproduced successfully. As the boron content increases, the boron atoms tend to form B–C and B–Si bonds in the voids, which stretches the free carbon network in some instances, causing a increase in C–C distance. Young’s modulus remains stable at high temperature for the high-carbon case, which indicates that the free carbon network plays a critical role in the structural and thermal stability of SiBCO.
Graphical Abstract Three major typical structures in the cooling down process for silicon boron oxycarbide (Si5BC2O8). Bonds: red Si–O, blue Si–C, black C–C, green B–C, purple Si–B
Efficient n-type dye-sensitized solar cells are known since the seminal work of O’Reagan and Grätzel in 1991. However, highly efficient p-type dye-sensitized solar cells have not been developed so far. This hinders the construction of tandem dye-sensitized solar cells, which can surpass the performance of n-type devices. Within this work, we investigate if a temporary coordination of transition metal-based redox mediators at a sensitizer can increase the efficiency of p-type dye-sensitized solar cells. Based on a computational screening, diverse Cu, Ni, and Co redox mediators were selected to construct p-type dye-sensitized solar cells. Unfortunately, the efficiency of the investigated devices does not surpass analogous cells with iodide-triiodide as redox mediator. While Ni and Cu complexes might be reduced to Ni(0) and Cu(0), respectively, the investigated Co-complex quenches the excited state efficiently. As a result, the necessary electron injection from the semiconductor is too slow, which hinders the construction of a highly efficient p-type dye-sensitized solar cell.
Graphical Abstract Comparison of the mode of action of p-type dye-sensitized solar cells. While top shows the traditional one, bottom shows the investigated devices where a temporary link between dye and redox mediator plays a crucial role.
An all-atom force field consistent with the general AMBER force field (GAFF) format for poly(ethylene glycol) dimethyl ether (diglyme or G2) was developed by fitting to experimental liquid densities and dielectric constants. Not surprisingly, the new force field gives excellent agreement with experimental liquid phase densities and dielectric constants over a wide temperature range. Other dynamic and thermodynamic properties of liquid G2 such as its self-diffusion coefficient, shear viscosity, and vaporization enthalpy were also calculated and compared to experimental data. For all of the properties studied, the performance of the proposed new force field is better than that of the standard GAFF force field. The force field parameters were transferred to model two other poly(ethylene glycol) ethers: monoglyme (G1) and tetraglyme (G4). The predictive ability of the modified force field for G1 and G4 was significantly better than that of the original GAFF force field. The proposed force field provides an alternative option for the simulation of mixtures containing glymes using GAFF-compatible force fields, particularly for electrochemical applications. The accuracy of a previously published force field based on the OPLS-AA format and the accuracies of two modified versions of that force field were also examined for G1, G2, and G4. It was found that the original OPLS-AA force field is superior to the modified versions of it, and that it has a similar accuracy to the proposed new GAFF-compatible force field.
Density functional theory (DFT) calculations were carried out to explore the geometric, spectroscopic, and electronic properties of three magic silver clusters Agn (n?=?8, 18, and 20) in detail. The computed results show that the global minima of these clusters are compact, near-spherical structures, while other low-lying isomers exhibit oblate or prolate shapes. Vertical ionization energies for the low-lying isomers were also computed and assigned with respect to available experimental values. Although several isomers were predicted to have similar energies, their electronic and vibrational signatures were quite distinctive, meaning that they could be used as fingerprint signals to distinguish between isomers. In addition, the electronic structures of these systems were explored using the phenomenological shell model. Calculations for the coinage metal clusters M20 (M?=?Cu, Ag, Au) indicated that the structures and properties of the Ag cluster are similar to those of the Cu cluster in that both Cu20 and Ag20 prefer a compact structure whereas Au20 prefers to adopt a tetrahedral form.
Cinnamaldehyde (CAL) belongs to the group of aromatic α,β-unsaturated aldehydes; the selective hydrogenation of CAL plays an important role in the fine chemical and pharmaceutical industries. Using Ptn clusters as catalytic models, we studied the selective hydrogenation reaction mechanism for CAL catalyzed by Ptn (n?=?6, 10, 14, 18) clusters by means of B3LYP in density functional theory at the 6-31+?G(d) level (the LanL2DZ extra basis set was used for the Pt atom). The rationality of the transition state was proved by vibration frequency analysis and intrinsic reaction coordinate computation. Moreover, atoms in molecules theory and nature bond orbital theory were applied to discuss the interaction among orbitals and the bonding characteristics. The results indicate that three kinds of products, namely 3-phenylpropyl aldehyde, 3-phenyl allyl alcohol and cinnamyl alcohol, are produced in the selective hydrogenation reaction catalyzed by Ptn clusters; each pathway possesses two reaction channels. Ptn clusters are more likely to catalyze the activation and hydrogenation of the C?=?O bond in CAL molecules, eventually producing cinnamic alcohol, which proves that Ptn clusters have a strong reaction selectivity to catalyze CAL. The reaction selectivity of the catalyzer cluster is closely related to the size of the Ptn cluster, with Pt14 clusters having the greatest reaction selectivity.
Graphical Abstract The reaction mechanism for the selective hydrogenation reaction ofcinnamaldehyde catalyzed by Ptn clusters was studied by densityfunctional theory. The reactionselectivity of cluster catalyzer was concluded to be closely related to the size of Ptn clusters, with Pt14 clusters having the greatest reaction selectivity
The present paper reports the analysis of surface decoration on the structural, electronic, and optical properties of (n,0) ZnO nanotubes, performed by means of a density function theory based ab-initio approach. Fe functionalization induced buckling in ZnO nanotubes affects its electronic and optical properties. Increase in Fe functionalization leads to better stability of ZnO nanotube and shows enhanced metallic character. The possibility of its use in optoelectronics has been analyzed in terms of dielectric constant, absorption coefficient, and refractive index. In another observation, the high sensitivity of the HCN molecule for the Fe-incorporated ZnO nanotube suggests it as a potential gas sensor.
Utilizing first-principles calculations, we studied the electronic and optical properties of C24, C12X6Y6, and X12Y12 fullerenes (X?=?B, Al; Y?=?N, P). These fullerenes are energetically stable, as demonstrated by their negative cohesive energies. The energy gap of C24 may be tuned by doping, and the B12N12 fullerene was found to have the largest energy gap. All of the fullerenes had finite optical gaps, suggesting that they are optical semiconductors, and they strongly absorb UV radiation, so they could be used in UV light protection devices. They could also be used in solar cells and LEDs due to their low reflectivities.
The catalytic coupling reaction mechanism of the transformation from 4-nitrobenzenethiol (4-NBT) to 4,4′-dimercaptoazobenzene (4,4′-DMAB) on a silver cluster was studied by density functional theory. Reactants, intermediates, transition states and products were optimized with the B3LYP method using the 6-311 + G(d,p) basis set (Ag using the pseudo potential basis set of LanL2DZ). Transition states and intermediates were confirmed by the corresponding vibration analysis and intrinsic reaction coordinates (IRC). Consistent with literature reports, the key point of the transformation from 4-NBT absorbed on the surface of Ag5 clusters to 4,4′-DMAB is the elimination of two O atoms on the amino group. Meanwhile, the catalytic coupling reaction of 4-nitrobenzenethiol on a silver cluster is easy to carry out under irradiation. The possibility of “inter system channeling” (ISC) between different potential energy surfaces in the coupling reaction of 4-NBT is further discussed. The irradiation has an auxiliary catalytic effect on the coupling reaction. Our research results can explain the observed experimental phenomena.
Graphical abstract Catalytic coupling reaction mechanism of the transformation from 4-nitrothiophenol (4-NBT) to 4,4′-dimercaptoazobenzene (4,4′-DMAB) on silver clusters studied by density functional theory
In this article, we explore the capacity of formed Schiff base complexes to trap metal atoms or ions, using their aromatic ends. The intrinsic geometry of each complex defines the process of substitution. Two cases were studied; one involving a trans Schiff base complex and the other considering how a salen ligand, with nickel systems traps chromium. We also assessed the nature of the new bonds and the frontier molecular orbitals.
Seven models that related the features of molecular surface electrostatic potentials (ESPs) above the bond midpoints and rings, statistical parameters of ESPs to the experimental impact sensitivities h50 of eight strained cyclic explosives with the C–NO2 bonds were theoretically predicted at the DFT-B3LYP/6-311++G** level. One of the models was used to investigate the changes of h50 for the nitrocyclohydrocarbon frameworks in the H-bonded complexes of HF with nitrocyclopropane, nitrocyclobutane, nitrocyclopentane, and nitrocyclohexane. The results show that the correlation coefficients of the obtained models are small. When adding the effect of ring strain, the value of correlation coefficient is increased. According to the calculated h50, the sensitivities in the frameworks are increased after hydrogen bonding. As a global feature of molecules, surface electrostatic potential is more available to judge the sensitivity change than the trigger bond dissociation energy or ring strain energy in H-bonded complex.
Graphical Abstract A theoretical prediction of the relationships between the impact sensitivity and electrostatic potential in strained cyclic explosive and application to H-bonded complex of nitrocyclohydrocarbon?
In this article, we explore, both theoretically and experimentally, the general reactivity of alkyl hydrogeno-phenylphosphinates with alcohols. We show that alcohol molecules act exclusively as nucleophilic species, and add to alkyl hydrogeno-phenylphosphinates, leading to pentacoordinated intermediates. These intermediates are shown to subsequently competitively undergo alcohol eliminations and/or Berry pseudorotations. This offers several possible routes for racemizations and/or alcohol exchange reactions. Transition standard Gibbs free energies predicted from DFT calculations for the overall alcohol exchange mechanism are shown to be compatible with those experimentally measured in case ethanol reacts with ethyl hydrogeno-phenylphosphinate (134.5～136.0 kJ mol?1 at 78 °C).
Bond critical points (BCPs) in the quantum theory of atoms in molecules (QTAIM) are shown to be a consequence of the molecular topology, symmetry, and the Poincaré-Hopf relationship, which defines the numbers of critical points of different types in a scalar field. BCPs can be induced by a polarizing field or by addition of a single non-bonded atom to a molecule. BCPs and their associated bond paths are therefore suggested not to be a suitable means of identifying chemical bonds, or even attractive intermolecular interactions.
Vitamin C is one of the most abundant exogenous antioxidants in the cell, and it is of the utmost importance to elucidate its mechanism of action against radicals. In this study, the reactivity of vitamin C toward OH and \( {HO}_2/{O}_2^{-} \) radicals in aqueous medium was analyzed by ab initio molecular dynamics using CPMD code. The simulations led to results similar to those of static studies or experiments for the pair of \( {HO}_2/{O}_2^{-} \) radicals but bring new insights for the reactivity with hydroxyl radical: the reaction takes place before the formation of an adduct and consists of two steps: first an electron is transferred to hydroxyl radical and then the ascorbyl radical loses a proton.
