FROG: Exploiting all-atom molecular dynamics trajectories to calculate linear and non-linear optical responses of molecular liquids within Dalton's QM/MM polarizable embedding scheme.

Quantum mechanical/molecular mechanics (QM/MM) methods are interesting to model the impact of a complex environment on the spectroscopic properties of a molecule. In this context, a FROm molecular dynamics to second harmonic Generation (FROG) code is a tool to exploit molecular dynamics trajectories to perform QM/MM calculations of molecular optical properties. FROG stands for "FROm molecular dynamics to second harmonic Generation" since it was developed for the calculations of hyperpolarizabilities. These are relevant to model non-linear optical intensities and compare them with those obtained from second harmonic scattering or second harmonic generation experiments. FROG's specificity is that it is designed to study simple molecular liquids, including solvents or mixtures, from the bulk to the surface. For the QM/MM calculations, FROG relies on the Dalton package: its electronic-structure models, response theory, and polarizable embedding schemes. FROG helps with the global workflow needed to deal with numerous QM/MM calculations: it permits the user to separate the system into QM and MM fragments, to write Dalton's inputs, to manage the submission of QM/MM calculations, to check whether Dalton's calculation finished successfully, and finally to perform averages on relevant QM observables. All molecules within the simulation box and several time steps are tackled within the same workflow. The platform is written in Python and installed as a package. Intermediate data such as local electric fields or individual molecular properties are accessible to the users in the form of Python object arrays. The resulting data are easily extracted, analyzed, and visualized using Python scripts that are provided in tutorials.

Liquid Water: When Hyperpolarizability Fluctuations Boost and Reshape the Second Harmonic Scattering Intensities.

Second harmonic scattering (SHS) is a method of choice to investigate the molecular structure of liquids. While a clear interpretation of SHS intensity exists for diluted solutions of dyes, the scattering due to solvents remains difficult to interpret quantitatively. Here, we report a quantum mechanics/molecular mechanics (QM/MM) approach to model the polarization-resolved SHS intensity of liquid water, quantifying different contributions to the signal. We point out that the molecular hyperpolarizability fluctuations and correlations cannot be neglected. The intermolecular orientational and hyperpolarizability correlations up to the third solvation layer strongly increase the scattering intensities and modulate the polarization-resolved oscillation that is predicted here by QM/MM without fitting parameters. Our approach can be generalized to other pure liquids to provide a quantitative interpretation of SHS intensities in terms of short-range molecular ordering.

First hyperpolarizability of water in bulk liquid phase: long-range electrostatic effects included via the second hyperpolarizability.

The molecular first hyperpolarizability ß contributes to second-order optical non-linear signals collected from molecular liquids. For the Second Harmonic Generation (SHG) response, the first hyperpolarizability ß(2ω, ω, ω) often depends on the molecular electrostatic environment. This is especially true for water, due to its large second hyperpolarizability γ(2ω, ω, ω,0). In this study we compute the electronic γ(2ω, ω, ω,0) and ß(2ω, ω, ω) for water molecules in liquid water using QM/MM calculations. The average value of γ(2ω, ω, ω,0) is smaller than the one for the gaz phase, and its standard deviation among the molecules is relatively small. In addition, we demonstrate that the average bulk second hyperpolarizability 〈γ(2ω, ω, ω,0)〉 can be used to describe the electrostatic effects of the distant neighborhood on the first hyperpolarizability ß(2ω, ω, ω). In comparison with more complex schemes to take into account long-range effects, the approximation is simple, and does not require any modifications of the QM/MM implementation. The long-range correction can be added explicitly, using an average value of γ for water in the condensed phase. It can also be easily added implicitly in QM/MM calculations through an additional embedding electric field, without the explicit calculation of γ.

First hyperpolarizability of water at the air-vapor interface: a QM/MM study questions standard experimental approximations.

Surface Second-Harmonic Generation (S-SHG) experiments provide a unique approach to probe interfaces. One important issue for S-SHG is how to interpret the S-SHG intensities at the molecular level. Established frameworks commonly assume that each molecule emits light according to an average molecular hyperpolarizability tensor ß(-2ω,ω,ω). However, for water molecules, this first hyperpolarizability is known to be extremely sensitive to their environment. We have investigated the molecular first hyperpolarizability of water molecules within the liquid-vapor interface, using a quantum description with explicit, inhomogeneous electrostatic embedding. The resulting average molecular first hyperpolarizability tensor depends on the distance relative to the interface, and it practically respects the Kleinman symmetry everywhere in the liquid. Within this numerical approach, based on the dipolar approximation, the water layer contributing to the Surface Second Harmonic Generation (S-SHG) intensity is less than a nanometer. The results reported here question standard interpretations based on a single, averaged hyperpolarizability for all molecules at the interface. Not only the molecular first hyperpolarizability tensor significantly depends on the distance relative to the interface, but it is also correlated to the molecular orientation. Such hyperpolarizability fluctuations may impact the S-SHG intensity emitted by an aqueous interface.

Electrokinetic sweeping of colloids at a reactive magnesium oxide interface.

Investigating the electrokinetic (EK) response in the vicinity of interfaces has regained interest due to the development of new membrane based processes for energy harvesting or soil depollution. However, the case of reactive interfaces, ubiquitous in these processes, remains scarcely explored. Here we experimentally investigate the EK response of a model interface between an aqueous electrolyte and a bulk MgO crystal surface (100), for different pH. For that purpose, we use a lab-scale non invasive method to monitor the zeta potential of the interface versus time, by confocal fluorescent particle tracking. An unexpected motion of the particles, repelled and then attracted again by the interface is observed. We attributed this motion to the surface reactivity, inducing ion concentration gradients perpendicular to the interface and subsequent diffusiophoresis of the charged particle. Accordingly, we could describe at a semi-quantitative level the particle dynamics by solving numerically the Poisson-Nernst-Planck equations to establish concentration profile in the system and subsequent diffusiophoretic motion. These experiments open the way to the characterization of both the EK response and the reaction rate in the vicinity of reactive interfaces.

