Search

We explore the application of an extrapolative method that yields very accurate total and relative energies from variational and diffusion quantum Monte Carlo (VMC and DMC) results. For a trial wave function consisting of a small configuration interaction (CI) wave function obtained from full CI quantum Monte Carlo and reoptimized in the presence of a Jastrow factor and an optional backflow transformation, we find that the VMC and DMC energies are smooth functions of the sum of the squared coefficients of the initial CI wave function and that quadratic extrapolations of the non-backflow VMC and backflow DMC energies intersect within uncertainty of the exact total energy. With adequate statistical treatment of quasi-random fluctuations, the extrapolate and intersect with polynomials of order two method is shown to yield results in agreement with benchmark-quality total and relative energies for the C2, N2, CO2, and H2O molecules, as well as for the C2 molecule in its first electronic singlet excited state, using only small CI expansion sizes. [ABSTRACT FROM AUTHOR]

The potential energy landscape (PEL) formalism is a tool within statistical mechanics that has been used in the past to calculate the equation of states (EOS) of classical rigid model liquids at low temperatures, where computer simulations may be challenging. In this work, we use classical molecular dynamics (MD) simulations and the PEL formalism to calculate the EOS of the flexible q-TIP4P/F water model. This model exhibits a liquid–liquid critical point (LLCP) in the supercooled regime, at (Pc = 150 MPa, Tc = 190 K, and ρc = 1.04 g/cm3) [using the reaction field technique]. The PEL-EOS of q-TIP4P/F water and the corresponding location of the LLCP are in very good agreement with the MD simulations. We show that the PEL of q-TIP4P/F water is Gaussian, which allows us to calculate the configurational entropy of the system, Sconf. The Sconf of q-TIP4P/F water is surprisingly similar to that reported previously for rigid water models, suggesting that intramolecular flexibility does not necessarily add roughness to the PEL. We also show that the Adam–Gibbs relation, which relates the diffusion coefficient D with Sconf, holds for the flexible q-TIP4P/F water model. Overall, our results indicate that the PEL formalism can be used to study molecular systems that include molecular flexibility, the common case in standard force fields. This is not trivial since the introduction of large bending/stretching mode frequencies is problematic in classical statistical mechanics. For example, as shown previously, we find that such high frequencies lead to unphysical (negative) entropy for q-TIP4P/F water when using classical statistical mechanics (yet, the PEL formalism can be applied successfully). [ABSTRACT FROM AUTHOR]

Electronic relaxation dynamics of solution-phase redox-exfoliated molybdenum disulfide (MoS2) monolayer and multilayer ensembles are described. MoS2 was exfoliated using polyoxometalate (POM) reductants. This process yields a colloidal heterostructure consisting of MoS2 2D sheet multilayers with surface-bound POM complexes. Using two-dimensional electronic spectroscopy, transient bleaching and photoinduced absorption signals were detected at excitation/detection energies of 1.82/1.87 and 1.82/1.80 eV, respectively. Approximate 100-fs bandgap renormalization (BGR) and subsequent defect- and phonon-mediated relaxation on the picosecond timescale were resolved for several MoS2 thicknesses spanning from 1 to 2 L to ∼20 L. BGR rates were independent of sample thickness and slightly slower than observations for chemical vapor deposition-grown MoS2 monolayers. However, defect-mediated relaxation accelerated ∼10-fold with increased sample thicknesses. The relaxation rates increased from 0.33 ± 0.05 to 1.2 ± 0.1 and 3.1 ± 0.4 ps−1 for 1–2 L, 3–4 L, and 20 L fractions. The thicknesses-dependent relaxation rates for POM-MoS2 heterostructures were modeled using a saturating exponential function that showed saturation at thirteen MoS2 layers. The results suggest that the increased POM surface coverage leads to larger defect density in the POM-MoS2 heterostructure. These are the first descriptions of the influence of sample thickness on electronic relaxation rates in solution-phase redox-exfoliated POM-MoS2 heterostructures. Outcomes of this work are expected to impact the development of solution-phase exfoliation of 2D metal-chalcogenide heterostructures. [ABSTRACT FROM AUTHOR]

We report threshold collision induced dissociation experiments on protonated water clusters thermalized at low temperature for sizes n = 19–23. Fragmentation cross sections are recorded as a function of the collision energy and analyzed with a statistical model. This model allows us to account for dissociation cascades and provides values for the dissociation energies of each cluster. These values, averaging around 0.47 eV, are in good agreement with theoretical predictions at various levels of theory. Furthermore, the dissociation energies show a trend for the n = 21 magic and n = 22 anti-magic numbers relative to their neighbours, which is also in agreement with theory. These results provide further evidence to resolve the disagreement between previously published experimental values. A careful quantitative treatment of cascade dissociation in this model introduces interdependence between the dissociation energies of neighboring sizes, which reduces the number of free fitting parameters and improves both reliability and uncertainties on absolute dissociation energies deduced from experiments. [ABSTRACT FROM AUTHOR]

Recent vibrational sum frequency generation spectroscopic experiments [Sengupta et al., J. Phys. Chem. Lett. 13, 11391–11397 (2022)] demonstrated synergistic interfacial adsorption effects between the anionic dodecyl sulfate (DS−) and the polar, but charge-neutral hexaethylene glycol monododecyl ether (C12E6), surfactants. In this study, the interfacial adsorption thermodynamics underlying these synergistic effects are analyzed through free energy decompositions. A general decomposition method utilizing alchemical intermediate states is outlined. Combining free energy decompositions with the potential distribution theorem illuminates the statistical interpretations of correlated effects between different system components. This approach allows for the identification of the physical effects leading to synergistic adsorption thermodynamics of DS− binding to the air-C12E6-water interface. The binding properties are found to result from a combination of effects predominantly including energetic van der Waals stabilization between DS− and C12E6, as well as competing energetic and entropic effects due to changes in the interfacial water structure as a result of introducing a C12E6 monolayer into the bare air–water interface. [ABSTRACT FROM AUTHOR]

We investigate the wetting properties of PDMS (Polydimethylsiloxane) pseudo-brush anchored on glass substrates. These PDMS pseudo-brushes exhibit a significantly lower contact angle hysteresis compared to hydrophobic silanized substrates. The effect of different molar masses of the used PDMS on the wetting properties seems negligible. The surface roughness and thickness of the PDMS pseudo-brush are measured by atomic force microscopy and x-ray reflectivity. The outcome shows that these surfaces are extremely smooth (topologically and chemically), which explains the reduction in contact angle hysteresis. These special features make this kind of surfaces very useful for wetting experiments. Here, the dynamics of the four-phase contact point are studied on these surfaces. The four-phase contact point dynamics on PDMS pseudo-brushes deviate substantially from its dynamics on other substrates. These changes depend only a little on the molar mass of the used PDMS. In general, PDMS pseudo-brushes increase the traveling speed of four-phase contact point on the surface and change the associated power law of position vs time. [ABSTRACT FROM AUTHOR]

