1
|
Wang R, Tiwary P. Atomic scale insights into NaCl nucleation in nanoconfined environments. Chem Sci 2024:d4sc04042b. [PMID: 39234215 PMCID: PMC11367593 DOI: 10.1039/d4sc04042b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2024] [Accepted: 08/23/2024] [Indexed: 09/06/2024] Open
Abstract
In this work we examine the nucleation from NaCl aqueous solutions within nano-confined environments, employing enhanced sampling molecular dynamics simulations integrated with machine learning-derived reaction coordinates. Through our simulations, we successfully induce phase transitions between solid, liquid, and a hydrated phase, typically observed at lower temperatures in bulk environments. Interestingly, while generally speaking nano-confinement serves to stabilize the solid phase and elevate melting points, there are subtle variations in the thermodynamics of competing phases with the precise extent of confinement. Our simulations explain these findings by underscoring the significant role of water, alongside ion aggregation and subtle, anisotropic dielectric behavior, in driving nucleation within nano-confined environments. This report thus provides a framework for sampling, analyzing and understanding nucleation processes under nano-confinement.
Collapse
Affiliation(s)
- Ruiyu Wang
- Institute for Physical Science and Technology, University of Maryland College Park MD 20742 USA
| | - Pratyush Tiwary
- Institute for Physical Science and Technology, University of Maryland College Park MD 20742 USA
- Department of Chemistry and Biochemistry, University of Maryland College Park MD 20742 USA
- University of Maryland Institute for Health Computing Bethesda Maryland 20852 USA
| |
Collapse
|
2
|
Al-Awad AS, Batet L, Rives R, Sedano L. Stochastic computer experiments of the thermodynamic irreversibility of bulk nanobubbles in supersaturated and weak gas-liquid solutions. J Chem Phys 2024; 161:024503. [PMID: 38984961 DOI: 10.1063/5.0204665] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2024] [Accepted: 06/24/2024] [Indexed: 07/11/2024] Open
Abstract
Spontaneous gas-bubble nucleation in weak gas-liquid solutions has been a challenging topic in theory, experimentation, and computer simulations. In analogy with recent advances in crystallization and droplet formation studies, the diffusive-shielding stabilization and thermodynamic irreversibility of bulk nanobubble (bNB) mechanisms are revisited and deployed to characterize nucleation processes in a stochastic framework of computer experiments using the large-scale atomic/molecular massively parallel simulator code. Theoretical bases, assumptions, and limitations underlying the irreversibility hypothesis of bNBs, and their computational counterparts, are extensively described and illustrated. In essence, it is established that the irreversibility hypothesis can be numerically investigated by converging the system volume (due to the finiteness of interatomic forces) and the initial dissolved-gas concentration in the solution (due to the single-bNB limitation). Helium nucleation in liquid Pb17Li alloy is selected as a representative case study, where it exhibits typical characteristics of noble-gas/liquid-metal systems. The proposed framework lays down the bases on which the stability of gas-bNBs in weak and supersaturated gas-liquid solutions can be inferred and explained from a novel perspective. In essence, it stochastically marches toward a unique irreversible state along out-of-equilibrium nucleation/growth trajectories. Moreover, it does not attempt to characterize the interface or any interface-related properties, neither theoretically nor computationally. It was concluded that bNBs of a few tens of He-atoms are irreversible when dissolved-He concentrations in the weak gas-liquid solution are at least ∼50 and ∼105 mol m-3 at 600 and 1000 K (and ∼80 MPa), respectively, whereas classical molecular dynamics -estimated solubilities are at least two orders of magnitude smaller.
Collapse
Affiliation(s)
- Abdulrahman S Al-Awad
- Department of Physics, Universitat Politècnica de Catalunya-BarcelonaTech (UPC), Barcelona 08028, Spain
| | - Lluis Batet
- Department of Physics, Universitat Politècnica de Catalunya-BarcelonaTech (UPC), Barcelona 08028, Spain
| | - Ronny Rives
- Department of Physics, Universitat Politècnica de Catalunya-BarcelonaTech (UPC), Barcelona 08028, Spain
| | - Luis Sedano
- Department of Physics, Universitat Politècnica de Catalunya-BarcelonaTech (UPC), Barcelona 08028, Spain
- Instituto de Ciencia de Materiales de Barcelona (ICMAB/CSIC), Bellaterra 08193, Spain
| |
Collapse
|
3
|
Zou Z, Tiwary P. Enhanced Sampling of Crystal Nucleation with Graph Representation Learnt Variables. J Phys Chem B 2024. [PMID: 38502931 DOI: 10.1021/acs.jpcb.4c00080] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/21/2024]
Abstract
In this study, we present a graph neural network (GNN)-based learning approach using an autoencoder setup to derive low-dimensional variables from features observed in experimental crystal structures. These variables are then biased in enhanced sampling to observe state-to-state transitions and reliable thermodynamic weights. In our approach, we used simple convolution and pooling methods. To verify the effectiveness of our protocol, we examined the nucleation of various allotropes and polymorphs of iron and glycine in their molten states. Our graph latent variables, when biased in well-tempered metadynamics, consistently show transitions between states and achieve accurate thermodynamic rankings in agreement with experiments, both of which are indicators of dependable sampling. This underscores the strength and promise of our GNN variables for improved sampling. The protocol shown here should be applicable for other systems and other sampling methods.
Collapse
Affiliation(s)
- Ziyue Zou
- Department of Chemistry and Biochemistry, University of Maryland, College Park 20742, Maryland, United States
| | - Pratyush Tiwary
- Department of Chemistry and Biochemistry, University of Maryland, College Park 20742, Maryland, United States
- Institute for Physical Science and Technology, University of Maryland, College Park 20742, Maryland, United States
- University of Maryland Institute for Health Computing, Rockville, Maryland 20852, United States
| |
Collapse
|
4
|
Wang R, Mehdi S, Zou Z, Tiwary P. Is the Local Ion Density Sufficient to Drive NaCl Nucleation from the Melt and Aqueous Solution? J Phys Chem B 2024; 128:1012-1021. [PMID: 38262436 DOI: 10.1021/acs.jpcb.3c06735] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2024]
Abstract
Even though nucleation is ubiquitous in different science and engineering problems, investigating nucleation is extremely difficult due to the complicated ranges of time and length scales involved. In this work, we simulate NaCl nucleation in both molten and aqueous environments using enhanced sampling of all-atom molecular dynamics with deep-learning-based estimation of reaction coordinates. By incorporating various structural order parameters and learning the reaction coordinate as a function thereof, we achieve significantly improved sampling relative to traditional ad hoc descriptions of what drives nucleation, particularly in an aqueous medium. Our results reveal a one-step nucleation mechanism in both environments, with reaction coordinate analysis highlighting the importance of local ion density in distinguishing solid and liquid states. However, although fluctuations in the local ion density are necessary to drive nucleation, they are not sufficient. Our analysis shows that near the transition states, descriptors such as enthalpy and local structure become crucial. Our protocol proposed here enables robust nucleation analysis and phase sampling and could offer insights into nucleation mechanisms for generic small molecules in different environments.
Collapse
Affiliation(s)
- Ruiyu Wang
- Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, United States
| | - Shams Mehdi
- Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, United States
- Biophysics Program, University of Maryland, College Park, Maryland 20742, United States
| | - Ziyue Zou
- Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742, United States
| | - Pratyush Tiwary
- Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, United States
- Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742, United States
| |
Collapse
|
5
|
Ray D, Parrinello M. Kinetics from Metadynamics: Principles, Applications, and Outlook. J Chem Theory Comput 2023; 19:5649-5670. [PMID: 37585703 DOI: 10.1021/acs.jctc.3c00660] [Citation(s) in RCA: 21] [Impact Index Per Article: 21.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/18/2023]
Abstract
Metadynamics is a popular enhanced sampling algorithm for computing the free energy landscape of rare events by using molecular dynamics simulation. Ten years ago, Tiwary and Parrinello introduced the infrequent metadynamics approach for calculating the kinetics of transitions across free energy barriers. Since then, metadynamics-based methods for obtaining rate constants have attracted significant attention in computational molecular science. Such methods have been applied to study a wide range of problems, including protein-ligand binding, protein folding, conformational transitions, chemical reactions, catalysis, and nucleation. Here, we review the principles of elucidating kinetics from metadynamics-like approaches, subsequent methodological developments in this area, and successful applications on chemical, biological, and material systems. We also highlight the challenges of reconstructing accurate kinetics from enhanced sampling simulations and the scope of future developments.
Collapse
Affiliation(s)
- Dhiman Ray
- Atomistic Simulations, Italian Institute of Technology, Via Enrico Melen 83, 16152 Genova, Italy
| | - Michele Parrinello
- Atomistic Simulations, Italian Institute of Technology, Via Enrico Melen 83, 16152 Genova, Italy
| |
Collapse
|
6
|
Bal KM, Neyts EC. Extending and validating bubble nucleation rate predictions in a Lennard-Jones fluid with enhanced sampling methods and transition state theory. J Chem Phys 2022; 157:184113. [DOI: 10.1063/5.0120136] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
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.
Collapse
Affiliation(s)
- Kristof M. Bal
- Department of Chemistry and NANOlab Center of Excellence, University of Antwerp, Universiteitsplein 1, 2610 Antwerp, Belgium
| | - Erik C. Neyts
- Department of Chemistry and NANOlab Center of Excellence, University of Antwerp, Universiteitsplein 1, 2610 Antwerp, Belgium
| |
Collapse
|
7
|
Exploring the Energy Landscape of Riboswitches Using Collective Variables Based on Tertiary Contacts. J Mol Biol 2022; 434:167788. [PMID: 35963460 PMCID: PMC10042644 DOI: 10.1016/j.jmb.2022.167788] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2022] [Revised: 08/03/2022] [Accepted: 08/07/2022] [Indexed: 12/24/2022]
Abstract
Messenger RNA regulatory elements, such as riboswitches, can display a high degree of flexibility. By characterizing their energy landscapes, and corresponding distributions of 3D configurations, structure-function relationships can be elucidated. Molecular dynamics simulation with enhanced sampling is an important strategy used to computationally access free energy landscapes characterizing the accessible 3D conformations of RNAs. While tertiary contacts are thought to play important roles in RNA dynamics, it is difficult, in explicit solvent, to sample the formation and breakage of tertiary contacts, such as helix-helix interactions, pseudoknot interactions, and junction interactions, while maintaining intact secondary structure elements. To this end, we extend previously developed collective variables and metadynamics efforts, to establish a simple metadynamics protocol, which utilizes only one collective variable, based on multiple tertiary contacts, to characterize the underlying free energy landscape of any RNA molecule. We develop a modified collective variable, the tertiary contacts distance (QTC), which can probe the formation and breakage of all or selectively chosen tertiary contacts of the RNA. The SAM-I riboswitch in the presence of three ionic and substrate conditions was investigated and validated against the structure ensemble previously generated using SAXS experiments. This efficient and easy to implement all-atom MD simulation based approach incorporating metadynamics to study RNA conformational dynamics can also be transferred to any other type of biomolecule.
Collapse
|
8
|
Beyerle ER, Mehdi S, Tiwary P. Quantifying Energetic and Entropic Pathways in Molecular Systems. J Phys Chem B 2022; 126:3950-3960. [PMID: 35605180 DOI: 10.1021/acs.jpcb.2c01782] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
When examining dynamics occurring at nonzero temperatures, both energy and entropy must be taken into account to describe activated barrier crossing events. Furthermore, good reaction coordinates need to be constructed to describe different metastable states and the transition mechanisms between them. Here we use a physics-based machine learning method called state predictive information bottleneck (SPIB) to find nonlinear reaction coordinates for three systems of varying complexity. SPIB is able to correctly predict an entropic bottleneck for an analytical flat-energy double-well system and identify the entropy- and energy-dominated pathways for an analytical four-well system. Finally, for a simulation of benzoic acid permeation through a lipid bilayer, SPIB is able to discover the entropic and energetic barriers to the permeation process. Given these results, we thus establish that SPIB is a reasonable and robust method for finding the important entropy, energy, and enthalpy barriers in physical systems, which can then be used to enhance the understanding and sampling of different activated mechanisms.
Collapse
Affiliation(s)
- Eric R Beyerle
- Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20740, United States
| | - Shams Mehdi
- Biophysics Program and Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, United States
| | - Pratyush Tiwary
- Department of Chemistry and Biochemistry and Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, United States
| |
Collapse
|
9
|
Ketkaew R, Creazzo F, Luber S. Machine Learning-Assisted Discovery of Hidden States in Expanded Free Energy Space. J Phys Chem Lett 2022; 13:1797-1805. [PMID: 35171614 DOI: 10.1021/acs.jpclett.1c04004] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Collective variables (CVs) are crucial parameters in enhanced sampling calculations and strongly impact the quality of the obtained free energy surface. However, many existing CVs are unique to and dependent on the system they are constructed with, making the developed CV non-transferable to other systems. Herein, we develop a non-instructor-led deep autoencoder neural network (DAENN) for discovering general-purpose CVs. The DAENN is used to train a model by learning molecular representations upon unbiased trajectories that contain only the reactant conformers. The prior knowledge of nonconstraint reactants coupled with the here-introduced topology variable and loss-like penalty function are only required to make the biasing method able to expand its configurational (phase) space to unexplored energy basins. Our developed autoencoder is efficient and relatively inexpensive to use in terms of a priori knowledge, enabling one to automatically search for hidden CVs of the reaction of interest.
Collapse
Affiliation(s)
- Rangsiman Ketkaew
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Fabrizio Creazzo
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Sandra Luber
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| |
Collapse
|
10
|
Zou Z, Tsai ST, Tiwary P. Toward Automated Sampling of Polymorph Nucleation and Free Energies with the SGOOP and Metadynamics. J Phys Chem B 2021; 125:13049-13056. [PMID: 34788047 DOI: 10.1021/acs.jpcb.1c07595] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Understanding the driving forces behind the nucleation of different polymorphs is of great importance for material sciences and the pharmaceutical industry. This includes understanding the reaction coordinate that governs the nucleation process and correctly calculating the relative free energies of different polymorphs. Here, we demonstrate, for the prototypical case of urea nucleation from the melt, how one can learn such a one-dimensional reaction coordinate as a function of prespecified order parameters and use it to perform efficient biased all-atom molecular dynamics simulations. The reaction coordinate is learnt as a function of the generic thermodynamic and structural order parameters using the "spectral gap optimization of order parameters (SGOOP)" approach [Tiwary, P. and Berne, B. J. Proc. Natl. Acad. Sci. U.S.A. (2016)] and is biased using well-tempered metadynamics simulations. The reaction coordinate gives insights into the role played by different structural and thermodynamics order parameters, and the biased simulations obtain accurate relative free energies for different polymorphs. This includes an accurate prediction of the approximate pressure at which urea undergoes a phase transition and one of the metastable polymorphs becomes the most stable conformation. We believe the ideas demonstrated in this work will facilitate efficient sampling of nucleation in complex, generic systems.
Collapse
Affiliation(s)
- Ziyue Zou
- Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742, United States
| | - Sun-Ting Tsai
- Department of Physics, University of Maryland, College Park, Maryland 20742, United States.,Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, United States
| | - Pratyush Tiwary
- Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742, United States.,Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, United States
| |
Collapse
|
11
|
Tsai ST, Smith Z, Tiwary P. SGOOP-d: Estimating Kinetic Distances and Reaction Coordinate Dimensionality for Rare Event Systems from Biased/Unbiased Simulations. J Chem Theory Comput 2021; 17:6757-6765. [PMID: 34662516 DOI: 10.1021/acs.jctc.1c00431] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Understanding kinetics including reaction pathways and associated transition rates is an important yet difficult problem in numerous chemical and biological systems, especially in situations with multiple competing pathways. When these high-dimensional systems are projected on low-dimensional coordinates, which are often needed for enhanced sampling or for interpretation of simulations and experiments, one can end up losing the kinetic connectivity of the underlying high-dimensional landscape. Thus, in the low-dimensional projection, metastable states might appear closer or further than they actually are. To deal with this issue, in this work, we develop a formalism that learns a multidimensional yet minimally complex reaction coordinate (RC) for generic high-dimensional systems. When projected along this RC, all possible kinetically relevant pathways can be demarcated and the true high-dimensional connectivity is maintained. One of the defining attributes of our method lies in that it can work on long unbiased simulations as well as biased simulations often needed for rare event systems. We demonstrate the utility of the method by studying a range of model systems including conformational transitions in a small peptide Ace-Ala3-Nme, where we show how two-dimensional and three-dimensional RCs found by our previously published spectral gap optimization method "SGOOP" [Tiwary, P. and Berne, B. J. Proc. Natl. Acad. Sci. 2016, 113, 2839] can capture the kinetics for 23 and all 28 out of the 28 dominant state-to-state transitions, respectively.
Collapse
Affiliation(s)
- Sun-Ting Tsai
- Department of Physics and Institute for Physical Science and Technology, University of Maryland, College Park 20742, Maryland, United States
| | - Zachary Smith
- Biophysics Program and Institute for Physical Science and Technology, University of Maryland, College Park 20742, Maryland, United States
| | - Pratyush Tiwary
- Department of Chemistry and Biochemistry and Institute for Physical Science and Technology, University of Maryland, College Park 20742, Maryland, United States
| |
Collapse
|
12
|
Bal KM. Nucleation rates from small scale atomistic simulations and transition state theory. J Chem Phys 2021; 155:144111. [PMID: 34654300 DOI: 10.1063/5.0063398] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
The evaluation of nucleation rates from molecular dynamics trajectories is hampered by the slow nucleation time scale and impact of finite size effects. Here, we show that accurate nucleation rates can be obtained in a very general fashion relying only on the free energy barrier, transition state theory, and a simple dynamical correction for diffusive recrossing. In this setup, the time scale problem is overcome by using enhanced sampling methods, in casu metadynamics, whereas the impact of finite size effects can be naturally circumvented by reconstructing the free energy surface from an appropriate ensemble. Approximations from classical nucleation theory are avoided. We demonstrate the accuracy of the approach by calculating macroscopic rates of droplet nucleation from argon vapor, spanning 16 orders of magnitude and in excellent agreement with literature results, all from simulations of very small (512 atom) systems.
Collapse
Affiliation(s)
- Kristof M Bal
- Department of Chemistry and NANOlab Center of Excellence, University of Antwerp, Universiteitsplein 1, 2610 Antwerp, Belgium
| |
Collapse
|
13
|
Ediger MD, Jensen L, Manolopoulos DE, Martinez TJ, Michaelides A, Reichman DR, Sherrill CD, Shi Q, Straub JE, Vega C, Wang LS, Brigham EC, Lian T. JCP Emerging Investigator Special Collection 2019. J Chem Phys 2020; 153:110402. [PMID: 32962387 DOI: 10.1063/5.0021946] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Affiliation(s)
- Mark D Ediger
- Department of Chemistry, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Lasse Jensen
- Department of Chemistry, Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | - David E Manolopoulos
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - Todd J Martinez
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Angelos Michaelides
- Thomas Young Centre and London Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom and Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - David R Reichman
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Qiang Shi
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China; University of Chinese Academy of Sciences, Beijing 100049, China; and Physical Science Laboratory, Huairou National Comprehensive Science Center, Beijing 101407, China
| | - John E Straub
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, USA
| | - Carlos Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Lai-Sheng Wang
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| | | | - Tianquan Lian
- Department of Chemistry, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|