Accurate Quantum-Mechanically Derived Force-Fields through a Fragment-Based Approach: Balancing Specificity and Transferability in the Prediction of Self-Assembly in Soft Matter

The wide range of time/length scales covered by self-assembly in soft matter makes molecular dynamics (MD) the ideal candidate for simulating such a supramolecular phenomenon at an atomistic level. However, the reliability of MD outcomes heavily relies on the accuracy of the adopted force-field (FF). The spontaneous re-ordering in liquid crystalline materials stands as a clear example of such collective self-assembling processes, driven by a subtle and delicate balance between supramolecular interactions and single-molecule flexibility. General-purpose transferable FFs often dramatically fail to reproduce such complex phenomena, for example, the error on the transition temperatures being larger than 100 K. Conversely, quantum-mechanically derived force-fields (QMD-FFs), specifically tailored for the target system, were recently shown (J. Phys. Chem. Lett.2022,13, 243) to allow for the required accuracy as they not only well reproduced transition temperatures but also yielded a quantitative agreement with the experiment on a wealth of structural, dynamic, and thermodynamic properties. The main drawback of this strategy stands in the computational burden connected to the numerous quantum mechanical (QM) calculations usually required for a target-specific parameterization, which has undoubtedly hampered the routine application of QMD-FFs. In this work, we propose a fragment-based strategy to extend the applicability of QMD-FFs, in which the amount of QM calculations is significantly reduced, being a single-molecule-optimized geometry and its Hessian matrix the only QM information required. To validate this route, a new FF is assembled for a large mesogen, exploiting the parameters obtained for two smaller liquid crystalline molecules, in this and previous work. Lengthy MD simulations are carried out with the new transferred QMD-FF, observing again a spontaneous re-orientation in the correct range of temperatures, with good agreement with the available experimental measures. The present results strongly suggest that a partial transfer of QMD-FF parameters can be invoked without a significant loss of accuracy, thus paving the way to exploit the method's intrinsic predictive capabilities in the simulation of novel soft materials.

Complexes of Alkaline and Ammonium Cations with Dopamine and Eumelanin Precursors: Dissecting the Role of Noncovalent Cation-π and Cation-Lone Pair (σ-Type) Interactions

Cation-π interactions and their possible competition with other noncovalent interactions (NCI) might play a key role in both dopamine- and eumelanin-based bioinspired materials. In this contribution, to unravel the delicate interplay between cation-π interactions and other possible competing forces, the configurational space of noncovalent complexes formed by dopamine or eumelanin precursors (o-benzoquinone, DHI and a semiquinone dimer) and three different cations (Na⁺, K⁺, and NH₄⁺) is sampled by means of accurate ab initio calculations. To this end, we resort to the mp2^mod method, recently validated by us for benzene-, phenol-, and catechol-cation complexes, whose computational convenience allows for an extensive exploration of the cation-molecule interaction energy surface, by sampling a total of more than 10⁴ arrangements. The mp2^mod interaction energy landscapes reveal that, besides the expected cation-π driven arrangements, for all considered molecule-cation pairs the most stable complexes are found when the cation lies within the plane containing the six-membered ring, thus maximizing the σ-type interaction with the oxygen's lone pairs. Due to the loss of aromaticity, the σ-type/cation-π strength ratio is remarkably large in o-benzoquinone, where cation-π complexes seem unlikely to be formed. The above features are shared among all considered cations but are significantly larger when considering the smaller Na⁺. Besides delivering a deeper insight onto the NCI network established by the considered precursors in the presence of ions, the present results can serve as a reference database to validate or refine lower level methods, as, for instance, the force fields employed in classical simulations.

Automated Parameterization of Quantum Mechanically Derived Force Fields for Soft Materials and Complex Fluids: Development and Validation

The reliability of molecular dynamics (MD) simulations in predicting macroscopic properties of complex fluids and soft materials, such as liquid crystals, colloidal suspensions, or polymers, relies on the accuracy of the adopted force field (FF). We present an automated protocol to derive specific and accurate FFs, fully based on ab initio quantum mechanical (QM) data. The integration of the Joyce and Picky procedures, recently proposed by our group to provide an accurate description of simple liquids, is here extended to larger molecules, capable of exhibiting more complex fluid phases. While the standard Joyce protocol is employed to parameterize the intramolecular FF term, a new automated procedure is here proposed to handle the computational cost of the QM calculations required for the parameterization of the intermolecular FF term. The latter is thus obtained by integrating the old Picky procedure with a fragmentation reconstruction method (FRM) that allows for a reliable, yet computationally feasible sampling of the intermolecular energy surface at the QM level. The whole FF parameterization protocol is tested on a benchmark liquid crystal, and the performances of the resulting quantum mechanically derived (QMD) FF were compared with those delivered by a general-purpose, transferable one, and by the third, "hybrid" FF, where only the bonded terms were refined against QM data. Lengthy atomistic MD simulations are carried out with each FF on extended 5CB systems in both isotropic and nematic phases, eventually validating the proposed protocol by comparing the resulting macroscopic properties with other computational models and with experiments. The QMD-FF yields the best performances, reproducing both phases in the correct range of temperatures and well describing their structure, dynamics, and thermodynamic properties, thus providing a clear protocol that may be explored to predict such properties on other complex fluids or soft materials.

The phenoxyl group-modulated interplay of cation-π and σ-type interactions in the alkali metal series

The interaction potential energy surfaces (IPESs) of four alkaline metal cations (Na+, K+, Rb+ and Cs+) complexed with phenol and catechol were explored by accurate ab initio calculations to investigate the interplay of different noncovalent interactions and their behavior along the alkali metal series and upon -OH substitution. Selected one-dimensional interaction energy curves revealed two different minimum energy configurations for all phenol- and catechol-metal complexes, characterized either by cation-π or σ-type interactions. For each investigated complex several two-dimensional IPES maps were also computed, exploiting the computational advantages of the MP2mod approach. The size of the alkali cation was found to play a similar role in modulating both kinds of complexes, as the interaction strength always decreases along the metal series, from Na+ to Cs+. Conversely, the number of hydroxyl substituents markedly affected cation-π complexes vs. σ-type ones. As a most relevant finding, in catechol-metal complexes the strength of cation-π interactions is around half that of the σ-type ones. It is argued that the combined effect of cation dimensions and hydroxyl substitution in catechol-Na+ complexes makes σ-type configurations remarkably more stable and easily accessible than cation-π ones. Besides shedding new light on the origin of biological phenomena connected with underwater adhesion, the quantum mechanical interaction energy database provided herein may offer a useful reference for tuning accurate force fields, suitable for molecular dynamics simulations, where environmental effects might be also taken into account.

Total Synthesis of (±)-Liphagal via Organic-Redox-Driven Palladium-Catalyzed Hydroxybenzofuran Formation

A synthetic route to liphagal, a natural PI3Kα inhibitor isolated from Aka coralliphaga, was established. The present route features an organic redox process where an alkynylquinone undergoes reductive cyclization in the presence of a hydroquinone derivative such as hydroxyquinol (1,2,4-benzenetriol) and catalytic PdCl₂ to provide a substituted benzofuran suitable for accessing the natural product. The benzofuran formation takes place via the redox transformation between the alkynylquinone and the electron-rich hydroquinones followed by the concomitant Pd(II)-catalyzed oxycyclization of the resultant alkynylhydroquinone.

Emerging investigator series: critical review of photophysical models for the optical and photochemical properties of dissolved organic matter

Optical measurements (absorbance and fluorescence) are widely used to track dissolved organic matter (DOM) quantity and quality in natural and engineered systems. Despite many decades of research on the optical properties of DOM, there is a lack of understanding with regards to the underlying photophysical model that is the basis for these optical properties. This review both summarizes advances to date on the photophysical properties of DOM and seeks to critically evaluate the photophysical models for DOM optical properties. Recent studies have refined the quantitative understanding of DOM photophysical properties such as excited state lifetimes and energies, rates of different photophysical processes, and quantum yields. Considering fundamental models, more clarity is needed on whether DOM photophysical processes are due to a superposition of non-interacting components (superposition model), or whether a portion of optical signals can be ascribed to electronically interacting moieties, for example in the form of electron donor-acceptor complexes (charge transfer model). Multiple studies over more than two decades have provided evidence for the charge transfer model. Questions have been raised, however, about the broad applicability of the charge transfer model. The charge transfer and superposition model are critically reviewed in light of this current research. Recommendations are given for future studies to help clarify the accuracy of these competing photophysical models.

Redox-Active Guanidines in Proton-Coupled Electron-Transfer Reactions: Real Alternatives to Benzoquinones?

Guanidino-functionalized aromatics (GFAs) are readily available, stable organic redox-active compounds. In this work we apply one particular GFA compound, 1,2,4,5-tetrakis(tetramethylguanidino)benzene, in its oxidized form in a variety of oxidation/oxidative coupling reactions to demonstrate the scope of its proton-coupled electron transfer (PCET) reactivity. Addition of an excess of acid boosts its oxidation power, enabling the oxidative coupling of substrates with redox potentials of at least +0.77 V vs. Fc⁺ /Fc. The green recyclability by catalytic re-oxidation with dioxygen is also shown. Finally, a direct comparison indicates that GFAs are real alternatives to toxic halo- or cyano-substituted benzoquinones.

UV-Vis spectrophotometry of quinone flow battery electrolyte for in situ monitoring and improved electrochemical modeling of potential and quinhydrone formation