With a transcorrelated Hamiltonian, we perform a many body perturbation calculation on the uniform electron gas in the high density regime. By using a correlation factor optimized for a single determinant Jastrow ansatz, the second order correlation energy is calculated as 1 − ln ⁡ 2 π 2 ln ( r s ) − 0.050 75. This already reproduces the exact logarithmic term of the random phase approximation (RPA) result, while the constant term is roughly 7% larger than the RPA one. The close agreement with the RPA method demonstrates that the transcorrelated method offers a viable and potentially efficient method for treating metallic systems. [ABSTRACT FROM AUTHOR]

We perform path integral molecular dynamics (PIMD) simulations of a monatomic liquid that exhibits a liquid–liquid phase transition and liquid–liquid critical point. PIMD simulations are performed using different values of Planck's constant h, allowing us to study the behavior of the liquid as nuclear quantum effects (NQE, i.e., atoms delocalization) are introduced, from the classical liquid (h = 0) to increasingly quantum liquids (h > 0). By combining the PIMD simulations with the ring-polymer molecular dynamics method, we also explore the dynamics of the classical and quantum liquids. We find that (i) the glass transition temperature of the low-density liquid (LDL) is anomalous, i.e., T g LDL (P) decreases upon compression. Instead, (ii) the glass transition temperature of the high-density liquid (HDL) is normal, i.e., T g HDL (P) increases upon compression. (iii) NQE shift both T g LDL (P) and T g HDL (P) toward lower temperatures, but NQE are more pronounced on HDL. We also study the glass behavior of the ring-polymer systems associated with the quantum liquids studied (via the path-integral formulation of statistical mechanics). There are two glass states in all the systems studied, low-density amorphous ice (LDA) and high-density amorphous ice (HDA), which are the glass counterparts of LDL and HDL. In all cases, the pressure-induced LDA–HDA transformation is sharp, reminiscent of a first-order phase transition. In the low-quantum regime, the LDA–HDA transformation is reversible, with identical LDA forms before compression and after decompression. However, in the high-quantum regime, the atoms become more delocalized in the final LDA than in the initial LDA, raising questions on the reversibility of the LDA–HDA transformation. [ABSTRACT FROM AUTHOR]

Over the last several decades, a range of experimental techniques from x-ray crystallography and atomic force microscopy to nuclear magnetic resonance and small angle x-ray scattering have probed nucleic acid structure and conformation with high resolution both in the condensed state and in solution. We present a computational study that examines the prospect of using electrostatic free energy measurements to detect 3D conformational properties of nucleic acid molecules in solution. As an example, we consider the conformational difference between A- and B-form double helices whose structures differ in the values of two key parameters—the helical radius and rise per basepair. Mapping the double helix onto a smooth charged cylinder reveals that electrostatic free energies for molecular helices can, indeed, be described by two parameters: the axial charge spacing and the radius of a corresponding equivalent cylinder. We show that electrostatic free energies are also sensitive to the local structure of the molecular interface with the surrounding electrolyte. A free energy measurement accuracy of 1%, achievable using the escape time electrometry (ETe) technique, could be expected to offer a measurement precision on the radius of the double helix of approximately 1 Å. Electrostatic free energy measurements may, therefore, not only provide information on the structure and conformation of biomolecules but could also shed light on the interfacial hydration layer and the size and arrangement of counterions at the molecular interface in solution. [ABSTRACT FROM AUTHOR]

We design a geometric Brownian information engine by considering overdamped Brownian particles inside a two-dimensional monolobal confinement with irregular width along the transport direction. Under such detention, particles experience an effective entropic potential which has a logarithmic form. We employ a feedback control protocol as an outcome of error-free position measurement. The protocol comprises three stages: measurement, feedback, and relaxation. We reposition the center of the confinement to the measurement distance (xp) instantaneously when the position of the trapped particle crosses xp for the first time. Then, the particle is allowed for thermal relaxation. We calculate the extractable work, total information, and unavailable information associated with the feedback control using this equilibrium probability distribution function. We find the exact analytical value of the upper bound of extractable work as ( 5 3 − 2 l n 2) k B T. We introduce a constant force G downward to the transverse coordinate (y). A change in G alters the effective potential of the system and tunes the relative dominance of entropic and energetic contributions in it. The upper bound of the achievable work shows a crossover from ( 5 3 − 2 l n 2) k B T to 1 2 k B T when the system changes from an entropy-dominated regime to an energy-dominated one. Compared to an energetic analog, the loss of information during the relaxation process is higher in the entropy-dominated region, which accredits the less value in achievable work. Theoretical predictions are in good agreement with the Langevin dynamics simulation studies. [ABSTRACT FROM AUTHOR]

We calculate the distribution probability of hexagonal six-site rings in the disordered state of a cubic or hexagonal ice lattice model with ice rules perfectly obeyed. The mean-field distribution obtained is in significant agreement with those, slightly different among them, obtained by Monte Carlo simulations of the cubic or hexagonal model. The results are discussed in connection with the equilibrium and non-equilibrium transition from disorder to ferroelectric proton order. [ABSTRACT FROM AUTHOR]

Graphene is an ideal material for biosensors due to the large surface area for multiple bonding sites, the high electrical conductivity allowing for high sensitivity, and the high tensile strength providing durability in fabricated sensor devices. For graphene to be successful as a biosensing platform, selectivity must be achieved through functionalization with specific chemical groups. However, the device performance and sensor sensitivity must still be maintained after functionalization, which can be challenging. We compare phenyl amine and 1,5-diaminonaphthalene functionalization methods for chemical vapor deposition grown graphene, both used to obtain graphene modified with amine groups—which is required for surface attachment of highly selective antibody bio-receptors. Through atomic force microscopy (AFM), Raman spectroscopy, and time-of-flight secondary ion mass spectrometry imaging of co-located areas, the chemistry, thickness, and coverage of the functional groups bound to the graphene surface have been comprehensively analyzed. We demonstrate the modification of functionalized graphene using AFM, which unexpectedly suggests the removal of covalently bonded functional groups, resulting in a "recovered" graphene structure with reduced disorder, confirmed with Raman spectroscopy. This removal explains the decrease in the ID/IG ratio observed in Raman spectra from other studies on functionalized graphene after mechanical strain or a chemical reaction and reveals the possibility of reverting to the non-functionalized graphene structure. Through this study, preferred functionalization processes are recommended to maintain the performance properties of graphene as a biosensor. [ABSTRACT FROM AUTHOR]

