Wiley Interdisciplinary Reviews: Computational Molecular Science最新文献

Ferroelectrics with electrically switchable spontaneous polarizations can be used for information storage, where a low switching barrier is favorable to reduce the energy cost and enhance the speed for data writing. Meanwhile their robustness at working temperature should be ensured, which is a challenge for the designs of low-barrier ferroelectrics. Here we review several types of ferroelectric mechanisms that may render both low switching barriers and room-temperature robustness, which have been theoretically proposed in previous studies. (1) The prediction of sliding ferroelectricity with ultralow switching barriers has been experimentally confirmed in a series of van der Waals layers, which may enable convenient electrical control of various physical properties in 2D materials, like magnetic, photovoltaic, valleytronic and topological properties. (2) Hydrogen-bonded ferroelectricity spontaneously formed by head-to-tail chains can be switched by proton-transfer crossing a low barrier, and a mechanism of ultra-high piezoelectricity utilizing the specific features of hydrogen bonding has been proposed. (3) High-ionicity ferroelectricity induced by covalent-like ionic bondings may entail high polarizations and low barriers during switching, which is attributed to the features of long-range Coulomb interaction, and the long ion-displacements crossing unitcell may give rise to unconventional ferroelectricity with quantized polarizations even in crystals of non-ferroelectric point groups. Those low-barrier ferroelectric mechanisms may bring in both new physics and technological advances, which are to be further explored.

This article is categorized under:

Biological membranes (biomembranes) are one of the most complicated structures that allow life to exist. Investigating their structure, dynamics, and function is crucial for advancing our knowledge of cellular mechanisms and developing novel therapeutic strategies. However, experimental investigation of many biomembrane phenomena is challenging due to their compositional and structural complexity, as well as the inherently multi-scalar features. Computational approaches, particularly molecular dynamics (MD) simulations, have emerged as powerful tools for addressing the atomic details of biomembrane systems, driving breakthroughs in our understanding of biomembranes and their roles in cellular function. This review presents an overview of the latest advancements in related computational approaches, from force fields and model construction to MD simulations and trajectory analysis. We also discussed current hot research topics and challenges. Finally, we outline future directions, emphasizing the integration of force field development, enhanced sampling techniques, and data-driven approaches to accelerate the growth of this field in the years to come. We aim to equip readers with an understanding of the promise and limitations of emerging computational technologies in biomembrane systems and offer valuable recommendations for future research endeavors.

This article is categorized under:

Hybrid halide perovskite solar cells have been recognized as one of the most promising future photovoltaic technologies due to their demonstrated high-power conversion efficiency, versatile stoichiometry and low cost. However, degradation caused by environmental exposure and structural instability due to ionic defect migration hinders the commercialization of this technology. While the experimental studies try to understand the phenomenology of the degradation mechanisms and devise practical measures to improve the stability of these materials, theoretical studies have attempted to bridge the gaps in our understanding of the fundamental degradation mechanisms at different time and length scales. A deeper understanding of the physical and chemical phenomena at an atomic level through multiscale materials modeling is going to be crucial for the knowledge-based prognosis and design of future halide perovskites. There have been increased efforts in this direction in the last few years. However, the instability fundamentals explored through atomistic modeling and simulation methods have not been reviewed comprehensively in the literature yet. Therefore, this paper is an attempt to present a critical review, while identifying the existing gaps and opportunities in the investigation of the degradation and instability issues of the halide perovskites using computational methods. The review will primarily focus on the instability caused due to the intrinsic ionic defect migration and degradation due to thermal, moisture and light exposure. The findings from the simulation studies conducted primarily using density functional theory, ab initio molecular dynamics, classical molecular dynamics and machine learning methods will be presented.

This article is categorized under:

