Molecular theory of underdamped dielectric relaxation: understanding collective effects in dipolar liquids

A molecular theory of underdamped dielectric relaxation of a dense dipolar liquid is presented. This theory properly takes into account the collective effects that are present (due to strong intermolecular correlations) in a dipolar liquid. For small rigid molecules, the theory again leads to a three-variable description which, however, is somewhat different from the traditional version. In particular, two of the three parameters are collective in nature and are determined by the orientational pair correlation function. A detailed comparison between the theory and the computer simulation results of Neria and Nitzan is performed and an excellent agreement is obtained without the use of any adjustable or free parameter - the calculation is fully microscopic. The theory can also provide a systematic description of the Poley absorption often observed in dipolar liquids in the high-frequency regime.

Microscopic origin of the chirality driven morphologies of the amphiphilic monolayers and bilayers

It is widely known that the compressed monolayers and bilayers of chiral lipids or fatty acids form helical morphologies, while the corresponding racemic modification gives only flat platelets without twist. No molecular explanation of this phenomenon is yet available, although subtle interactions at the chiral centers have often been proposed as the driving force behind the morphology of the aggregate to form a particular shape. In the present study, the morphologies of the chiral amphiphilic assemblies have been predicted on the basis of an effective pair potential between the molecules, which depends on the relative sizes of the groups attached to the chiral centers, the orientation of the amphiphilic molecules and also on the distance between them. It is shown that fur a pair of same kind of enantiomers, the minimum energy conformation favours a twist angle between them. This twist between the neighbouring molecules gives rise to the helicity of the aggregate. The present theory also shows from the molecular considerations that for a pair of mirror-image isomers (i.e. the racemic modification) the minimum energy conformation corresponds to the zero angle between the molecules, thus giving rise to flat platelets as observed in experiments. Another fascinating aspect of such chirality driven helical structures is that the sense (or the handedness) of the helix is highly specific about the chirality of the monomer concerned. The molecular theory shows, for the first time, that the sense of the helical structures in many cases is determined by the sizes of the groups attached to the chiral centers and the effective potential between them. The predicted senses of the helical structures are in complete agreement with the experimental results.

Isomerization dynamics in viscous liquids: Microscopic investigation of the coupling and decoupling of the rate to and from solvent viscosity and dependence on the intermolecular potential

A detailed investigation of viscosity dependence of the isomerization rate is carried out for continuous potentials by using a fully microscopic, self-consistent mode-coupling theory calculation of both the friction on the reactant and the viscosity of the medium. In this calculation we avoid approximating the short time response by the Enskog limit, which overestimates the friction at high frequencies. The isomerization rate is obtained by using the Grote-Hynes formula. The viscosity dependence of the rate has been investigated for a large number of thermodynamic state points. Since the activated barrier crossing dynamics probes the high-frequency frictional response of the liquid, the barrier crossing rate is found to be sensitive to the nature of the reactant-solvent interaction potential. When the solute-solvent interaction is modeled by a 6-12 Lennard-Jones potential, we find that over a large variation of viscosity (eta), the rate (k) can indeed be fitted very well to a fractional viscosity dependence: (k similar to eta(-alpha)), with the exponent alpha in the range 1 greater than or equal to alpha >0. The calculated values of the exponent appear to be in very good agreement with many experimental results. In particular, the theory, for the first time, explains the experimentally observed high value of alpha even at the barrier frequency, omega(b). similar or equal to 9 X 10(12) s(-1) for the isomerization reaction of 2-(2'-propenyl)anthracene in liquid eta-alkanes. The present study can also explain the reason for the very low value of vb observed in another study for the isomerization reaction of trans-stilbene in liquid n-alkanes. For omega(b) greater than or equal to 2.0 X 10(13) s(-1), we obtain alpha similar or equal to 0, which implies that the barrier crossing rate becomes identical to the transition-state theory predictions. A careful analysis of isomerization reaction dynamics involving large amplitude motion suggests that the barrier crossing dynamics itself may become irrelevant in highly viscous liquids and the rate might again be coupled directly to the viscosity. This crossover is predicted to be strongly temperature dependent and could be studied by changing the solvent viscosity by the application of pressure. (C) 1999 American Institute of Physics. [S0021-9606(9950514-X].

Non-Fermi-liquid theory for disordered metals near two dimensions

We consider the Finkelstein action describing a system of spin-polarized or spinless electrons in 2+2epsilon dimensions, in the presence of disorder as well as the Coulomb interactions. We extend the renormalization-group analysis of our previous work and evaluate the metal-insulator transition of the electron gas to second order in an epsilon expansion. We obtain the complete scaling behavior of physical observables like the conductivity and the specific heat with varying frequency, temperature, and/or electron density. We extend the results for the interacting electron gas in 2+2epsilon dimensions to include the quantum critical behavior of the plateau transitions in the quantum Hall regime. Although these transitions have a very different microscopic origin and are controlled by a topological term in the action (theta term), the quantum critical behavior is in many ways the same in both cases. We show that the two independent critical exponents of the quantum Hall plateau transitions, previously denoted as nu and p, control not only the scaling behavior of the conductances sigma(xx) and sigma(xy) at finite temperatures T, but also the non-Fermi-liquid behavior of the specific heat (c(v)proportional toT(p)). To extract the numerical values of nu and p it is necessary to extend the experiments on transport to include the specific heat of the electron gas.

Anisotropic isometric fluctuation relations in experiment and theory on a self-propelled rod

The isometric fluctuation relation (IFR) P. I. Hurtado et al., Proc. Natl. Acad. Sci. USA 108, 7704 (2011)] relates the relative probability of current fluctuations of fixed magnitude in different spatial directions. We test its validity in an experiment on a tapered rod, rendered motile by vertical vibration and immersed in a sea of spherical beads. We analyze the statistics of the velocity vector of the rod and show that they depart significantly from the IFR of Hurtado et al. Aided by a Langevin-equation model we show that our measurements are largely described by an anisotropic generalization of the IFR R. Villavicencio et al., Europhys. Lett. 105, 30009 (2014)], with no fitting parameters, but with a discrepancy in the prefactor whose origin may lie in the detailed statistics of the microscopic noise. The experimentally determined large-deviation function of the velocity vector has a kink on a curve in the plane.