Dielectric interfaces are crucial to the behavior of charged membranes, from graphene to synthetic and biological lipid bilayers. Understanding electrolyte behavior near these interfaces remains a challenge, especially in the case of rough dielectric surfaces. A lack of analytical solutions consigns this problem to numerical treatments. We report an analytic method for determining electrostatic potentials near curved dielectric membranes in a two-dimensional periodic "slab" geometry using a periodic summation of Green's functions. This method is amenable to simulating arbitrary groups of charges near surfaces with two-dimensional deformations. We concentrate on one-dimensional undulations. We show that increasing membrane undulation increases the asymmetry of interfacial charge distributions due to preferential ionic repulsion from troughs. In the limit of thick membranes, we recover results mimicking those for electrolytes near a single interface. Our work demonstrates that rough surfaces generate charge patterns in electrolytes of charged molecules or mixed-valence ions. [ABSTRACT FROM AUTHOR]

We investigate the polarization of disk electrodes immersed in an electrolyte solution and subjected to a small external AC voltage over a wide range of frequencies. A mathematical model is developed based on the Debye–Falkenhagen approximation to the coupled Poisson–Nernst–Planck equations. Analytical techniques are used for predicting the spatial distribution of the electric potential and the complex impedance of the system. Scales for impedance and frequency are identified, which lead to a self-similar behavior for a range of frequencies. Experiments are conducted with gold electrodes of sizes in the range 100–350 μm immersed in a high-conductivity KCl solution over five orders of magnitude in frequency. A collapse of data on impedance magnitude and phase angle onto universal curves is observed with scalings motivated by the mathematical model. A direct comparison with the approximate analytical formula for impedance is made without any fitting parameters, and a good agreement is found for the range of frequencies where the analytical model is valid. [ABSTRACT FROM AUTHOR]

Heat transfer between a silver nanoparticle and surrounding water has been studied using molecular dynamics (MD) simulations. The thermal conductance (Kapitza conductance) at the interface between a nanoparticle and surrounding water has been calculated using four different approaches: transient with/without temperature gradient (internal thermal resistance) in the nanoparticle, steady-state non-equilibrium, and finally equilibrium simulations. The results of steady-state non-equilibrium and equilibrium are in agreement but differ from the transient approach results. MD simulation results also reveal that in the quenching process of a hot silver nanoparticle, heat dissipates into the solvent over a length-scale of ∼2 nm and over a time scale of less than 5 ps. By introducing a continuum solid-like model and considering a heat conduction mechanism in water, it is observed that the results of the temperature distribution for water shells around the nanoparticle agree well with the MD results. It is also found that the local water thermal conductivity around the nanoparticle is greater by about 50% than that of bulk water. These results have important implications for understanding heat transfer mechanisms in nanofluid systems and also for cancer photothermal therapy, wherein an accurate local description of heat transfer in an aqueous environment is crucial. [ABSTRACT FROM AUTHOR]

High-resolution anion photoelectron spectroscopy of the ZrO3H2− and ZrO3D2− anions and complementary electronic structure calculations are used to investigate the reaction between zirconium dioxide and a single water molecule, ZrO20/− + H2O. Experimental spectra of ZrO3H2− and ZrO3D2− were obtained using slow photoelectron velocity-map imaging of cryogenically cooled anions, revealing the presence of two dissociative adduct conformers and yielding insight into the vibronic structure of the corresponding neutral species. Franck–Condon simulations for both the cis- and trans-dihydroxide structures are required to fully reproduce the experimental spectrum. Additionally, it was found that water-splitting is stabilized more by ZrO2 than TiO2, suggesting Zr-based catalysts are more reactive toward hydrolysis. [ABSTRACT FROM AUTHOR]

To investigate the chemical isotope-exchange reactions within a system composed of a mixture of hydrogen and deuterium (H/D) in the plasma media, the ReaxFFHD potential was parameterized against an appropriate quantum mechanics (QM)-based training set. These QM data involve structures and energies related to bond dissociation, angle distortion, and an exchange reaction of the tri-atomic molecular ions, H3+, D3+, H2D+, and D2H+, produced in the hydrogen plasma. Using the ReaxFFHD potential, a range of reactive molecular dynamics simulations were performed on different mixtures of H/D systems. Analysis of the reactions involved in the production of these tri-atomic molecular ions was carried out over 1 ns simulations. The results show that the ReaxFFHD potential can properly model isotope-exchange reactions of tri-atomic molecular ions and that it also has a perfect transferability to reactions taking place in these systems. In our simulations, we observed some intermediate molecules (H2, D2, and HD) that undergo secondary reactions to form the tri-atomic molecular ions as the most likely products in the hydrogen plasma. Moreover, there remains a preference for D in the produced molecular ions, which is related to the lower zero-point energy of the D-enriched species, showing the isotope effects at the heart of the ReaxFFHD potential. [ABSTRACT FROM AUTHOR]

We study the steady-state behavior of active, dipolar, Brownian spheroids in a planar channel subjected to an imposed Couette flow and an external transverse field, applied in the "downward" normal-to-flow direction. The field-induced torque on active spheroids (swimmers) is taken to be of magnetic form by assuming that they have a permanent magnetic dipole moment, pointing along their self-propulsion (swim) direction. Using a continuum approach, we show that a host of behaviors emerges over the parameter space spanned by the particle aspect ratio, self-propulsion and shear/field strengths, and the channel width. The cross-stream migration of the model swimmers is shown to involve a regime of linear response (quantified by a linear-response factor) in weak fields. For prolate swimmers, the weak-field behavior crosses over to a regime of full swimmer migration to the bottom half of the channel in strong fields. For oblate swimmers, a counterintuitive regime of reverse migration arises in intermediate fields, where a macroscopic fraction of swimmers reorient and swim to the top channel half at an acute "upward" angle relative to the field axis. The diverse behaviors reported here are analyzed based on the shear-induced population splitting (bimodality) of the swim orientation, giving two distinct, oppositely polarized, swimmer subpopulations (albeit very differently for prolate/oblate swimmers) in each channel half. In strong fields, swimmers of both types exhibit net upstream currents relative to the laboratory frame. The onsets of full migration and net upstream current depend on the aspect ratio, enabling efficient particle separation strategies in microfluidic setups. [ABSTRACT FROM AUTHOR]