Recent efforts to synthetically expand drug-like chemical libraries have led to the emergence of unprecedently large virtual databases. This surge of make-on-demand molecular datasets has been received enthusiastically across the drug discovery community as a new paradigm. In several recent studies, virtual screening (VS) of larger make-on-demand collections resulted in the identification of novel molecules with higher potency and specificity compared to more conventional VS campaigns relying on smaller in-stock libraries. These results inspired ultra-large VS against various clinically relevant targets, including key proteins of the SARS-CoV-2 virus. As library sizes rapidly surpassed the billion compounds mark, new computational screening strategies emerged, shifting from conventional docking to fragment-based and machine learning-accelerated methods. These approaches significantly reduce computational demands of ultra-large screenings by lowering the number of molecules explicitly docked onto a target. Such strategies already demonstrated promise in evaluating libraries of tens of billions of molecules at relatively low computational cost. Herein, we review recent advancements in structure-based methods for ultra-large virtual screening that drug discovery practitioners have adopted to explore the ever-expanding chemical universe.

This article is categorized under:

We describe the development, implementation, and application of a polarizable QM/MM strategy, based on the AMOEBA polarizable force field, for calculating molecular properties and performing dynamics of molecular systems embedded in complex matrices. We show that polarizable QM/MM is a well-understood, mature technology that can be deployed using a state-of-the-art implementation that combines efficient numerical methods and linear scaling techniques. Thanks to these numerical advances and to the availability of parameters for a wide manifold of systems in the AMOEBA force field, polarizable QM/AMOEBA can be used for advanced production applications, that range from the prediction of spectroscopies to ground- and excited-state multiscale ab initio molecular dynamics simulations.

This article is categorized under:

We present MultiPsi, an open-source MCSCF program for the calculation of ground and excited states properties of strongly correlated systems. The program currently implements a general MCSCF code with excited states available using either state-averaging or linear response. It is written in a highly modular fashion using Python/C++ which makes it well suited as a development platform, enabling easy prototyping of novel methods, and as a teaching tool using interactive notebooks. The code is also very efficient and designed for modern high-performance computing environments using hybrid OpenMP/MPI parallelization. This efficiency is demonstrated with the calculation of the CASSCF energy and linear response of a molecule with more than 700 atoms as well as a fully optimized conventional CI calculation on more than 400 billion determinants.

This article is categorized under:

Drug discovery is a daunting and failure-prone task. A critical process in this research field is represented by the biological target and pocket identification steps as they heavily determine the subsequent efforts in selecting a putative ligand, most often a small molecule. Finding “ligandable” pockets, namely protein cavities that may accept a drug-like binder is instrumental to the more general and drug discovery oriented “druggability” estimation process. While high-throughput experimental techniques exist to identify putative binding sites other than the orthosteric one, these techniques are relatively expensive and not so commonly available in labs. In this regard, computational means of detecting ligandable pockets are advisable for their inexpensiveness and speed. These methods can become, in principle, particularly predictive when supported by machine learning methodologies that provide the modeling framework. As with any data-driven effort, the outcome critically depends on the input data, its featurization process and possible associated biases. Also, the machine learning task, (supervised/unsupervised) the learning method, and the possible usage of molecular dynamics data considerably shape the inherent assumptions of the modeling step. Defining a proper quantitative thermodynamic and/or kinetic score (or label) is key to the modeling process; here we revise literature and propose residence time as a novel ideal indicator of ligandability. Interestingly the vast majority of the methods does not keep into consideration kinetics nor thermodynamics when devising predictors.

This article is categorized under:

The dramatic progress of experimental attosecond science has called for the development of new theoretical and computational tools capable of accurately model the correlated electron dynamics triggered by attosecond molecular photoionization. We describe the nature and the main outcome of this development, with particular focus on the B-spline ADC and RCS-ADC ab initio methods.

This article is categorized under:

Aerosol particles are important for our global climate, but the mechanisms and especially the relative importance of various vapors for new particles formation (NPF) remain uncertain. Quantum chemical (QC) studies on organic enhanced nucleation has for the past couple of decades attracted immense attention, but very little remains known about the exact organic compounds that potentially are important for NPF. Here we comprehensively review the QC literature on atmospheric cluster formation involving organic compounds. We outline the potential cluster systems that should be further investigated. Cluster formation involving complex multi-functional organic accretion products warrant further investigations, but such systems are out of reach with currently applied methodologies. We suggest a “cluster of functional groups” approach to address this issue, which will allow for the identification of the potential structure of organic compounds that are involved in atmospheric NPF.

This article is categorized under: