Artificial Molecular Ratchets: Tools Enabling Endergonic Processes.

Non-equilibrium chemical systems underpin multiple domains of contemporary interest, including supramolecular chemistry, molecular machines, systems chemistry, prebiotic chemistry, and energy transduction. Experimental chemists are now pioneering the realization of artificial systems that can harvest energy away from equilibrium. In this tutorial Review, we provide an overview of artificial molecular ratchets: the chemical mechanisms enabling energy absorption from the environment. By focusing on the mechanism type-rather than the application domain or energy source-we offer a unifying picture of seemingly disparate phenomena, which we hope will foster progress in this fascinating domain of science.

A [3]-catenane non-autonomous molecular motor model: Geometric phase, no-pumping theorem, and energy transduction.

We study a model of a synthetic molecular motor-a [3]-catenane consisting of two small macrocycles mechanically interlocked with a bigger one-subjected to time-dependent driving using stochastic thermodynamics. The model presents nontrivial features due to the two interacting small macrocycles but is simple enough to be treated analytically in limiting regimes. Among the results obtained, we find a mapping into an equivalent [2]-catenane that reveals the implications of the no-pumping theorem stating that to generate net motion of the small macrocycles, both energies and barriers need to change. In the adiabatic limit (slow driving), we fully characterize the motor's dynamics and show that the net motion of the small macrocycles is expressed as a surface integral in parameter space, which corrects previous erroneous results. We also analyze the performance of the motor subjected to step-wise driving protocols in the absence and presence of an applied load. Optimization strategies for generating large currents and maximizing free energy transduction are proposed. This simple model provides interesting clues into the working principles of non-autonomous molecular motors and their optimization.

Kinetic Barrier Diagrams to Visualize and Engineer Molecular Nonequilibrium Systems.

Molecular nonequilibrium systems hold great promises for the nanotechnology of the future. Yet, their development is slowed by the absence of an informative representation. Indeed, while potential energy surfaces comprise in principle all the information, they hide the dynamic interplay of multiple reaction pathways underlying nonequilibrium systems, i.e., the degree of kinetic asymmetry. To offer an insightful visual representation of kinetic asymmetry, we extended an approach pertaining to catalytic networks, the energy span model, by focusing on system dynamics - rather than thermodynamics. Our approach encompasses both chemically and photochemically driven systems, ranging from unimolecular motors to simple self-assembly schemes. The obtained diagrams give immediate access to information needed to guide experiments, such as states' population, rate of machine operation, maximum work output, and effects of design changes. The proposed kinetic barrier diagrams offer a unifying graphical tool for disparate nonequilibrium phenomena.

Using Catalysis to Drive Chemistry Away from Equilibrium: Relating Kinetic Asymmetry, Power Strokes, and the Curtin-Hammett Principle in Brownian Ratchets.

Chemically fueled autonomous molecular machines are catalysis-driven systems governed by Brownian information ratchet mechanisms. One fundamental principle behind their operation is kinetic asymmetry, which quantifies the directionality of molecular motors. However, it is difficult for synthetic chemists to apply this concept to molecular design because kinetic asymmetry is usually introduced in abstract mathematical terms involving experimentally inaccessible parameters. Furthermore, two seemingly contradictory mechanisms have been proposed for chemically driven autonomous molecular machines: Brownian ratchet and power stroke mechanisms. This Perspective addresses both these issues, providing accessible and experimentally useful design principles for catalysis-driven molecular machinery. We relate kinetic asymmetry to the Curtin-Hammett principle using a synthetic rotary motor and a kinesin walker as illustrative examples. Our approach describes these molecular motors in terms of the Brownian ratchet mechanism but pinpoints both chemical gating and power strokes as tunable design elements that can affect kinetic asymmetry. We explain why this approach to kinetic asymmetry is consistent with previous ones and outline conditions where power strokes can be useful design elements. Finally, we discuss the role of information, a concept used with different meanings in the literature. We hope that this Perspective will be accessible to a broad range of chemists, clarifying the parameters that can be usefully controlled in the design and synthesis of molecular machines and related systems. It may also aid a more comprehensive and interdisciplinary understanding of biomolecular machinery.

Information thermodynamics for deterministic chemical reaction networks.

Information thermodynamics relates the rate of change of mutual information between two interacting subsystems to their thermodynamics when the joined system is described by a bipartite stochastic dynamics satisfying local detailed balance. Here, we expand the scope of information thermodynamics to deterministic bipartite chemical reaction networks, namely, composed of two coupled subnetworks sharing species but not reactions. We do so by introducing a meaningful notion of mutual information between different molecular features that we express in terms of deterministic concentrations. This allows us to formulate separate second laws for each subnetwork, which account for their energy and information exchanges, in complete analogy with stochastic systems. We then use our framework to investigate the working mechanisms of a model of chemically driven self-assembly and an experimental light-driven bimolecular motor. We show that both systems are constituted by two coupled subnetworks of chemical reactions. One subnetwork is maintained out of equilibrium by external reservoirs (chemostats or light sources) and powers the other via energy and information flows. In doing so, we clarify that the information flow is precisely the thermodynamic counterpart of an information ratchet mechanism only when no energy flow is involved.