Mode coupling theory analysis of electrolyte solutions: Time dependent diffusion, intermediate scattering function, and ion solvation dynamics

A self-consistent mode coupling theory (MCT) with microscopic inputs of equilibrium pair correlation functions is developed to analyze electrolyte dynamics. We apply the theory to calculate concentration dependence of (i) time dependent ion diffusion, (ii) intermediate scattering function of the constituent ions, and (iii) ion solvation dynamics in electrolyte solution. Brownian dynamics with implicit water molecules and molecular dynamics method with explicit water are used to check the theoretical predictions. The time dependence of ionic self-diffusion coefficient and the corresponding intermediate scattering function evaluated from our MCT approach show quantitative agreement with early experimental and present Brownian dynamic simulation results. With increasing concentration, the dispersion of electrolyte friction is found to occur at increasingly higher frequency, due to the faster relaxation of the ion atmosphere. The wave number dependence of intermediate scattering function, F(k, t), exhibits markedly different relaxation dynamics at different length scales. At small wave numbers, we find the emergence of a step-like relaxation, indicating the presence of both fast and slow time scales in the system. Such behavior allows an intriguing analogy with temperature dependent relaxation dynamics of supercooled liquids. We find that solvation dynamics of a tagged ion exhibits a power law decay at long times-the decay can also be fitted to a stretched exponential form. The emergence of the power law in solvation dynamics has been tested by carrying out long Brownian dynamics simulations with varying ionic concentrations. The solvation time correlation and ion-ion intermediate scattering function indeed exhibit highly interesting, non-trivial dynamical behavior at intermediate to longer times that require further experimental and theoretical studies. (c) 2015 AIP Publishing LLC.

The kinetic-theory for dilute solid liquid 2-phase flow

On the basis of a brief review of the continuum theory for macroscopic descriptions and the kinetic theory for microscopic descriptions in solid/liquid two-phase flows, some suggestions are presented, i.e. the solid phase may be described by the Boltzmann equation and the liquid phase still be described by conservation laws in the continuum theory. Among them the action force on the particles by the liquid fluid is a coupling factor which connects the phases. For dilute steady solid/liquid two-phase flows, the particle velocity distribution function can be derived by analogy with the procedures in the kinetic theory of gas molecules for the equilibrium state instead of being assumed, as previous investigators did. This done, more detailed information, such as the velocity probability density distribution, mean velocity distribution and fluctuating intensity etc. can be obtained directly from the particle velocity distribution function or from its integration. Experiments have been performed for dilute solid/liquid two-phase flow in a 4 x 6 cm2 sized circulating square pipe system by means of laser Doppler anemometry so that the theories can be examined. The comparisons show that the theories agree very well with all the measured data.