In the framework of thermodynamic perturbation theory (TPT), a new perturbed-chain equation of state (EOS) is presented for a fully flexible Lennard-Jones (LJ) chain system. The EOS is the sum of repulsive and perturbation contributions. The reference term of the EOS is derived based on first- and second-order TPT of Wertheim for the chains interacting with each other through the Weeks–Chandler–Anderson potential model. In order to derive the perturbation term, we have used the radial distribution function of the hard-chain system with a chain range of m = 2–10 and packing fraction range of η = 0.10–0.72, which cover the entire density range from vapor to solid phases. The performance of the EOS is tested against simulation data of the compressibility factor, residual internal energy, and phase equilibrium. A close agreement was observed across all cases. The EOS has three pure component parameters and is able to describe the global vapor–liquid–solid phase diagram of the LJ chain. [ABSTRACT FROM AUTHOR]

Ionizing interactions between charged particles and molecules of biological relevance have attracted considerable interest in the last decade due to its importance in medical radiation therapy. We have previously calculated triply differential cross sections for five biomolecules in collaboration with experimental groups. We used the molecular 3-body distorted wave approximation for these calculations. For ionization of biomolecules, experimentalists are unable to determine the orientation of the molecule at the time of ionization, which means that the calculated cross sections need to be averaged over all molecular orientations. At the time the calculations were performed, it was not numerically feasible for us to perform proper averaging over orientations, so we introduced the orientation averaged molecular orbital approximation to make the calculations possible. We now have the computational capability to properly perform this average, so, here, we present new results with a proper average over orientations and compare with the previous calculations and experiment. Since the original calculations, results from two different distorted-wave models have also been published and the new results will also be compared with those calculations. Overall, the new results are in better agreement with the experiment. [ABSTRACT FROM AUTHOR]

As a classic ferromagnetic material, nickel has been an important research candidate used to study dynamics and interactions of electron, spin, and lattice degrees of freedom. In this study, we specifically chose a thick, 150 nm ferromagnetic nickel (111) single crystal rather than 10–20 nm thin crystals that are typically used in ultrafast studies, and we revealed both the ultrafast heating within the skin depth and the heat transfer from the surface (skin) layer to the bulk of the crystal. The lattice deformation after femtosecond laser excitation was investigated by means of 8.04 keV subpicosecond x-ray pulses, generated from a table-top laser-plasma based source. The temperature evolution of the electron, spin, and lattice was determined using a three temperature model. In addition to coherent phonon oscillations, the blast force and sonic waves, induced by the hot electron temperature gradient, were also observed by monitoring the lattice contractions during the first couple of picoseconds after laser irradiation. This study further revealed the tens of picoseconds time required for heating the hundred nanometer bulk of the Ni (111) single crystals. [ABSTRACT FROM AUTHOR]

We report an experimental and theoretical investigation into the dynamics of electron-impact ionization of R-carvone. Experimental triple differential cross sections are obtained in asymmetric coplanar kinematic conditions for the ionization of the unresolved combination of the three outermost molecular orbitals (41a-39a) of R-carvone. These cross sections are compared with theoretical cross sections calculated within a molecular 3-body distorted wave (M3DW) framework employing either a proper orientation average or orbital average to account for the random orientation of the molecule probed in the experiment. Here, we observe that the overall scattering behavior observed in the experiment is fairly well reproduced within the M3DW framework when implementing the proper average over orientations. The character of the ionized orbitals also provides some qualitative explanation for the observed scattering behavior. This represents substantial progress when trying to describe the scattering dynamics observed for larger molecules under intermediate-impact energy and asymmetric energy sharing scattering conditions. [ABSTRACT FROM AUTHOR]

We report an experimental and theoretical investigation of electron-impact single ionization of the highest occupied molecular orbital 1t2 and the next highest occupied molecular orbital 2a1 states of CH4 at an incident electron energy of 250 eV. Triple differential cross sections measured in two different laboratories were compared with results calculated within the molecular 3-body distorted wave and generalized Sturmian function theoretical models. For ionization of the 1t2 state, the binary peak was observed to have a single maximum near the momentum transfer direction that evolved into a double peak for increasing projectile scattering angles, as has been seen for ionization of atomic p-states. A detailed investigation of this evolution was performed. As expected because of its s-type character, for ionization of the 2a1 state, only a single binary peak was observed. Overall, good agreement was found between experiment and theory. [ABSTRACT FROM AUTHOR]

Poly-glycine-alanine (poly-GA) proteins are widely believed to be one of the main toxic dipeptide repeat molecules associated with amyotrophic lateral sclerosis (ALS) and frontotemporal dementia diseases. Using discontinuous molecular dynamics simulation and an all-atom model of the proteins, we study folding, stability, and aggregation of poly-GA. The results demonstrate that poly-GA is an aggregation-prone protein that, after a long enough time, forms β-sheet-rich aggregates that match recent experiment data and that two unique helical structures are formed very frequently, namely, β-helix and double-helix. The details of the two structures are analyzed. The analysis indicates that such helical structures are stable and share the characteristics of both α-helices and β-sheets. Molecular simulations indicate that identical phenomena also occur in the aggregation of poly-glycine-arginine (poly-GR). Therefore, we hypothesize that proteins of type (GX)n in which X may be any non-glycine amino acid and n is the repeat length may share the same folding structures of β-helix and double-helix and that it is the glycine in the repeat that contributes the most to this characteristic. Molecular dynamics simulation with continuous interaction potentials and explicit water molecules as the solvent supports the hypothesis. To our knowledge, this is the first molecular dynamics simulation of the phenomena involving poly-GA and poly-GR proteins. [ABSTRACT FROM AUTHOR]

We develop a general framework for transport of polyions, solvent and salt, with intended application to Layer-by-Layer (LbL) assembly of polyelectrolyte monolayers (PEMs). The formulation for the first time includes electrostatics, chemical potential gradients, and mechanical stress gradients as driving forces for mass transport. The general model allows all species to be mobile throughout the process and avoids the assumptions of stepwise instantaneous equilibrium and/or immobilized structures typical of previous approaches, while reducing to these models in appropriate limits. A simple constitutive equation is derived for a mixture of oppositely charged polyelectrolytes that accounts for network strand dilution and cross-chain ion pairing by appending reactive terms to the Smoluchowski probability diffusion equation for network strand end-to-end vectors. The resulting general framework encompasses the Poisson equation describing the electrostatic potential distribution, an osmotic pressure balance, a stress constitutive equation, and a generalized flux law of polymer transport. The computational domain is split into a PEM phase and an external solution phase with an appropriate boundary condition derived for the interface between the two. The mobile species (water and small salt ions) are taken to be in a state of dynamic equilibrium with their distributions enslaved to the perturbations in the two polyion compositions. The proposed model captures the swelling response of PEM films to external solutions. For the first time, we studied the effects of the temporal evolution of electrostatic and stress distribution on the rate of chain loss and absorption during rinsing and dipping of an idealized and arbitrarily selected and rigid brush layer into external solutions. The temporal evolution provides a kinetic basis for the ability of LbL films to grow under conditions that thermodynamics alone suggests would cause them to be washed away and to account for partial desorption during washing. The proposed transport framework constitutes a solid basis for eventual quantitative modeling of LbL assembly and transport in polyion networks more generally. [ABSTRACT FROM AUTHOR]