Kinetic and energetic insights into the dissipative non-equilibrium operation of an autonomous light-powered supramolecular pump.

Natural and artificial autonomous molecular machines operate by constantly dissipating energy coming from an external source to maintain a non-equilibrium state. Quantitative thermodynamic characterization of these dissipative states is highly challenging as they exist only as long as energy is provided. Here we report on the detailed physicochemical characterization of the dissipative operation of a supramolecular pump. The pump transduces light energy into chemical energy by bringing self-assembly reactions to non-equilibrium steady states. The composition of the system under light irradiation was followed in real time by 1H NMR for four different irradiation intensities. The experimental composition and photon flow were then fed into a theoretical model describing the non-equilibrium dissipation and the energy storage at the steady state. We quantitatively probed the relationship between the light energy input and the deviation of the dissipative state from thermodynamic equilibrium in this artificial system. Our results provide a testing ground for newly developed theoretical models for photoactivated artificial molecular machines operating away from thermodynamic equilibrium.

Sample size dependence of tagged molecule dynamics in steady-state networks with bimolecular reactions: Cycle times of a light-driven pump.

Here, steady-state reaction networks are inspected from the viewpoint of individual tagged molecules jumping among their chemical states upon the occurrence of reactive events. Such an agent-based viewpoint is useful for selectively characterizing the behavior of functional molecules, especially in the presence of bimolecular processes. We present the tools for simulating the jump dynamics both in the macroscopic limit and in the small-volume sample where the numbers of reactive molecules are of the order of few units with an inherently stochastic kinetics. The focus is on how an ideal spatial "compartmentalization" may affect the dynamical features of the tagged molecule. Our general approach is applied to a synthetic light-driven supramolecular pump composed of ring-like and axle-like molecules that dynamically assemble and disassemble, originating an average ring-through-axle directed motion under constant irradiation. In such an example, the dynamical feature of interest is the completion time of direct/inverse cycles of tagged rings and axles. We find a surprisingly strong robustness of the average cycle times with respect to the system's size. This is explained in the presence of rate-determining unimolecular processes, which may, therefore, play a crucial role in stabilizing the behavior of small chemical systems against strong fluctuations in the number of molecules.

Insights from an information thermodynamics analysis of a synthetic molecular motor.

Information is physical, a realization that has transformed the physics of measurement and communication. However, the flow between information, energy and mechanics in chemical systems remains largely unexplored. Here we analyse a minimalist autonomous chemically driven molecular motor in terms of information thermodynamics, a framework that quantitatively relates information to other thermodynamic parameters. The treatment reveals how directional motion is generated by free energy transfer from chemical to mechanical (conformational and/or co-conformational) processes by 'energy flow' and 'information flow'. It provides a thermodynamic level of understanding of molecular motors that is general, complements previous analyses based on kinetics and has practical implications for machine design. In line with kinetic analysis, we find that power strokes do not affect the directionality of chemically driven machines. However, we find that power strokes can modulate motor velocity, the efficiency of free energy transfer and the number of fuel molecules consumed per cycle. This may help explain the role of such (co-)conformational changes in biomachines and illustrates the interplay between energy and information in chemical systems.

Nonequilibrium thermodynamics of light-induced reactions.

Current formulations of nonequilibrium thermodynamics of open chemical reaction networks only consider chemostats as free-energy sources sustaining nonequilibrium behaviors. Here, we extend the theory to include incoherent light as a source of free energy. We do so by relying on a local equilibrium assumption to derive the chemical potential of photons relative to the system they interact with. This allows us to identify the thermodynamic potential and the thermodynamic forces driving light-reacting chemical systems out-of-equilibrium. We use this framework to treat two paradigmatic photochemical mechanisms describing light-induced unimolecular reactions-namely, the adiabatic and diabatic mechanisms-and highlight the different thermodynamics they lead to. Furthermore, using a thermodynamic coarse-graining procedure, we express our findings in terms of commonly measured experimental quantities, such as quantum yields.

Nonequilibrium thermodynamics of non-ideal chemical reaction networks.

