Introduction. In: Chipot C, Pohorille A, editors. Free Energy Calculations. Berlin: Springer Berlin Heidelberg; 2007. pp. 1-31. [DOI: 10.1007/978-3-540-38448-9_1]

Self-assembly of nanocrystals into strongly electronically coupled all-inorganic supercrystals

Colloidal nanocrystals of metals, semiconductors, and other functional materials can self-assemble into long-range ordered crystalline and quasicrystalline phases, but insulating organic surface ligands prevent the development of collective electronic states in ordered nanocrystal assemblies. We reversibly self-assembled colloidal nanocrystals of gold, platinum, nickel, lead sulfide, and lead selenide with conductive inorganic ligands into supercrystals exhibiting optical and electronic properties consistent with strong electronic coupling between the constituent nanocrystals. The phase behavior of charge-stabilized nanocrystals can be rationalized and navigated with phase diagrams computed for particles interacting through short-range attractive potentials. By finely tuning interparticle interactions, the assembly was directed either through one-step nucleation or nonclassical two-step nucleation pathways. In the latter case, the nucleation was preceded by the formation of two metastable colloidal fluids.

Carbon dioxide capture in 2,2'-iminodiethanol aqueous solution from ab initio molecular dynamics simulations

The reaction of carbon dioxide (CO₂) with aqueous 2,2'-iminodiethanol (trivial name is diethanolamine: DEA) has been investigated using both blue moon ensemble and metadynamics approaches combined with ab initio molecular dynamics (AIMD) simulations. A spontaneous direct proton transfer from DEA zwitterion (DEAZW) to DEA but not to H₂O has been observed in straightforward AIMD simulation in the time scale of ps. The ab initio free-energy calculations reproduced the overall free-energy difference, predicting the ionic products DEA carbamate ion (DEAC) and the protonated DEA (DEAH). The computed free-energy barrier for the first reaction step, which is the CO₂ binding (48 kJ mol^-1), is found to agree reasonably well with the available experimental data (52-56 kJ mol^-1). By contrast, the barriers for the next step, the deprotonation of zwitterion realized either via reaction with DEA or H₂O, are underestimated by 25-35 kJ mol^-1 compared to the experimental reference. A part of this error is attributed to the neglected reversible work needed to bring two reactants together, which might significantly contribute to the free-energy of activation of bimolecular reactions in a dilute solution. The computed free-energy profile is compared with our results [Y. Kubota et al., J. Chem. Phys. 146, 094303 (2017)] for the same reaction in 2-aminoethanol (trivial name is monoethanolamine: MEA).

Combine umbrella sampling with integrated tempering method for efficient and accurate calculation of free energy changes of complex energy surface

Umbrella sampling is an efficient method for the calculation of free energy changes of a system along well-defined reaction coordinates. However, when there exist multiple parallel channels along the reaction coordinate or hidden barriers in directions perpendicular to the reaction coordinate, it is difficult for conventional umbrella sampling to reach convergent sampling within limited simulation time. Here, we propose an approach to combine umbrella sampling with the integrated tempering sampling method. The umbrella sampling method is applied to chemically more relevant degrees of freedom that possess significant barriers. The integrated tempering sampling method is used to facilitate the sampling of other degrees of freedom which may possess statistically non-negligible barriers. The combined method is applied to two model systems, butane and ACE-NME molecules, and shows significantly improved sampling efficiencies as compared to standalone conventional umbrella sampling or integrated tempering sampling approaches. Further analyses suggest that the enhanced performance of the new method come from the complemented advantages of umbrella sampling with a well-defined reaction coordinate and integrated tempering sampling in orthogonal space. Therefore, the combined approach could be useful in the simulation of biomolecular processes, which often involves sampling of complex rugged energy landscapes.