Journal of Chemical Physics最新文献

It is still critical to prepare a high-quality absorber layer for high-performance Cu2ZnSnSe4 (CZTSe) multi-component thin film solar cell. The gas pressure during the selenization process is commonly referred to as the pressure of inert gas in the tube furnace, while the exact selenium partial pressure is difficult to be controlled. Therefore, the grain growth under different selenium partial pressures cannot be made clear, and the film quality cannot be controlled as well. In this work, we use a sealed quartz tube as the selenization vessel, which can provide a relatively high and controllable selenium partial pressure during the selenization process. To further tailor the grain growth, lithium doping is also utilized. We find that lithium can greatly promote the growth of CZTSe films as the selenium partial pressure is controlled near the selenium saturation vapor pressure. Combined with ALD-Al2O3, the crystallization quality of CZTSe absorber films is significantly enhanced and the efficiency of CZTSe solar cells achieved a significant improvement. This work clarifies the effect of controllable Se pressure on CZTSe film growth and can lead to better results in CZTSe and other multi-compound thin film solar cells.

Path-like collective variables (CVs) can be very effective for accurately modeling complex biomolecular processes in molecular dynamics simulations. Recently, we have introduced DeepLNE (deep-locally non-linear-embedding), a machine learning-based path-like CV that provides a progression variable s along the path as a non-linear combination of several descriptors. We have demonstrated the effectiveness of DeepLNE by showing that for simple models such as the Müller-Brown potential and alanine dipeptide, the progression along the path variable closely approximates the ideal reaction coordinate. However, DeepLNE is computationally expensive for realistic systems needing many descriptors and limited in its ability to handle multi-state reactions. Here, we present DeepLNE++, which uses a knowledge distillation approach to significantly accelerate the evaluation of DeepLNE, making it feasible to compute free energy landscapes for large and complex biomolecular systems. In addition, DeepLNE++ encodes system-specific knowledge within a supervised multitasking framework, enhancing its versatility and effectiveness.

A novel statistical mechanical methodology is applied to clusters of N ≤ 7 atoms. Exact statistical analogs for any energy derivative of entropy ∂mS/∂Em are used in rigorous microcanonical Monte Carlo simulations to vastly enlarge the pool of measurable thermodynamic properties relative to previous work. All analogs are given for two alternative partition functions of the microcanonical ensemble. Coarse grained energy distributions are used to establish the existence of melting transitions. LJ7, LJ5, and LJ4 are found to exhibit trimodal distributions, a feature not being reported before. Varieties of combinations of entropy derivatives are tested for a direct detection of the melting region. It is shown that for such a purpose, derivatives of at least fourth order are necessary.

Lipid rafts are nanoscopic assemblies of sphingolipids, cholesterol, and specific membrane proteins. They are believed to underlie the experimentally observed lateral heterogeneity of eukaryotic plasma membranes and implicated in many cellular processes, such as signaling and trafficking. Ternary model membranes consisting of saturated lipids, unsaturated lipids, and cholesterol are common proxies because they exhibit phase coexistence between a liquid-ordered (lo) and liquid-disordered (ld) phase and an associated critical point. However, plasma membranes are also asymmetric in terms of lipid type, lipid abundance, leaflet tension, and corresponding cholesterol distribution, suggesting that rafts cannot be examined separately from questions about elasticity, curvature torques, and internal mechanical stresses. Unfortunately, it is challenging to capture this wide range of physical phenomenology in a single model that can access sufficiently long length- and time scales. Here we extend the highly coarse-grained Cooke model for lipids, which has been extensively characterized on the curvature-elastic front, to also represent raft-like lo/ld mixing thermodynamics. In particular, we capture the shape and tie lines of a coexistence region that narrows upon cholesterol addition, terminates at a critical point, and has coexisting phases that reflect key differences in membrane order and lipid packing. We furthermore examine elasticity and lipid diffusion for both phase separated and pure systems and how they change upon the addition of cholesterol. We anticipate that this model will enable significant insight into lo/ld phase separation and the associated question of lipid rafts for membranes that have compositionally distinct leaflets that are likely under differential stress-like the plasma membrane.

In spintronics, a perennial goal has been the generation of organic spin-bearing semiconductor materials with magnetic ordering stable at room temperature. To this end, the class of transition metal phthalocyanines has shown much promise in fulfilling this ambition. In particular, alpha-phase cobalt (II) phthalocyanine (α-CoPc) exhibits strong antiferromagnetic exchange interactions producing a long range order up to ∼100 K. However, the underlying mechanism by which this magnetic interaction proceeds is not well understood. In this report, a simple mechanism has been proposed based on the Hubbard Hamiltonian, which elucidates the exchange coupling in α-CoPc. The mechanism provides stipulations for increasing the magnetic coupling, and this directs to a proposal that substitution of the central cobalt ion for rhodium will lead to a significant increase in coupling strength. The strength of this exchange interaction has been evaluated using broken symmetry hybrid exchange density functional theory and indicates that the novel rhodium (II) phthalocyanine system is indeed predicted to exhibit significantly stronger magnetic ordering. This study, therefore, identifies the coupling mechanism in α-CoPc as primarily attributable to kinetic exchange, explains its previously reported strong coupling relative to its first-row transition metal counterparts, and suggests that rhodium (II) phthalocyanine is likely to exhibit stable magnetic ordering at room temperature.

