No Arabic abstract
We present the extension of the quantum/classical polarizable fluctuating charge model to the calculation of single residues of quadratic response functions, as required for the computational modeling of two-photon absorption cross-sections. By virtue of a variational formulation of the quantum/classical polarizable coupling, we are able to exploit an atomic orbital-based quasienergy formalism to derive the additional coupling terms in the response equations. Our formalism can be extended to the calculation of arbitrary order response functions and their residues. The approach has been applied to the challenging problem of one- and two-photon spectra of rhodamine 6G (R6G) in aqueous solution. Solvent effects on one- and two-photon spectra of R6G in aqueous solution have been analyzed by considering three different approaches, from a continuum (QM/PCM) to two QM/MM models (non-polarizable QM/TIP3P and polarizable QM/FQ). Both QM/TIP3P and QM/FQ simulated OPA and TPA spectra show that the inclusion of discrete water solvent molecules is essential to increase the agreement between theory and experiment. QM/FQ has been shown to give the best agreement with experiments.
A method for study of charge-transfer interactions between solute molecules and solvent based on the comparison of the ratios of spectral shifts of different electronic transitions in solute molecules in chemically inert solvent is proposed. The method is applicable to molecules that do not change their dipole moment on excitation. As an example, a presence of charge transfer interactions in higher electronic states of aromatic hydrocarbons (benzene, phenanthrene, and naphthalene) dissolved in water and alcohols was demonstrated.
Entangled photon pairs have been promised to deliver a substantial quantum advantage for two-photon absorption spectroscopy. However, recent work has challenged the previously reported magnitude of quantum enhancement in two-photon absorption. Here, we present an experimental comparison of sum-frequency generation and molecular absorption, each driven by isolated photon pairs. We establish an upper bound on the enhancement for entangled-two-photon absorption in Rhodamine-6G, which lies well below previously reported values.
Entangled two-photon absorption (ETPA) has recently become a topic of lively debate, mainly due to the apparent inconsistencies in the experimentally-reported ETPA cross sections of organic molecules. In this work, we provide a thorough experimental study of ETPA in the organic molecules Rhodamine B (RhB) and Zinc Tetraphenylporphirin (ZnTPP). The goal of this contribution is twofold: on one hand, it seeks to reproduce the results of previous experimental reports and, on the other, it aims to determine the effects of different temporal correlations -- introduced as a controllable time-delay between the photons to be absorbed -- on the strength of the ETPA signal. In our experiment, the samples are excited by entangled pairs produced by type-I SPDC, with a spectral distribution centered at 810 nm. Surprisingly, and contrary to what was expected, the time delay did not produce in our experiments any systematic change in the cross-sections when monitoring the ETPA signal using a transmission measurement scheme. As a plausible cause of this unexpected result, we argue that the photon-pair flux, typically-used in these experiments, is not sufficient to promote the two-photon absorption process in these molecules. This suggests that the actual absorption cross-section values are lower than those previously reported, and therefore do not lead to a measurable ETPA effect for the transmission method.
Thermal gradients induce concentration gradients in alkali halide solutions, and the salt migrates towards hot or cold regions depending on the average temperature of the solution. This effect has been interpreted using the heat of transport, which provides a route to rationalize thermophoretic phenomena. Early theories provide estimates of the heat of transport at infinite dilution. These values are used to interpret thermodiffusion (Soret) and thermoelectric (Seebeck) effects. However, accessing heats of transport of individual ions at finite concentration remains an outstanding question both theoretically and experimentally. Here we discuss a computational approach to calculate heats of transport of aqueous solutions at finite concentrations, and apply our method to study lithium chloride solutions at concentrations $>0.5$~M. The heats of transport are significantly different for Li$^+$ and Cl$^-$ ions, unlike what is expected at infinite dilution. We find theoretical evidence for the existence of minima in the Soret coefficient of LiCl, where the magnitude of the heat of transport is maximized. The Seebeck coefficient obtained from the ionic heats of transport varies significantly with temperature and concentration. We identify thermodynamic conditions leading to a maximization of the thermoelectric response of aqueous solutions.
Molecular dynamics simulations of biomolecules have been widely adopted in biomedical studies. As classical point-charge models continue to be used in routine biomolecular applications, there have been growing demands on developing polarizable force fields for handling more complicated biomolecular processes. Here we focus on a recently proposed polarizable Gaussian Multipole (pGM) model for biomolecular simulations. A key benefit of pGM is its screening of all short-range electrostatic interactions in a physically consistent manner, which is critical for stable charge-fitting and is needed to reproduce molecular anisotropy. Another advantage of pGM is that each atoms multipoles are represented by a single Gaussian function or its derivatives, allowing for more efficient electrostatics than other Gaussian-based models. In this study we present an efficient formulation for the pGM model defined with respect to a local frame formed with a set of covalent basis vectors. The covalent basis vectors are chosen to be along each atoms covalent bonding directions. The new local frame allows molecular flexibility during molecular simulations and facilitates an efficient formulation of analytical electrostatic forces without explicit torque computation. Subsequent numerical tests show that analytical atomic forces agree excellently with numerical finite-difference forces for the tested system. Finally, the new pGM electrostatics algorithm is interfaced with the PME implementation in Amber for molecular simulations under the periodic boundary conditions. To validate the overall pGM/PME electrostatics, we conducted an NVE simulation for a small water box of 512 water molecules. Our results show that, to achieve energy conservation in the polarizable model, it is important to ensure enough accuracy on both PME and induction iteration.