Macroscopically Dissipative Systems with Underlying Microscopic Dynamics : Properties and Limits of Measurement

While some of the deepest results in nature are those that give explicit bounds between important physical quantities, some of the most intriguing and celebrated of such bounds come from fields where there is still a great deal of disagreement and confusion regarding even the most fundamental aspects of the theories. For example, in quantum mechanics, there is still no complete consensus as to whether the limitations associated with Heisenberg's Uncertainty Principle derive from an inherent randomness in physics, or rather from limitations in the measurement process itself, resulting from phenomena like back action. Likewise, the second law of thermodynamics makes a statement regarding the increase in entropy of closed systems, yet the theory itself has neither a universally-accepted definition of equilibrium, nor an adequate explanation of how a system with underlying microscopically Hamiltonian dynamics (reversible) settles into a fixed distribution.

Motivated by these physical theories, and perhaps their inconsistencies, in this thesis we use dynamical systems theory to investigate how the very simplest of systems, even with no physical constraints, are characterized by bounds that give limits to the ability to make measurements on them. Using an existing interpretation, we start by examining how dissipative systems can be viewed as high-dimensional lossless systems, and how taking this view necessarily implies the existence of a noise process that results from the uncertainty in the initial system state. This fluctuation-dissipation result plays a central role in a measurement model that we examine, in particular describing how noise is inevitably injected into a system during a measurement, noise that can be viewed as originating either from the randomness of the many degrees of freedom of the measurement device, or of the environment. This noise constitutes one component of measurement back action, and ultimately imposes limits on measurement uncertainty. Depending on the assumptions we make about active devices, and their limitations, this back action can be offset to varying degrees via control. It turns out that using active devices to reduce measurement back action leads to estimation problems that have non-zero uncertainty lower bounds, the most interesting of which arise when the observed system is lossless. One such lower bound, a main contribution of this work, can be viewed as a classical version of a Heisenberg uncertainty relation between the system's position and momentum. We finally also revisit the murky question of how macroscopic dissipation appears from lossless dynamics, and propose alternative approaches for framing the question using existing systematic methods of model reduction.

A molecular dynamics study of the microscopic properties of simple dense fluids

The microscopic properties of a two-dimensional model dense fluid of Lennard-Jones disks have been studied using the so-called "molecular dynamics" method. Analyses of the computer-generated simulation data in terms of "conventional" thermodynamic and distribution functions verify the physical validity of the model and the simulation technique.

The radial distribution functions g(r) computed from the simulation data exhibit several subsidiary features rather similar to those appearing in some of the g(r) functions obtained by X-ray and thermal neutron diffraction measurements on real simple liquids. In the case of the model fluid, these "anomalous" features are thought to reflect the existence of two or more alternative configurations for local ordering.

Graphical display techniques have been used extensively to provide some intuitive insight into the various microscopic phenomena occurring in the model. For example, "snapshots" of the instantaneous system configurations for different times show that the "excess" area allotted to the fluid is collected into relatively large, irregular, and surprisingly persistent "holes". Plots of the particle trajectories over intervals of 2.0 to 6.0 x 10^-12 sec indicate that the mechanism for diffusion in the dense model fluid is "cooperative" in nature, and that extensive diffusive migration is generally restricted to groups of particles in the vicinity of a hole.

A quantitative analysis of diffusion in the model fluid shows that the cooperative mechanism is not inconsistent with the statistical predictions of existing theories of singlet, or self-diffusion in liquids. The relative diffusion of proximate particles is, however, found to be retarded by short-range dynamic correlations associated with the cooperative mechanism--a result of some importance from the standpoint of bimolecular reaction kinetics in solution.

A new, semi-empirical treatment for relative diffusion in liquids is developed, and is shown to reproduce the relative diffusion phenomena observed in the model fluid quite accurately. When incorporated into the standard Smoluchowski theory of diffusion-controlled reaction kinetics, the more exact treatment of relative diffusion is found to lower the predicted rate of reaction appreciably.

Finally, an entirely new approach to an understanding of the liquid state is suggested. Our experience in dealing with the simulation data--and especially, graphical displays of the simulation data--has led us to conclude that many of the more frustrating scientific problems involving the liquid state would be simplified considerably, were it possible to describe the microscopic structures characteristic of liquids in a concise and precise manner. To this end, we propose that the development of a formal language of partially-ordered structures be investigated.

Theory and simulation in neuroscience.