Reversible ion binding equilibria in polyelectrolyte solutions are strongly affected by interactions between dissociated ionic species. We examine how the structural correlations between ionic groups on polyelectrolytes impact the counterion binding. Treating the electrostatic correlation free energy using the classical Debye-Hückel expression leads to complete counterion dissociation in the concentrated regime. This unphysical behavior is shown to stem from improper regularization of the self-energy of dissociated ions and polyions and is mitigated by smearing point-like charges across a finite width. The influence of the self-energy on counterion binding is elaborated on by generalizing the Debye-Hückel free energy to polyelectrolytes with variable fractal dimension and stiffness. In the dilute regime, a greater propensity for binding is found for chains with more compact architectures, which in turn reduces the harsh self-repulsions of tightly packed arrangements of charge. In the concentrated regime, the effects of electrostatic correlations weaken due to screening and the extent of binding is governed by a balance of short-ranged interactions and the translational entropy of ions. [ABSTRACT FROM AUTHOR]

We study the spatial and orientational distribution of charged nanorods (rodlike counterions) as well as the effective interaction mediated by them between two plane-parallel surfaces that carry fixed (quenched) heterogeneous charge distributions. The nanorods are assumed to have an internal charge distribution, specified by a multivalent monopolar moment and a finite quadrupolar moment, and the quenched surface charge is assumed to be randomly distributed with equal mean and variance on the two surfaces. While equally charged surfaces are known to repel within the traditional mean-field theories, the presence of multivalent counterions has been shown to cause attractive interactions between uniformly charged surfaces due to the prevalence of strong electrostatic couplings that grow rapidly with the counterion valency. We show that the combined effects due to electrostatic correlations (caused by the coupling between the mean surface field and the multivalent, monopolar, charge valency of counterions) as well as the disorder-induced interactions (caused by the coupling between the surface disorder field and the quadrupolar moment of counterions) lead to much stronger attractive interactions between two randomly charged surfaces. The interaction profile turns out to be a nonmonotonic function of the intersurface separation, displaying an attractive minimum at relatively small separations, where the ensuing attraction can exceed the maximum strong-coupling attraction (produced by multivalent monopolar counterions between uniformly charged surfaces) by more than an order of magnitude. [ABSTRACT FROM AUTHOR]

Amorphous TiO2 (a-TiO2) could offer an attractive alternative to conventional crystalline TiO2 phases for photocatalytic applications. However, the atomic structure of a-TiO2 remains poorly understood with respect to that of its crystalline counterparts. Here, we conduct some classical molecular dynamics simulations of a-TiO2 based on a selection of empirical potentials. We show that, on account of its ability to dynamically assign the charge of each atom based on its local environment, the secondmoment tight-binding charge equilibration potential yields an unprecedented agreement with available experimental data. Based on these simulations, we investigate the degree of order and disorder in a-TiO2. Overall, the results suggest that a-TiO2 features a large flexibility in its local topology, which may explain the high sensitivity of its structure to the synthesis method being used. [ABSTRACT FROM AUTHOR]

Transport of saltwater through pristine and positively charged single-layer graphene nanoporous membranes is investigated using molecular dynamics simulations. Pressure-driven flows are induced by motion of specular reflecting boundaries at feed and permeate sides with constant speed. Unlike previous studies in the literature, this method induces a desired flow rate and calculates the resulting pressure difference in the reservoirs. Due to the hexagonal structure of graphene, the hydraulic diameters of nano-pores are used to correlate flow rate and pressure drop data. Simulations are performed for three different pore sizes and flow rates for the pristine and charged membrane cases. In order to create better statistical averages for salt rejection rates, ten different initial conditions of Na+ and Cl− distribution in the feed side are used for each simulation case. Using data from 180 distinct simulation cases and utilizing the Buckingham Pi theorem, we develop a functional relationship between the volumetric flow rate, pressure drop, pore diameter, and the dynamic viscosity of saltwater. A linear relationship between the volumetric flow rate and pressure drop is observed. For the same flow rate and pore size, charged membranes exhibit larger pressure drops. Graphene membranes with 9.90 Å pore diameter results in 100% salt rejection with 163.2 l/h cm2 water flux, requiring a pressure drop of 35.02 MPa. [ABSTRACT FROM AUTHOR]

Slow electron velocity-map imaging spectroscopy of cryogenically cooled TiO3H2− anions is used to probe the simplest titania/water reaction, TiO20/− + H2O. The resultant spectra show vibrationally resolved structure assigned to detachment from the cis-dihydroxide TiO(OH)2− geometry based on density functional theory calculations, demonstrating that for the reaction of the anionic TiO2− monomer with a single water molecule, the dissociative adduct (where the water is split) is energetically preferred over a molecularly adsorbed geometry. This work represents a significant improvement in resolution over previous measurements, yielding an electron affinity of 1.2529(4) eV as well as several vibrational frequencies for neutral TiO(OH)2. The energy resolution of the current results combined with photoelectron angular distributions reveals Herzberg-Teller coupling-induced transitions to Franck-Condon forbidden vibrational levels of the neutral ground state. A comparison to the previously measured spectrum of bare TiO2− indicates that reaction with water stabilizes neutral TiO2 more than the anion, providing insight into the fundamental chemical interactions between titania and water. [ABSTRACT FROM AUTHOR]

