Shear viscosity prediction of alcohols, hydrocarbons, halogenated, carbonyl, nitrogen-containing, and sulfur compounds using the variable force fields

Machine Learning-Guided Adaptive Parametrization for Coupling Terms in a Mixed United-Atom/Coarse-Grained Model for Diphenylalanine Self-Assembly in Aqueous Ionic Liquids

Precise regulation of the peptide self-assembly into ordered nanostructures with intriguing properties has attracted intense attention. However, predicting peptide assembly at atomic resolution is a challenge due to both the structural flexibility of peptides and the associated huge computational costs. A machine learning-guided adaptive parametrization method was proposed for developing a mixed atomic and coarse-grained (CG) model through a multiobjective optimization strategy. Our model incorporates the united-atom (UA) model for diphenylalanine (P) and the polarizable electrostatic-variable coarse-grained (VaCG) model for aqueous ionic liquid [BMIM]+[BF4]- solution. In this mixed model, the coupling van der Waals (vdW) interaction is addressed by introducing virtual sites (VS) in the UA model to interact with solvent CG beads. The coupling parameters, including the electrostatic parameter and vdW parameters, are automatically optimized through ML-guided adaptive parametrization. The performance of this model was tested by some microstructural properties, e.g., the average number of P-P intermolecular hydrogen bonds (HBs) and radius distribution functions (RDFs) between P and different fragments of IL, in comparison with all-atom (AA) simulations. The computational cost is significantly reduced using such a parametrization scheme, which could search tens of thousands of force-field parameter sets, while needing only a small fraction of them to be assessed with molecular dynamics (MD) simulations. We used such a mixed resolution model to investigate the self-assembly in IL-water mixtures with variants of IL concentration (X). The long-range-ordered fibril structure is formed in a pure water system (X = 0). With an increase of IL concentrations, the formation of an ordered self-assembly nanostructure is prohibited, instead forming branched fibril at X = 2 mol % or amorphous aggregates when X > 10 mol %, resulting from the interplay between π-stacking and HB interactions between P and IL. The qualitative agreement between the simulated structures and the observed morphologies in experiments indicates the applicability of ML-guided parametrization strategy in the study of complex systems, such as polymers, lipid bilayers, and polysaccharides.

Treating Polarization Effects in Charged and Polar Bio-Molecules Through Variable Electrostatic Parameters

Polarization plays important roles in charged and hydrogen bonding containing systems. Much effort ranging from the construction of physics-based models to quantum mechanism (QM)-based and machine learning (ML)-assisted models have been devoted to incorporating the polarization effect into the conventional force fields at different levels, such as atomic and coarse grained (CG). The application of polarizable force fields or polarization models was limited by two aspects, namely, computational cost and transferability. Different from physics-based models, no predetermining parameters were required in the QM-based approaches. Taking advantage of both the accuracy of QM calculations and efficiency of molecular mechanism (MM) and ML, polarization effects could be treated more efficiently while maintaining the QM accuracy. The computational cost could be reduced with variable electrostatic parameters, such as the charge, dipole, and electronic dielectric constant with the help of linear scaling fragmentation-based QM calculations and ML models. Polarization and entropy effects on the prediction of partition coefficient of druglike molecules are demonstrated by using both explicit or implicit all-atom molecular dynamics simulations and machine learning-assisted models. Directions and challenges for future development are also envisioned.

Molecular partition coefficient from machine learning with polarization and entropy embedded atom-centered symmetry functions

Efficient prediction of the partition coefficient (log P) between polar and non-polar phases could shorten the cycle of drug and materials design. In this work, a descriptor, named 〈q - ACSFs〉_conf, is proposed to take the explicit polarization effects in the polar phase and the conformation ensemble of energetic and entropic significance in the non-polar phase into consideration. The polarization effects are involved by embedding the partial charge directly derived from force fields or quantum chemistry calculations into the atom-centered symmetry functions (ACSFs), together with the entropy effects, which are averaged according to the Boltzmann distribution of different conformations taken from the similarity matrix. The model was trained with high-dimensional neural networks (HDNNs) on a public dataset PhysProp (with 41 039 samples). Satisfactory log P prediction performance was achieved on three other datasets, namely, Martel (707 molecules), Star & Non-Star (266) and Huuskonen (1870). The present 〈q - ACSFs〉_conf model was also applicable to n-carboxylic acids with the number of carbons ranging from 2 to 14 and 54 kinds of organic solvent. It is easy to apply the present method to arbitrary sized systems and give a transferable atom-based partition coefficient.

Computational and data driven molecular material design assisted by low scaling quantum mechanics calculations and machine learning

Electronic structure methods based on quantum mechanics (QM) are widely employed in the computational predictions of the molecular properties and optoelectronic properties of molecular materials. The computational costs of these QM methods, ranging from density functional theory (DFT) or time-dependent DFT (TDDFT) to wave-function theory (WFT), usually increase sharply with the system size, causing the curse of dimensionality and hindering the QM calculations for large sized systems such as long polymer oligomers and complex molecular aggregates. In such cases, in recent years low scaling QM methods and machine learning (ML) techniques have been adopted to reduce the computational costs and thus assist computational and data driven molecular material design. In this review, we illustrated low scaling ground-state and excited-state QM approaches and their applications to long oligomers, self-assembled supramolecular complexes, stimuli-responsive materials, mechanically interlocked molecules, and excited state processes in molecular aggregates. Variable electrostatic parameters were also introduced in the modified force fields with the polarization model. On the basis of QM computational or experimental datasets, several ML algorithms, including explainable models, deep learning, and on-line learning methods, have been employed to predict the molecular energies, forces, electronic structure properties, and optical or electrical properties of materials. It can be conceived that low scaling algorithms with periodic boundary conditions are expected to be further applicable to functional materials, perhaps in combination with machine learning to fast predict the lattice energy, crystal structures, and spectroscopic properties of periodic functional materials.