Modeling work in neuroscience can be classified using two different criteria. The first one is the complexity of the model, ranging from simplified conceptual models that are amenable to mathematical analysis to detailed models that require simulations in order to understand their properties. The second criterion is that of direction of workflow, which can be from microscopic to macroscopic scales (bottom-up) or from behavioral target functions to properties of components (top-down). We review the interaction of theory and simulation using examples of top-down and bottom-up studies and point to some current developments in the fields of computational and theoretical neuroscience.

Superfluidity in neutron matter and symmetric nuclear matter and the effect of a microscopic three-body force

The neutron (PF2)-P-3 pairing gap in pure neutron matter, neutron (PF2)-P-3 gap and neutron-proton (SD1)-S-3 gap in symmetric nuclear matter have been studied by using the Brueckner-Hartree-Fock(BHF) approach and the BCS theory. We have concentrated on investigating and discussing the three-body force effect on the nucleon superfluidity. The calculated results indicate that the three-body force enhances remaxkably the (PF2)-P-3 superfluidity in neutron matter. It also enhances the (PF2)-P-3 superfluidity in symmetric nuclear matter and its effect increases monotonically as the Fermi-momentum k(F) increases, whereas the three-body force is shown to influence only weakly the neutron-proton (SD1)-S-3 gap in symmetric nuclear matter.

Neutron (PF2)-P-3 superfluidity in neutron matter and the effect of microscopic three-body force

The neutron (PF2)-P-3 pairing gap in pure neutron matter has been studied by using the Brueckner-Hartree-Fock( BHF) approach and the BCS theory. We have concentrated our attention on investigating the three-body force effect on the neutron superfluidity in the (PF2)-P-3 channel. The calculated results indicate that the three-body force enhances remarkably the (PF2)-P-3 superfluidity in neutron matter. When adopting the BHF single-particle spectrum, the three-body force turns out to increase the maximum value of the pairing gap from about 0.22 MeV to about 0.5 MeV.

Microscopic examination of NpNn dominance in the evolution of nuclear structure

The proton-neutron interaction in determining the evolution of nuclear structure has been studied by using the Brillouin-Wigner perturbation expansion. The particle-hole and particle-particle p-n interactions are unifiedly described in the theory. The obtained formulas of level energies and excitation energies scaled in the small- and large-NpNn limits can well explain the linearity of the extracted proton-neutron interaction energies and the attenuation of the 2(1)(+) excitation energies against the valence nucleon product NpNn for five mass regions from A = 100-200.

A microscopic study of structural and electronic properties of functionalized silicon surfaces based on first-principles