Simple and accurate expressions are presented for the equation of state (EOS) and absolute Helmholtz free energy of a system composed of simple atomic particles interacting through the repulsive Lennard- Jones potential model in the fluid and solid phases. The introduced EOS has 17 and 22 coefficients for fluid and solid phases, respectively, which are regressed to the Monte Carlo (MC) simulation data over the reduced temperature range of 0.6 ≤ T* ≤ 6.0 and the packing fraction range of 0.1 ≤ η ≤ 0.72. The average absolute relative percent deviation in fitting the EOS parameters to the MC data is 0.06 and 0.14 for the fluid and solid phases, respectively. The thermodynamic integration method is used to calculate the free energy using the MC simulation results. The Helmholtz free energy of the ideal gas is employed as the reference state for the fluid phase. For the solid phase, the values of the free energy at the reduced density equivalent to the close-packed of a hard sphere are used as the reference state. To check the validity of the predicted values of the Helmholtz free energy, the Widom particle insertion method and the Einstein crystal technique of Frenkel and Ladd are employed. The results obtained from the MC simulation approaches are well agreed to the EOS results, which show that the proposed model can reliably be utilized in the framework of thermodynamic theories. [ABSTRACT FROM AUTHOR]

We study effective two- and three-body interactions between non-active colloidal inclusions in an active bath of chiral or non-chiral particles, using Brownian dynamics simulations within a standard, two-dimensional model of disk-shaped inclusions and active particles. In a non-chiral active bath, we first corroborate previous findings on effective two-body repulsion mediated between the inclusions by elucidating the detailed non-monotonic features of the two-body force profiles, including a primary maximum and a secondary hump at larger separations that was not previously reported. We then show that these features arise directly from the formation, and sequential overlaps, of circular layers (or "rings") of active particles around the inclusions, as the latter are brought to small surface separations. These rings extend to radial distances of a few active-particle radii from the surface of inclusions, giving the hard-core inclusions relatively thick, soft, repulsive "shoulders," whose multiple overlaps then enable significant (non-pairwise) three-body forces in both non-chiral and chiral active baths. The resulting three-body forces can even exceed the two-body forces in magnitude and display distinct repulsive and attractive regimes at intermediate to large self-propulsion strengths. In a chiral active bath, we show that, while active particles still tend to accumulate at the immediate vicinity of the inclusions, they exhibit strong depletion from the intervening region between the inclusions and partial depletion from relatively thick, circular zones further away from the inclusions. In this case, the effective, predominantly repulsive interactions between the inclusions turn to active, chiralityinduced, depletion-type attractions, acting over an extended range of separations. [ABSTRACT FROM AUTHOR]

The magnitude of the pH of the surface of water continues to be a contentious topic in the physical chemistry of aqueous interfaces. Recent theoretical studies have shown little or no preference for the proton to be at the surface compared to the bulk. Using ab initio molecular dynamics simulations, we revisit the propensity of the excess proton for the air-water interface with a particular focus on the role of instantaneous liquid interfaces. We find a more pronounced presence for the proton to be at the air-water interface. The enhanced water structuring around the proton results in the presence of proton wires that run parallel to the surface as well as a hydrophobic environment made up of undercoordinated topological defect water molecules, both of which create favorable conditions for proton confinement at the surface. The Grotthuss mechanism within the structured water layer involves a mixture of both concerted and closely spaced stepwise proton hops. The proton makes excursions within the first solvation layer either in proximity to or along the instantaneous interface. [ABSTRACT FROM AUTHOR]

Crown-ether molecules are well known to selectively bind alkali atoms, so by incorporating these within wires, any change in electrical conductance of the wire upon binding leads to discriminating sensing. Using a density functional theory-based approach to quantum transport, we investigate the potential sensing capabilities of single-molecule junctions formed from crown ethers attached to anthraquinone units, which are in turn attached to gold electrodes via alkyl chains. We calculate the change in electrical conductance for binding of three different alkali ions (lithium, sodium, and potassium). Depending on the nature of the ionic analyte, the conductance is enhanced by different amounts. This change in electrical conductance is due to charge transfer from the ion to molecular wire causing the molecular resonances to shift closer to the electrode Fermi energy. [ABSTRACT FROM AUTHOR]

Ellipticine is a natural product that is currently being actively investigated for its inhibitory cancer and HIV properties. Here we use path-integral molecular dynamics coupled with excited state calculations to characterize the role of nuclear quantum effects on the structural and electronic properties of ellipticine in water, a common biological solvent. Quantum effects collectively enhance the fluctuations of both light and heavy nuclei of the covalent and hydrogen bonds in ellipticine. In particular, for the ellipticine-water system, where the proton donor and acceptor have different proton affinities, we find that nuclear quantum effects (NQEs) strengthen both the strong and the weak H bonds. This is in contrast to what is observed for the cases where the proton affinity of the donors and acceptors is same. These structural fluctuations cause a significant red-shift in the absorption spectra and an increase in the broadening, bringing it into closer agreement with the experiments. Our work shows that nuclear quantum effects alter both qualitatively and quantitatively the optical properties of this biologically relevant system and highlights the importance of the inclusion of these effects in the microscopic understanding of their optical properties. We propose that isotopic substitution will produce a blue shift and a reduction in the broadening of the absorption peak. [ABSTRACT FROM AUTHOR]

Triple differential cross sections (TDCSs) for the electron impact ionization of the unresolved combination of the 4 highest occupied molecular orbitals (4b3g, 5b2u, 1b1g, and 2b3u) of para-benzoquinone are reported. These were obtained in an asymmetric coplanar geometry with the scattered electron being observed at the angles -7.5°, -10.0°, -12.5° and -15.0°. The experimental cross sections are compared to theoretical calculations performed at the molecular 3-body distorted wave level, with a marginal level of agreement between them being found. The character of the ionized orbitals, through calculated momentum profiles, provides some qualitative interpretation for the measured angular distributions of the TDCS. [ABSTRACT FROM AUTHOR]

In the present work, we have measured the thermodiffusion coefficient of 51 binary liquid mixtures at 25 °C. These mixtures correspond to the series of the aromatics toluene and 1-methylnaphthalene with n-alkanes nCi (i = 6, 8, 10, 12, and 14) at different mass fractions in the whole range. For that, we have used the thermogravitational technique. It is shown that the thermodiffusion coefficient is a linear function of the mass fraction in all the mixtures. Extrapolating the lines, we obtain the thermodiffusion coefficient in dilute solutions of n-alkanes for both toluene and 1-methylnaphthalene. These limiting values show a linear dependence with the inverse of the product of the molecular weights. In addition, we have measured the molecular diffusion coefficient of all the mixtures at 0.5 of mass fraction and at 25 °C, by the sliding symmetric tubes technique. It is observed that the product of this coefficient with the viscosity at the same concentrations takes a constant value for each of the series considered. Finally, we have also determined the Soret coefficient of the equimass mixtures by the combination of the measurements of thermodiffusion and molecular diffusion coefficients. [ABSTRACT FROM AUTHOR]