Quinone-based aqueous flow batteries provide a potential opportunity for large-scale, low-cost energy storage due to their composition from earth abundant elements, high aqueous solubility, reversible redox kinetics and their chemical tunability such as reduction potential. In an operating flow battery utilizing 9,10-anthraquinone-2,7-disulfonic acid, the aggregation of an oxidized quinone and a reduced hydroquinone to form a quinhydrone dimer causes significant variations from ideal solution behavior and of optical absorption from the Beer-Lambert law. We utilize in situ UV-Vis spectrophotometry to establish (a), quinone, hydroquinone and quinhydrone molar attenuation profiles and (b), an equilibrium constant for formation of the quinhydrone dimer (K_QHQ) ∼ 80 M^-1. We use the molar optical attenuation profiles to identify the total molecular concentration and state of charge at arbitrary mixtures of quinone and hydroquinone. We report density functional theory calculations to support the quinhydrone UV-Vis measurements and to provide insight into the dimerization conformations. We instrument a quinone-bromine flow battery with a Pd-H reference electrode in order to demonstrate how complexation in both the negative (quinone) and positive (bromine) electrolytes directly impacts measured half-cell and full-cell voltages. This work shows how accounting for electrolyte complexation improves the accuracy of electrochemical modeling of flow battery electrolytes.

Noncovalent Interactions in the Catechol Dimer

Noncovalent interactions play a significant role in a wide variety of biological processes and bio-inspired species. It is, therefore, important to have at hand suitable computational methods for their investigation. In this paper, we report on the contribution of dispersion and hydrogen bonds in both stacked and T-shaped catechol dimers, with the aim of delineating the respective role of these classes of interactions in determining the most stable structure. By using second-order Møller⁻Plesset (MP2) calculations with a small basis set, specifically optimized for these species, we have explored a number of significant sections of the interaction potential energy surface and found the most stable structures for the dimer, in good agreement with the highly accurate, but computationally more expensive coupled cluster single and double excitation and the perturbative triples (CCSD(T))/CBS) method.

Eumelanin broadband absorption develops from aggregation-modulated chromophore interactions under structural and redox control

Eumelanins, the chief photoprotective pigments in man and mammals, owe their black color to an unusual broadband absorption spectrum whose origin is still a conundrum. Excitonic effects from the interplay of geometric order and disorder in 5,6-dihydroxyindole (DHI)-based oligomeric/polymeric structures play a central role, however the contributions of structural (scaffold-controlled) and redox (π-electron-controlled) disorder have remained uncharted. Herein, we report an integrated experimental-theoretical entry to eumelanin chromophore dynamics based on poly(vinyl alcohol)-controlled polymerization of a large set of 5,6-dihydroxyindoles and related dimers. The results a) uncover the impact of the structural scaffold on eumelanin optical properties, disproving the widespread assumption of a universal monotonic chromophore; b) delineate eumelanin chromophore buildup as a three-step dynamic process involving the rapid generation of oxidized oligomers, termed melanochromes (phase I), followed by a slow oxidant-independent band broadening (phase II) leading eventually to scattering (phase III); c) point to a slow reorganization-stabilization of melanochromes via intermolecular redox interactions as the main determinant of visible broadband absorption.

Experimental and Theoretical Study of the High-Temperature UV-Visible Spectra of Aqueous Hydroquinone and 1,4-Benzoquinone

UV-visible spectroscopic studies of aqueous hydroquinone (HQ) and 1,4-benzoquinone (BQ) have been carried out along with classical molecular dynamics (MD) and quantum calculations. The experimental results confirmed that HQ is stable in hot compressed water up to at least 523 K at 70 bar, but BQ decomposes at temperatures lower than 373 K, leading to the formation of HQ and other nonabsorbing products. Even though benzoquinone is not stable, our study significantly extended the temperature range of other spectroscopic studies, and the spectra of HQ up to 523 K can still be useful for other studies, particularly those related to organic species in deep ocean hydrothermal vents. Classical MD simulations at high temperatures show, as expected, a weakening of the solute-solvent H-bonding interactions. The dependence of the maximum absorption of BQ on temperature was also analyzed, although a significant degree of decomposition was observed in the time frame of our experiments. The shift of the maximum absorption peak of BQ with temperature was consistent with time-dependent density functional theory calculations.

Electronic State-Resolved Electron-Phonon Coupling in an Organic Charge Transfer Material from Broadband Quantum Beat Spectroscopy

The coupling of electron and lattice phonon motion plays a fundamental role in the properties of functional organic charge-transfer materials. In this Letter we extend the use of ultrafast vibrational quantum beat spectroscopy to directly elucidate electron-phonon coupling in an organic charge-transfer material. As a case study, we compare the oscillatory components of the transient reflection (TR) of a broadband probe pulse from single crystals of quinhydrone, a 1:1 cocrystal of hydroquinone and p-benzoquinone, after exciting nonresonant impulsive stimulated Raman scattering and resonant electronic transitions using ultrafast pulses. Spontaneous resonance Raman spectra confirm the assignment of these oscillations as coherent lattice phonon excitations. Fourier transforms of the vibrational quantum beats in our broadband TR measurements allow construction of spectra that we show report the ability of these phonons to directly modulate the electronic structure of quinhydrone. These results demonstrate how coherent ultrafast processes can characterize the complex interplay of charge transfer and lattice motion in materials of fundamental relevance to chemistry, materials sciences, and condensed matter physics.

Intermolecular interactions in eumelanins: a computational bottom-up approach. I. small building blocks

Building eumelanin: from basic units to spectral properties.