Longitudinal position dependence of the second-harmonic generation of optically trapped silica microspheres.

We report the design of a setup combining the simultaneous and independent optical trapping and second-harmonic generation (SHG) of 1 µm diameter silica microspheres with two independent laser beams. Optical trapping is achieved with a tightly focused continuous wave infrared laser beam whereas the SHG intensity from the trapped microparticles is obtained with a 810 nm femtosecond pulsed laser. The silica microparticles are dispersed in an aqueous solution, and a microfluidic channel flow is used to remove untrapped microparticles. We show by the perpendicular displacement of the optical trap from the microfluidic channel wall that it is possible to control the contribution of the channel wall/solution interface to the overall SHG intensity. Stable trapping and SHG detection of two microparticles is also demonstrated. Combining the independent trapping of centrosymmetrical silica microparticles with SHG offers new avenues for analytical studies like surface sensing or all-optical devices where the SHG intensity is controlled by the trapping beam.

Electroosmosis near surfactant laden liquid-air interfaces.

Generation of an electroosmostic (EO) flow near a liquid-gas interface covered with ionic surfactants is experimentally investigated. A combination of microscopic flow measurements with a molecular characterization of the interface by second harmonic generation (SHG) shows that under an electrical forcing, although a liquid flow is generated below the free surface, the surfactants remain immobile. The zeta potential was then determined and compared to the surfactant surface coverage. This combination of experimental techniques opens the route to simultaneously probe the liquid flow near a soapy interface and the corresponding surfactant repartition affecting the hydrodynamic boundary condition.

Electrokinetic transport in liquid foams.

Investigating electrokinetic transport in a liquid foam is at the confluence of two well developed research areas. On one hand, the study of electrokinetic flows (i.e. surface-driven flows generated close to a charged interface) is fairly well understood in regards the solid/liquid interface. On the other hand, the flow of liquid in a 3D deformable network, i.e a foam, under a volume force such as gravity has been thoroughly studied over the past decade. The overlapping zone of these two frameworks is of great interest for both communities as it gives rise to challenging new questions such as: what is the importance of the nature of the charged interface, created by mobile and soluble surfactants in the case of foam, on electrokinetic transport? How does a foam behave when submitted to a surface-driven flow? Can we compensate a volume-driven flow, i.e. gravity, by a surface-driven flow, i.e. electroosmosis? In this review, we will explore these questions on three different scales: a surfactant laden interface, a foam film and a macroscopic foam.

Surface conductivity measurements in nanometric to micrometric foam films.

Foam films (a liquid lamella in air covered by surfactants) are tools of choice for nanofluidic characterization as they are intrinsically nanometric. Their size is indeed fixed by a balance between external pressure and particular molecular interactions in the vicinity of interfaces. To probe the exact nature of these interfaces, different characterizations can be performed. Among them, conductivity in confined systems is a direct probe of the electrostatic environment in the vicinity of the surface. Therefore, we designed a dedicated experiment to measure this conductivity in a cylindrical bubble coupled to interferometry for film thickness characterization. We then show that this conductivity depends on the surfactant nature. These conductivity measurements have been performed in an extremely confined system, the so called Newton black foam films. Unexpectedly in this case, a conductivity close to surface conductivity is recovered.

Osteoid Osteoma with a Multicentric Nidus: Interstitial Laser Ablation under MRI Guidance.

Osteoid osteoma (OO) is a common benign tumor of the bone and is typically treated by thermal ablation with computed tomography (CT) guidance. Only a few cases of multicentric OO have been described. We here report the case of an 11-year-old boy with multicentric OO of the right femur treated with laser ablation under open high-field MRI guidance. The steps of the interventional MRI procedure are described, discussing the benefits and disadvantages of MRI versus CT guidance especially with regard to younger patients.

Soft nanofluidic transport in a soap film.

We investigate experimentally the electrokinetic properties of soft nanofluidic channels that consist in soap films with nanometric thickness, covered with charged surfactants. Both the electric and fluidic responses of the system are measured under an applied voltage drop along the film. The electric field is shown to induce an electro-osmotic hydrodynamic flow in the film. However, in contrast to systems confined between solid surfaces, the soft nature of the nanochannel results furthermore in a thickening of the film. This effect accordingly increases the total electro-osmotic flow rate, which behaves nonlinearly with the applied electric field. This behavior is rationalized in terms of an analogy with a Landau-Levich film withdrawn from a reservoir, with the driving velocity identified here with the electro-osmotic one.

Elastic instability in stratified core annular flow.

We study experimentally the interfacial instability between a layer of dilute polymer solution and water flowing in a thin capillary. The use of microfluidic devices allows us to observe and quantify in great detail the features of the flow. At low velocities, the flow takes the form of a straight jet, while at high velocities, steady or advected wavy jets are produced. We demonstrate that the transition between these flow regimes is purely elastic--it is caused by the viscoelasticity of the polymer solution only. The linear stability analysis of the flow in the short-wave approximation supplemented with a kinematic criterion captures quantitatively the flow diagram. Surprisingly, unstable flows are observed for strong velocities, whereas convected flows are observed for low velocities. We demonstrate that this instability can be used to measure the rheological properties of dilute polymer solutions that are difficult to assess otherwise.