We have performed 69Ga, 71 Ga, and 137 Ba NMR on Ba8Ga16 Ge30, a clathrate semiconductor which has been of considerable interest due to its large figure of merit for thermoelectric applications. In measurements from 4 K to 450 K, we used measurements on the two Ga nuclei to separate the magnetic and electric quadrupole hyperfine contributions and thereby gain information about the metallic and phonon behavior. The results show the presence of a pseudogap in the Ga electronic states within the conduction band, superposed upon a large Ba contribution to the conduction band. Meanwhile the phonon contributions to the Ga relaxation rates are large and increase more rapidly with temperature than typical semiconductors. These results provide evidence for enhanced anharmonicity of the propagative phonon modes over a wide range, providing experimental evidence for enhanced phonon-phonon scattering as a mechanism for the reduced thermal conductivity. [ABSTRACT FROM AUTHOR]

Starting from a microscopic model for a spherically symmetric active Janus particle, we study the interactions between two such active motors. The ambient fluid mediates a long range hydrodynamic interaction between two motors. This interaction has both direct and indirect hydrodynamic contributions. The direct contribution is due to the propagation of fluid flow that originated from a moving motor and affects the motion of the other motor. The indirect contribution emerges from the re-distribution of the ionic concentrations in the presence of both motors. Electric force exerted on the fluid from this ionic solution enhances the flow pattern and subsequently changes the motion of both motors. By formulating a perturbation method for very far separated motors, we derive analytic results for the translation and rotational dynamics of the motors. We show that the overall interaction at the leading order modifies the translational and rotational speeds of motors which scale as O([1/D]3) and O([1/D]4) with their separation, respectively. Our findings open up the way for studying the collective dynamics of synthetic micro-motors. [ABSTRACT FROM AUTHOR]

The development of autonomous micro/nanomotors driven by self-generated chemical gradients is a topic of high interest given their potential impact in medicine and environmental remediation. Although impressive functionalities of these devices have been demonstrated, a detailed understanding of the propulsion mechanism is still lacking. In this work, we perform a comprehensive numerical analysis of the key parameters governing the actuation of bimetallic catalytic micropumps. We show that the fluid motion is driven by self-generated electro-osmosis where the electric field originates by a proton current rather than by a lateral charge asymmetry inside the double layer. Hence, the surface potential and the electric field are the key parameters for setting the pumping strength and directionality. The proton flux that generates the electric field stems from the proton gradient induced by the electrochemical reactions taken place at the pump. Surprisingly the electric field and consequently the fluid flow are mainly controlled by the ionic strength and not by the conductivity of the solution, as one could have expected. We have also analyzed the influence of the chemical fuel concentration, electrochemical reaction rates, and size of the metallic structures for an optimized pump performance. Our findings cast light on the complex chemomechanical actuation of catalytic motors and provide important clues for the search, design, and optimization of novel catalytic actuators. [ABSTRACT FROM AUTHOR]

In this work, we use ab initio molecular dynamics coupled with metadynamics to explore and characterize the glassy potential energy landscape of the OH- in a 20 and 48 water cluster. The structural, energetic, and topological properties of OH- are characterized for both clusters and the molecular origins of the IR signatures are examined. We find that in both the small and large clusters, the OH- can donate or accept a varying number of hydrogen bonds confirming that the amphiphilic character does not depend on cluster size. However, we highlight some important differences found between the energetic and topological properties of both families of clusters which may have implications on understanding the changes in the solvation structure of OH- between bulk and interfacial environments. By studying the IR spectra of smaller subsets of molecules within the 20 water molecule cluster, we find that the IR spectrum of the bare OH- as well as the water molecule donating a strong hydrogen bond to it exhibits characteristic absorption along the amphiphilic band between 1500 and 3000 cm-1 at positions very similar to those found for the entire hydroxide cluster. The results presented here will be useful in the calibration and improvement of both ab initio and semi-empirical methods to model this complex anion. [ABSTRACT FROM AUTHOR]

We investigate the effective interaction between two randomly charged but otherwise net-neutral, planar dielectric slabs immersed in an asymmetric Coulomb fluid containing a mixture of mobile monovalent and multivalent ions. The presence of charge disorder on the apposed bounding surfaces of the slabs leads to substantial qualitative changes in the way they interact, as compared with the standard picture provided by the van der Waals and image-induced, ion-depletion interactions. While, the latter predict purely attractive interactions between strictly neutral slabs, we show that the combined effects from surface charge disorder, image depletion, Debye (or salt) screening, and also, in particular, their coupling with multivalent ions, give rise to a more diverse behavior for the effective interaction between net-neutral slabs at nano-scale separations. Disorder effects show large variation depending on the properly quantified strength of disorder, leading either to non-monotonic effective interaction with both repulsive and attractive branches when the surface charges are weakly disordered (small disorder variance) or to a dominating attractive interaction that is larger both in its range and magnitude than what is predicted from the van der Waals and image-induced, ion-depletion interactions, when the surfaces are strongly disordered (large disorder variance). [ABSTRACT FROM AUTHOR]

The He(I) photoelectron spectrum of furfural has been investigated, with its vibrational structure assigned for the first time. The ground and excited ionized states are assigned through ab initio calculations performed at the outer-valence Green's function level. Triple differential cross sections (TDCSs) for electron-impact ionization of the unresolved combination of the 4a"+ 21a' highest and next-highest occupied molecular orbitals have also been obtained. Experimental TDCSs are recorded in a combination of asymmetric coplanar and doubly symmetric coplanar kinematics. The experimental TDCSs are compared to theoretical calculations, obtained within a molecular 3-body distorted wave framework that employed either an orientation average or proper TDCS average. The proper average calculations suggest that they may resolve some of the discrepancies regarding the angular distributions of the TDCS, when compared to calculations employing the orbital average. [ABSTRACT FROM AUTHOR]

We propose a multireference linearized coupled cluster theory using matrix product states (MPSs- LCC) which provides remarkably accurate ground-state energies, at a computational cost that has the same scaling as multireference configuration interaction singles and doubles, for a wide variety of electronic Hamiltonians. These range from first-row dimers at equilibrium and stretched geometries to highly multireference systems such as the chromium dimer and lattice models such as periodic twodimensional 1-band and 3-band Hubbard models. The MPS-LCC theory shows a speed up of several orders of magnitude over the usual Density Matrix Renormalization Group (DMRG) algorithm while delivering energies in excellent agreement with converged DMRG calculations. Also, in all the benchmark calculations presented here, MPS-LCC outperformed the commonly used multi-reference quantum chemistry methods in some cases giving energies in excess of an order of magnitude more accurate. As a size-extensive method that can treat large active spaces, MPS-LCC opens up the use of multireference quantum chemical techniques in strongly correlated ab initio Hamiltonians, including twoand three-dimensional solids. [ABSTRACT FROM AUTHOR]