All current formulations of nonequilibrium thermodynamics of open chemical reaction networks rely on the assumption of non-interacting species. We develop a general theory that accounts for interactions between chemical species within a mean-field approach using activity coefficients. Thermodynamic consistency requires that rate equations do not obey standard mass-action kinetics but account for the interactions with concentration dependent kinetic constants. Many features of the ideal formulations are recovered. Crucially, the thermodynamic potential and the forces driving non-ideal chemical systems out of equilibrium are identified. Our theory is general and holds for any mean-field expression of the interactions leading to lower bounded free energies.

Thermodynamic efficiency in dissipative chemistry.

Chemical processes in closed systems inevitably relax to equilibrium. Living systems avoid this fate and give rise to a much richer diversity of phenomena by operating under nonequilibrium conditions. Recent experiments in dissipative self-assembly also demonstrated that by opening reaction vessels and steering certain concentrations, an ocean of opportunities for artificial synthesis and energy storage emerges. To navigate it, thermodynamic notions of energy, work and dissipation must be established for these open chemical systems. Here, we do so by building upon recent theoretical advances in nonequilibrium statistical physics. As a central outcome, we show how to quantify the efficiency of such chemical operations and lay the foundation for performance analysis of any dissipative chemical process.

Individual-Molecule Perspective Analysis of Chemical Reaction Networks: The Case of a Light-Driven Supramolecular Pump.

The first study in which stochastic simulations of a two-component molecular machine are performed in the mass-action regime is presented. This system is an autonomous molecular pump consisting of a photoactive axle that creates a directed flow of rings through it by exploiting light energy away from equilibrium. The investigation demonstrates that the pump can operate in two regimes, both experimentally accessible, in which light-driven steps can be rate-determining or not. The number of photons exploited by an individual molecular pump, as well as the precision of cycling and the overall efficiency, critically rely on the operating regime of the machine. This approach provides useful information not only to guide the chemical design of a self-assembling molecular device with desired features, but also to elucidate the effect of the environment on its performance, thus facilitating its experimental investigation.

Development and Implementation of Advanced Fitting Methods for the Calculation of Accurate Molecular Structures.

The determination of accurate equilibrium molecular structures plays a fundamental role for understanding many physical-chemical properties of molecules, ranging from the precise evaluation of the electronic structure to the analysis of dynamical and environmental effects in tuning their overall behavior. For this purpose the so-called semiexperimental approach, based on a nonlinear least-squares fit of the moments of inertia associated with a set of available isotopologues, allows one to obtain very accurate results, without the unfavorable computational cost characterizing high-level quantum chemical methods. In the present work the MSR (Molecular Structure Refinement) software for the determination of equilibrium structures by means of the semiexperimental approach is presented, and its implementation is discussed in some detail. The software, which is interfaced with a powerful graphical user interface, includes different optimization algorithms, an extended error analysis, and a number of advanced features, the most remarkable ones concerning the choice of internal coordinates and the method of predicate observations. In particular, a new black-box scheme for defining automatically a suitable set of nonredundant internal coordinates of A1 symmetry in place of the customary Z-matrix has been designed and tested. Finally, the implementation of the method of the predicate observations is discussed and validated for a set of test molecules. As an original application, the method is employed for the determination of the semiexperimental structure for the most stable conformer of glycine.

Structural features of the carbon-sulfur chemical bond: a semi-experimental perspective.

In this work semi-experimental and theoretical equilibrium geometries of 10 sulfur-containing organic molecules, as well as 4 oxygenated ones, are determined by means of a computational protocol based on density functional theory. The results collected in the present paper further enhance our online database of accurate semi-experimental equilibrium molecular geometries, adding 13 new molecules containing up to 8 atoms, for 12 of which the first semi-experimental equilibrium structure is reported, to the best of our knowledge. We focus in particular on sulfur-containing compounds, aiming both to provide new accurate data on some rather important chemical moieties, only marginally represented in the literature of the field, and to examine the structural features of carbon-sulfur bonds in the light of the previously presented linear regression approach. The structural changes issuing from substitution of oxygen by sulfur are discussed to get deeper insights on how modifications in electronic structure and nuclear potential can affect equilibrium geometries. With respect to our previous works, we perform non-linear constrained optimizations of equilibrium SE structures with a new general and user-friendly software under development in our group with updated definition of useful statistical indicators.

CC/DFT Route toward Accurate Structures and Spectroscopic Features for Observed and Elusive Conformers of Flexible Molecules: Pyruvic Acid as a Case Study.