Surface modification of silicon with organic monolayers tethered to the surface by different linkers is an important process in realizing future (opto-)electronic devices. Understanding the role played by the nature of the linking group and the chain length on the adsorption structures and electronic properties of these assemblies is vital to advance this technology. This Thesis is a study of such properties and contributes in particular to a microscopic understanding of induced changes in the work function of experimentally studied functionalized silicon surfaces. Using first-principles density functional theory (DFT), at the first step, we provide predictions for chemical trends in the work function of hydrogenated silicon (111) surfaces modified with various terminations. For nonpolar terminating atomic species such as F, Cl, Br, and I, the change in the work function is directly proportional to the amount of charge transferred from the surface, thus relating to the difference in electronegativity of the adsorbate and silicon atoms. The change is a monotonic function of coverage in this case, and the work function increases with increasing electronegativity. Polar species such as −TeH, −SeH, −SH, −OH, −NH2, −CH3, and −BH2 do not follow this trend due to the interaction of their dipole with the induced electric field at the surface. In this case, the magnitude and sign of the surface dipole moment need to be considered in addition to the bond dipole to generally describe the change in work function. Compared to hydrogenated surfaces, there is slight increase in the work function of H:Si(111)-XH, where X = Te, Se, and S, whereas reduction is observed for surfaces covered with −OH, −CH3, and −NH2. Next, we study the hydrogen passivated Si(111) surface modified with alkyl chains of the general formula H:Si–(CH2)n–CH2 and H:Si–X–(CH2)n–CH3, where X = NH, O, S and n = (0, 1, 3, 5, 7, 9, 11), at half coverage. For (X)–Hexyl and (X)–Dodecyl functionalization, we also examined various coverages up to full monolayer grafting in order to validate the result of half covered surface and the linker effect on the coverage. We find that it is necessary to take into account the van der Waals interaction between the alkyl chains. The strongest binding is for the oxygen linker, followed by S, N, and C, irrespective of chain length. The result revealed that the sequence of the stability is independent of coverage; however, linkers other than carbon can shift the optimum coverage considerably and allow further packing density. For all linkers apart from sulfur, structural properties, in particular, surface-linker-chain angles, saturate to a single value once n > 3. For sulfur, we identify three regimes, namely, n = 0–3, n = 5–7, and n = 9–11, each with its own characteristic adsorption structures. Where possible, our computational results are shown to be consistent with the available experimental data and show how the fundamental structural properties of modified Si surfaces can be controlled by the choice of linking group and chain length. Later we continue by examining the work function tuning of H:Si(111) over a range of 1.73 eV through adsorption of alkyl monolayers with general formula -[Xhead-group]-(CnH2n)-[Xtail-group], X = O(H), S(H), NH(2). The work function is practically converged at 4 carbons (8 for oxygen), for head-group functionalization. For tail-group functionalization and with both head- and tail-groups, there is an odd-even effect in the behavior of the work function, with peak-to-peak amplitudes of up to 1.7 eV in the oscillations. This behavior is explained through the orientation of the terminal-group's dipole. The shift in the work function is largest for NH2-linked and smallest for SH-linked chains and is rationalized in terms of interface dipoles. Our study reveals that the choice of the head- and/or tail-groups effectively changes the impact of the alkyl chain length on the work function tuning using self-assembled monolayers and this is an important advance in utilizing hybrid functionalized Si surfaces. Bringing together the understanding gained from studying single type functionalization of H:Si(111) with different alkyl chains and bearing in mind how to utilize head-group, tail-group or both as well as monolayer coverage, in the final part of this Thesis we study functionalized H:Si(111) with binary SAMs. Aiming at enhancing work function adjustment together with SAM stability and coverage we choose a range of terminations and linker-chains denoted as –X–(Alkyl) with X = CH3, O(H), S(H), NH(2) and investigate the stability and work function of various binary components grafted onto H:Si(111) surface. Using binary functionalization with -[NH(2)/O(H)/S(H)]-[Hexyl/Dodecyl] we show that work function can be tuned within the interval of 3.65-4.94 eV and furthermore, enhance the SAM’s stability. Although direct Si-C grafted SAMs are less favourable compared to their counterparts with O, N or S linkage, regardless of the ratio, binary functionalized alkyl monolayers with X-alkyl (X = NH, O) is always more stable than single type alkyl functionalization with the same coverage. Our results indicate that it is possible to go beyond the optimum coverage of pure alkyl functionalized SAMs (50%) by adding a linker with the correct choice of the linker. This is very important since dense packed monolayers have fewer defects and deliver higher efficiency. Our results indicate that binary anchoring can modify the charge injection and therefore bond stability while preserving the interface electronic structure.

Catalytic role of gold in gold-based catalysts: A density functional theory study on the CO oxidation on gold

Gold-based catalysts have been of intense interests in recent years, being regarded as a new generation of catalysts due to their unusually high catalytic performance. For example, CO oxidation on Au/TiO2 has been found to occur at a temperature as low as 200 K. Despite extensive studies in the field, the microscopic mechanism of CO oxidation on Au-based catalysts remains controversial. Aiming to provide insight into the catalytic roles of Au, we have performed extensive density functional theory calculations for the elementary steps in CO oxidation on Au surfaces. O atom adsorption, CO adsorption, O-2 dissociation, and CO oxidation on a series of Au surfaces, including flat surfaces, defects and small clusters, have been investigated in detail. Many transition states involved are located, and the lowest energy pathways are determined. We find the following: (i) the most stable site for O atom on Au is the bridge site of step edge, not a kink site; (ii) O-2 dissociation on Au (O-2-->20(ad)) is hindered by high barriers with the lowest barrier being 0.93 eV on a step edge; (iii) CO can react with atomic O with a substantially lower barrier, 0.25 eV, on Au steps where CO can adsorb; (iv) CO can react with molecular O-2 on Au steps with a low barrier of 0.46 eV, which features an unsymmetrical four-center intermediate state (O-O-CO); and (v) O-2 can adsorb on the interface of Au/TiO2 with a reasonable chemisorption energy. On the basis of our calculations, we suggest that (i) O-2 dissociation on Au surfaces including particles cannot occur at low temperatures; (ii) CO oxidation on Au/inactive-materials occurs on Au steps via a two-step mechanism: CO+O-2-->CO2+O, and CO+O-->CO2; and (iii) CO oxidation on Au/active-materials also follows the two-step mechanism with reactions occurring at the interface.