We numerically investigate the pressure drop reduction (PDR) performance of microchannels equipped with liquid-infused surfaces, along with determining the shape of the interface between the working fluid and lubricant within the microgrooves. The effects of different parameters, such as the Reynolds number of working fluid, density and viscosity ratios between the lubricant and working fluid, the ratio of the thickness of the lubricant layer over the ridges to the depth of the groove, and the Ohnesorge number as a representative of the interfacial tension, on the PDR and interfacial meniscus within the microgrooves are comprehensively studied. The results reveal that the density ratio and Ohnesorge number do not significantly affect the PDR. On the other hand, the viscosity ratio considerably affects the PDR, and a maximum PDR of 62% compared to a smooth non-lubricated microchannel is achieved for a viscosity ratio of 0.01. Interestingly, the higher the Reynolds number of the working fluid, the higher the PDR. The meniscus shape within the microgrooves is strongly affected by the Reynolds number of the working fluid. Despite the insignificant effect of interfacial tension on the PDR, the interface shape within the microgrooves is appreciably influenced by this parameter. [ABSTRACT FROM AUTHOR]

In this work, the transport coefficients of the ternary mixtures of the diffusion coefficient measurements in ternary mixtures 1 project were determined. The analyzed ternary mixtures are formed by 1,2,3,4-tetrahydronaphthalene, isobutylbenzene, and dodecane (nC12) at different compositions. In all cases, the analysis was carried out at 25 °C. The thermodiffusion coefficients were measured by a new thermogravitational column, and the molecular diffusion coefficients were determined by the sliding symmetric tubes technique. Finally, the Soret coefficients were ascertained from the measurements of the thermodiffusion and molecular diffusion coefficients. In addition, two new quantitative correlations which enable the prediction of the thermodiffusion and Soret coefficients of a ternary mixture are presented. The comparison between the experimental and the predicted data shows a good agreement. The presented results help to complete the lack of experimental data in ternary mixtures. In addition, this work improves the fundamental understanding of multicomponent mixtures. [ABSTRACT FROM AUTHOR]

Molecular dynamics simulations are carried out to study the translational and rotational diffusion of a single Janus particle immersed in a dense Lennard-Jones fluid. We consider a spherical particle with two hemispheres of different wettabilities. The analysis of the particle dynamics is based on the time-dependent orientation tensor, particle displacement, as well as the translational and angular velocity autocorrelation functions. It was found that both translational and rotational diffusion coefficients increase with decreasing surface energy at the nonwetting hemisphere, provided that the wettability of the other hemisphere remains unchanged. We also observed that in contrast to homogeneous particles, the nonwetting hemisphere of the Janus particle tends to rotate in the direction of the displacement vector during the rotational relaxation time. [ABSTRACT FROM AUTHOR]

We analyze van der Waals interactions between two rigid polymers with sequence-specific, anisotropic polarizabilities along the polymer backbones, so that the dipole moments fluctuate parallel to the polymer backbones. Assuming that each polymer has a quenched-in polarizability sequence which reflects, for example, the polynucleotide sequence of a double-stranded DNA molecule, we study the van derWaals interaction energy between a pair of such polymers with rod-like structure for the cases where their respective polarizability sequences are (i) distinct and (ii) identical, with both zero and non-zero correlation length of the polarizability correlator along the polymer backbones in the latter case. For identical polymers, we find a novel r-5 scaling behavior of the van der Waals interaction energy for small inter-polymer separation r, in contradistinction to the r-4 scaling behavior of distinct polymers, with furthermore a pronounced angular dependence favoring attraction between sufficiently aligned identical polymers. Such behavior can assist the molecular recognition between polymers. [ABSTRACT FROM AUTHOR]

The decomposition of electronic and nuclear motion presented in Abedi et al. [Phys. Rev. Lett. 105, 123002 (2010)] yields a time-dependent potential that drives the nuclear motion and fully accounts for the coupling to the electronic subsystem. Here, we show that propagation of an ensemble of independent classical nuclear trajectories on this exact potential yields dynamics that are essentially indistinguishable from the exact quantum dynamics for a model non-adiabatic charge transfer problem. We point out the importance of step and bump features in the exact potential that are critical in obtaining the correct splitting of the quasiclassical nuclear wave packet in space after it passes through an avoided crossing between two Born-Oppenheimer surfaces and analyze their structure. Finally, an analysis of the exact potentials in the context of trajectory surface hopping is presented, including preliminary investigations of velocity-adjustment and the force-induced decoherence effect. [ABSTRACT FROM AUTHOR]

Thermodynamic properties of the fluid in the hydrophobic pores of nanotubes are known to be different not only from the bulk phase but also from other conventional confinements. Here, we use a recently developed theoretical scheme of "two phase thermodynamic (2PT)" model to understand the driving forces inclined to spontaneous filling of carbon nanotubes (CNTs) with polar (water) and nonpolar (methane) fluids. The CNT confinement is found to be energetically favorable for both water and methane, leading to their spontaneous filling inside CNT(6,6). For both the systems, the free energy of transfer from bulk to CNT confinement is favored by the increased entropy (TAS), i.e., increased translational entropy and increased rotational entropy, which were found to be sufficiently high to conquer the unfavorable increase in enthalpy (AE) when they are transferred inside CNT. To the best of our knowledge, this is the first time when it has been established that the increase in translational entropy during confinement in CNT(6,6) is not unique to water-like H bonding fluid but is also observed in case of nonpolar fluids such as methane. The thermodynamic results are explained in terms of density, structural rigidity, and transport of fluid molecules inside CNT. The faster diffusion of methane over water in bulk phase is found to be reversed during the confinement in CNT(6,6). Studies reveal that though hydrogen bonding plays an important role in transport of water through CNT, but it is not the solitary driving factor, as the nonpolar fluids, which do not have any hydrogen bond formation capacity can go inside CNT and also can flow through it. The associated driving force for filling and transport of water and methane is enhanced translational and rotational entropies, which are attributed mainly by the strong correlation between confined fluid molecules and availability of more free space for rotation of molecule, i.e., lower density of fluid inside CNT due to their single file-like arrangement. To the best of our information, this is perhaps the first study of nonpolar fluid within CNT using 2PT method. Furthermore, the fast flow of polar fluid (water) over nonpolar fluid (methane) has been captured for the first time using molecular dynamic simulations. [ABSTRACT FROM AUTHOR]