The structures and relative stabilities as well as the rotational and vibrational spectra of the three low-energy conformers of pyruvic acid (PA) have been characterized using a state-of-the-art quantum-mechanical approach designed for flexible molecules. By making use of the available experimental rotational constants for several isotopologues of the most stable PA conformer, Tc-PA, the semiexperimental equilibrium structure has been derived. The latter provides a reference for the pure theoretical determination of the equilibrium geometries for all conformers, thus confirming for these structures an accuracy of 0.001 Å and 0.1 deg for bond lengths and angles, respectively. Highly accurate relative energies of all conformers (Tc-, Tt-, and Ct-PA) and of the transition states connecting them are provided along with the thermodynamic properties at low and high temperatures, thus leading to conformational enthalpies accurate to 1 kJ mol(-1). Concerning microwave spectroscopy, rotational constants accurate to about 20 MHz are provided for the Tt- and Ct-PA conformers, together with the computed centrifugal-distortion constants and dipole moments required to simulate their rotational spectra. For Ct-PA, vibrational frequencies in the mid-infrared region accurate to 10 cm(-1) are reported along with theoretical estimates for the transitions in the near-infrared range, and the corresponding infrared spectrum including fundamental transitions, overtones, and combination bands has been simulated. In addition to the new data described above, theoretical results for the Tc- and Tt-PA conformers are compared with all available experimental data to further confirm the accuracy of the hybrid coupled-cluster/density functional theory (CC/DFT) protocol applied in the present study. Finally, we discuss in detail the accuracy of computational models fully based on double-hybrid DFT functionals (mainly at the B2PLYP/aug-cc-pVTZ level) that avoid the use of very expensive CC calculations.

Semiexperimental equilibrium structures for building blocks of organic and biological molecules: the B2PLYP route.

The B2PLYP double hybrid functional, coupled with the correlation-consistent triple-ζ cc-pVTZ (VTZ) basis set, has been validated in the framework of the semiexperimental (SE) approach for deriving accurate equilibrium structures of molecules containing up to 15 atoms. A systematic comparison between new B2PLYP/VTZ results and several equilibrium SE structures previously determined at other levels, in particular B3LYP/SNSD and CCSD(T) with various basis sets, has put in evidence the accuracy and the remarkable stability of such model chemistry for both equilibrium structures and vibrational corrections. New SE equilibrium structures for phenylacetylene, pyruvic acid, peroxyformic acid, and phenyl radical are discussed and compared with literature data. Particular attention has been devoted to the discussion of systems for which lack of sufficient experimental data prevents a complete SE determination. In order to obtain an accurate equilibrium SE structure for these situations, the so-called templating molecule approach is discussed and generalized with respect to our previous work. Important applications are those involving biological building blocks, like uracil and thiouracil. In addition, for more general situations the linear regression approach has been proposed and validated.

Semi-experimental equilibrium structure determinations by employing B3LYP/SNSD anharmonic force fields: validation and application to semirigid organic molecules.

This work aims at extending the semi-experimental (SE) approach for deriving accurate equilibrium structures to large molecular systems of organic and biological interest. SE equilibrium structures are derived by a least-squares fit of the structural parameters to the experimental ground-state rotational constants of several isotopic species corrected by vibrational contributions computed by quantum mechanical (QM) methods. A systematic benchmark study on 21 small molecules (CCse set) is carried out to evaluate the performance of hybrid density functionals (in particular B3LYP) in the derivation of vibrational corrections to rotational constants. The resulting SE equilibrium structures show a very good agreement with the corresponding geometries obtained employing post-Hartree-Fock vibrational corrections. The use of B3LYP in conjunction with the double-ζ SNSD basis set strongly reduces the computational costs, thus allowing for the evaluation of accurate SE equilibrium structures for medium-sized molecular systems. On these grounds, an additional set of 26 SE equilibrium structures including the most common organic moieties has been set up by collecting the most accurate geometries available in the literature together with new determinations from the present work. The overall set of 47 SE equilibrium structures determined using B3LYP/SNSD vibrational corrections (B3se set) provides a high quality benchmark for validating the structural predictions of other experimental and/or computational approaches. Finally, we present a new strategy (referred to as the template approach) to deal with the cases for which it is not possible to fit all geometrical parameters due to the lack of experimental data.

Molecular structure and spectroscopic signatures of acrolein: theory meets experiment.

A comprehensive study of the molecular structure and IR spectrum of cis and trans acrolein has been performed by an integrated computational approach coupling methods rooted in the coupled-cluster ansatz and the density functional theory. From the one side, DFT anharmonic force fields allow us to determine very reliable semiexperimental structures for both isomers, which are in remarkable agreement with the geometries issuing from CCSD(T) computations accounting for the extrapolation to the complete basis set and core correlation. The same kind of coupled-cluster computations provide dipole moment, relative energies, and interconversion barrier in remarkable agreement with experiments. Finally, harmonic CCSD(T) results coupled to DFT evaluation of mechanical and electrical anharmonicity allow us, in the framework of second-order perturbative vibrational theory, to confirm most of the experimental assignments of IR spectra, and to suggest some additional interpretations for congested regions including fundamental bands together with overtones and combination bands.