We have implemented a cluster-continuum method using density functional theory to model water clusters and various charged species derived from water. The two aims of this study are to determine the minimal basis required for proper modeling of water autoionization and to determine the minimum number of explicit water molecules required to properly model the energetics of solvation. The thermodynamics of water autoionization converge following a modified power law to deliver chemically accurate values of the Gibbs energy change for autoionization with tractably small clusters. Convergence is slower and not exponential as assumed in previous work. We identify the n = 21 set of clusters as the first effectively bulk water like clusters that can capture the energetic influence of both the first and second solvation shells. In this cluster, a water molecule is encapsulated in the center of a closed shell of other water molecules that hydrogen bond to form five-membered rings. The total energy change for clusters with n ≥ 21 calculated using both the RPBE-D3 and ωB97X-D exchange-correlation functionals with the 6-311+G** basis set is shown to deliver good approximations to the free energy change at 298 K. This is true even though neither functional models the individual enthalpy or entropy contributions particularly well.

We compare two recently developed strategies, implemented in open source software packages, for computing linear optical spectra in condensed phase environments in the presence of nonadiabatic effects. Both approaches rely on computing excitation energy and transition dipole fluctuations along molecular dynamics (MD) trajectories, treating molecular and environmental degrees of freedom on the same footing. Spectra are then generated in two ways: in the recently developed Gaussian non-Condon theory, the linear response functions are computed in terms of independent adiabatic excited states, with non-Condon effects described through spectral densities of transition dipole fluctuations. For strongly coupled excited states, we instead parameterize a linear vibronic coupling Hamiltonian directly from spectral densities of energy fluctuations and diabatic couplings computed along the MD trajectory. The optical spectrum is then calculated using powerful, numerically exact tensor-network approaches. Both the electronic structure calculations to sample system fluctuations and the quantum dynamics simulations using tensor-network methods are carried out on graphics processing units, enabling rapid calculations on complex condensed phase systems. We assess the performance of the approaches using model systems in the presence of a conical intersection and the pyrazine molecule in different solvent environments.

Heterojunction catalysts in the field of hydrogen evolution reaction (HER) from electrocatalytic water splitting have recently become a hot research topic. In this paper, we systematically calculated the HER catalytic performance of a MoS2/CoS2 heterojunction for the first time, considering the effect of edge sites regulation, strain and electric field. The results indicate that the MoS2/CoS2 heterojunction exhibits synergistic catalytic performance compared to MoS2 and CoS2, the HER catalytic activity of which can be improved by exposing more edge sites or regulating the S content on the edges, with an optimized ratio of 25%. Surprisingly, applying strain has a slight effect on the catalytic activity of the edge, however, an obvious effect on the basal plane. For example, applying 2% tensile strain on the MoS2/CoS2 heterojunction can improve the edge catalytic performance by 13%, and for the basal plane, this value can reach 92%. In this case, the catalytic performance of the basal plane is better than that of the edge with 2% and without strain. Since the basal plane accounts for the majority of the two-dimensional catalysts, the catalytic performance of the basal plane is generally much lower than that of the edge. This discovery is of great significance, which means by adjusting strain, the catalytic performance of the heterojunction catalyst is likely to be improved by orders of magnitude. Moreover, considering the actual experimental process, we also calculated the effect of the electric field and found that 0.7 V/Å electric field can enhance the HER catalytic activity of the MoS2/CoS2 heterojunction by 23%.

Chiral structures, breaking spatial inversion symmetry, exhibit non-zero chiroptical activity (COA) due to the coupling between their electric and magnetic responses under external electromagnetic fields, an effect absent in achiral systems. Non-magnetic chiral structures also exhibit Chiral-Induced Spin Selectivity (CISS), primarily detected in two terminal measurements in the linear regime, where spin selection emerges without external magnetic influence. Despite the different origins of these physical phenomena, our model captures the relevant physics required to address CISS as an intrinsic molecular effect with the basic ingredients: (i) chirality/inversion asymmetry, (ii) meV atomic spin-orbit coupling, and (iii) decoherence as a source of reciprocity breaking. In this work, we derived how the electronic system couples with polarized electromagnetic radiation to yield a spin-dependent polarization rotation power, quantified through the Rosenfeld tensor, predicting characteristic spin signatures in the COA. The model also predicts that a net spin polarization manifests in the molecular terminations that have been surmised as an explanation for chiral species separation of racemic mixtures and interactions with surface magnetic domains. A recent sensitive spectroscopic measurement of electron transfer in donor-acceptor complexes is consistent with the standalone CISS effect.

By simulation and asymptotic theory, we investigate the transition-path time of a one-dimensional finite-mass reaction coordinate crossing a double-well potential in the presence of non-Markovian friction. First, we consider single-exponential memory kernels and demonstrate that memory accelerates transition paths compared to the Markovian case, especially in the low-mass/high-friction limit. Then, we generalize to multi-exponential kernels and construct an asymptotic formula for the transition-path time that compares well with simulation data.