Static structure and dynamics of the liquid Li-Na and Li-Mg alloys

We present calculations for the static structure and ordering properties of two lithium-based s-p bonded liquid alloys, Li-Na and Li-Mg. Our theoretical approach is based on the neutral pseudoatom method to derive the interatomic pair potentials, and on the modified-hypernetted-chain theory of liquids to obtain the liquid static structure, leading to a whole combination that is free of adjustable parameters. The study is complemented by performing molecular dynamics simulations which, besides checking the theoretical static structural results, also allow a calculation of some dynamical properties. The obtained results are compared with the available experimental data.

Self-dynamic structure factor of dense liquids: Theory and simulation

The self-intermediate dynamic structure factor Fs(k,t) of liquid lithium near the melting temperature is calculated by molecular dynamics. The results are compared with the predictions of several theoretical approaches, paying special attention to the Lovesey model and the Wahnstrm and Sjgren mode-coupling theory. To this end the results for the Fs(k,t) second memory function predicted by both models are compared with the ones calculated from the simulations.

The Impact of CAP Reforms on Farm Labour Structure. Factor Markets Working Document No. 63, August 2013

The labour force engaged in the agricultural sector is declining over time, and one can observe the reallocation of labour from family members to hired workers. Using farm-level data, this paper analyses the on-farm labour structure in Greece and assesses the factors driving its evolution over the period 1990-2008. The impact of agricultural policies and farm characteristics is examined in a dynamic panel analysis. Family and hired labour are found to be substitutes rather than complements, while agricultural support measures appear to negatively affect demand for both family and hired labour. Decoupled payments and subsidies on crops are found to have a significant impact on both sources of labour, as well as subsidies for rural development that do not favour on-farm labour use. The paper also finds that structural labour adjustments are the result of farm characteristics, such as farm size and location. The results are robust to various estimation techniques and specifications.

Molecular dynamics simulation of liquid trimethylphosphine

Structural and dynamical properties of liquid trimethylphosphine (TMP), (CH(3))(3)P, as a function of temperature is investigated by molecular dynamics (MD) simulations. The force field used in the MD simulations, which has been proposed from molecular mechanics and quantum chemistry calculations, is able to reproduce the experimental density of liquid TMP at room temperature. Equilibrium structure is investigated by the usual radial distribution function, g(r), and also in the reciprocal space by the static structure factor, S(k). On the basis of center of mass distances, liquid TMP behaves like a simple liquid of almost spherical particles, but orientational correlation due to dipole-dipole interactions is revealed at short-range distances. Single particle and collective dynamics are investigated by several time correlation functions. At high temperatures, diffusion and reorientation occur at the same time range as relaxation of the liquid structure. Decoupling of these dynamic properties starts below ca. 220 K, when rattling dynamics of a given TMP molecules due to the cage effect of neighbouring molecules becomes important. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3624408]

Einfluß der internen Architektur von Polymermikronetzwerken auf Struktur und Dynamik konzentrierter Kolloid-Dispersionen

Einfluß der internen Architektur von Polymermikronetzwerken auf Struktur und Dynamik konzentrierter Kolloid-Dispersionen Kugelförmige Polymermikronetzwerk-Kolloide gehören zur Klasse der sogenannten Mikrogele. Dabei handelt es sich um kolloidale Modellsysteme, die durch ihre interne Vernetzungsdichte charakterisiert werden.In dieser Arbeit sollte untersucht werden, ob sich die Wechselwirkungen zwischen den Mikrogel-Kolloiden über ein repulsives Potential der Form U(r) = 1/rn beschreiben lassen und ob der Poten-tialexponent n von der Vernetzungsdichte abhängt. Dazu wurden vor allem die innere Architektur, das Phasenverhalten und der statische Strukturfaktor 1:10, 1:50, 1:72 und 1:100 vernetzter Polymer-Mikronetzwerk-Kolloide bis in den Bereich hochkonzentrierter Dispersionen mit den Mitteln der Kleinwinkelneutronenstreuung, der Digitalphotographie und der statischen Lichtstreuung untersucht. Polymeranalytische Untersuchungen ergaben einen bei der Synthese anfallenden Anteil von unver-netztem, freiem Polymer innerhalb der Mikronetzwerke, welcher sich beim Lösen aus den Netzwerken herausbewegte. Das freie Polymer spielte vor allem beim Phasenverhalten der untersuchten Teilchen eine große Rolle und verursachte bei den Untersuchungen der statischen Strukturfaktoren Abweichun-gen vom 'harte Kugel'-Verhalten. Als Ergebnis der Kleinwinkel-Neutronenstreuung konnte eine ab-nehmende Verteilungsdichte der Vernetzer innerhalb der Polymermikronetzwerke in Richtung der -Teilchenoberfläche nachgewiesen werden. Die damit verbundene Konformationsfreiheit der Polymer-segmente auf der Teilchenoberfläche (bis hin zu 'mushroom'-Strukturen) wurde als Grund dafür an-gesehen, daß sich die Resultate der untersuchten Mikrogele aller Vernetzungsdichten im wesentlichen auf 'harte Kugeln' skalieren lassen.

Intrinsic compressibility and volume compression in solvated proteins by molecular dynamics simulation at high pressure.

Constant pressure and temperature molecular dynamics techniques have been employed to investigate the changes in structure and volumes of two globular proteins, superoxide dismutase and lysozyme, under pressure. Compression (the relative changes in the proteins' volumes), computed with the Voronoi technique, is closely related with the so-called protein intrinsic compressibility, estimated by sound velocity measurements. In particular, compression computed with Voronoi volumes predicts, in agreement with experimental estimates, a negative bound water contribution to the apparent protein compression. While the use of van der Waals and molecular volumes underestimates the intrinsic compressibilities of proteins, Voronoi volumes produce results closer to experimental estimates. Remarkably, for two globular proteins of very different secondary structures, we compute identical (within statistical error) protein intrinsic compressions, as predicted by recent experimental studies. Changes in the protein interatomic distances under compression are also investigated. It is found that, on average, short distances compress less than longer ones. This nonuniform contraction underlines the peculiar nature of the structural changes due to pressure in contrast with temperature effects, which instead produce spatially uniform changes in proteins. The structural effects observed in the simulations at high pressure can explain protein compressibility measurements carried out by fluorimetric and hole burning techniques. Finally, the calculation of the proteins static structure factor shows significant shifts in the peaks at short wavenumber as pressure changes. These effects might provide an alternative way to obtain information concerning compressibilities of selected protein regions.

Static and dynamic structure of liquid metals: Role of the different parts of the interaction potential

The influence of different parts of the interaction potential on the microscopic behavior of simple liquid metals is investigated by molecular dynamics simulation. The role of the soft-core repulsive, short-range attractive, and long-range oscillatory forces on the properties of liquid lithium close to the triple point is analyzed by comparing the results from simulations of identical systems but truncating the potential at different distances. Special attention is paid to dynamic collective properties such as the dynamic structure factors, transverse current correlation functions, and transport coefficients. It is observed that, in general, the effects of the short-range attractive forces are important. On the contrary, the influence of the oscillatory long-range interactions is considerably less, being the most pronounced for the dynamic structure factor at long wavelengths. The results of this work suggest that the influence of the attractive forces becomes less significant when temperature and density increase.

Electronic structure and properties of Cu2O

The structural and electronic properties of Cu2O have been investigated using the periodic Hartree-Fock method and a posteriori density-functional corrections. The lattice parameter, bulk modulus, and elastic constants have been calculated. The electronic structure of and bonding in Cu2O are analyzed and compared with x-ray photoelectron spectroscopy spectra, showing a good agreement for the valence-band states. To check the quality of the calculated electron density, static structure factors and Compton profiles have been calculated, showing a good agreement with the available experimental data. The effective electron and hole masses have been evaluated for Cu2O at the center of the Brillouin zone. The calculated interaction energy between the two interpenetrated frameworks in the cuprite structure is estimated to be around -6.0 kcal/mol per Cu2O formula. The bonding between the two independent frameworks has been analyzed using a bimolecular model and the results indicate an important role of d10-d10 type interactions between copper atoms.

Electronic structure and properties of Cu2O

The structural and electronic properties of Cu2O have been investigated using the periodic Hartree-Fock method and a posteriori density-functional corrections. The lattice parameter, bulk modulus, and elastic constants have been calculated. The electronic structure of and bonding in Cu2O are analyzed and compared with x-ray photoelectron spectroscopy spectra, showing a good agreement for the valence-band states. To check the quality of the calculated electron density, static structure factors and Compton profiles have been calculated, showing a good agreement with the available experimental data. The effective electron and hole masses have been evaluated for Cu2O at the center of the Brillouin zone. The calculated interaction energy between the two interpenetrated frameworks in the cuprite structure is estimated to be around -6.0 kcal/mol per Cu2O formula. The bonding between the two independent frameworks has been analyzed using a bimolecular model and the results indicate an important role of d10-d10 type interactions between copper atoms.

Electronic structure and properties of Cu2O

The structural and electronic properties of Cu2O have been investigated using the periodic Hartree-Fock method and a posteriori density-functional corrections. The lattice parameter, bulk modulus, and elastic constants have been calculated. The electronic structure of and bonding in Cu2O are analyzed and compared with x-ray photoelectron spectroscopy spectra, showing a good agreement for the valence-band states. To check the quality of the calculated electron density, static structure factors and Compton profiles have been calculated, showing a good agreement with the available experimental data. The effective electron and hole masses have been evaluated for Cu2O at the center of the Brillouin zone. The calculated interaction energy between the two interpenetrated frameworks in the cuprite structure is estimated to be around -6.0 kcal/mol per Cu2O formula. The bonding between the two independent frameworks has been analyzed using a bimolecular model and the results indicate an important role of d10-d10 type interactions between copper atoms.

Computer simulations of the structure of colloidal ferrofluids

The structure of a ferrofluid under the influence of an external magnetic field is expected to become anisotropic due to the alignment of the dipoles into the direction of the external field, and subsequently to the formation of particle chains due to the attractive head to tail orientations of the ferrofluid particles. Knowledge about the structure of a colloidal ferrofluid can be inferred from scattering data via the measurement of structure factors. We have used molecular-dynamics simulations to investigate the structure of both monodispersed and polydispersed ferrofluids. The results for the isotropic structure factor for monodispersed samples are similar to previous data by Camp and Patey that were obtained using an alternative Monte Carlo simulation technique, but in a different parameter region. Here we look in addition at bidispersed samples and compute the anisotropic structure factor by projecting the q vector onto the XY and XZ planes separately, when the magnetic field was applied along the z axis. We observe that the XY- plane structure factor as well as the pair distribution functions are quite different from those obtained for the XZ plane. Further, the two- dimensional structure factor patterns are investigated for both monodispersed and bidispersed samples under different conditions. In addition, we look at the scaling exponents of structure factors. Our results should be of value to interpret scattering data on ferrofluids obtained under the influence of an external field.

Amorphous ZrF4-a molecular dynamics study

Using an effective two-body interaction potential, a molecular dynamics study of the structural properties of amorphous ZrF4 phase is presented. The effective pair potential includes steric repulsion, Coulomb interaction due to charge transfer, and charge-dipole interaction due to the large electronic polarizability of anions. The results for structural correlations, such as pair distribution functions, coordination numbers, and bond angle distributions are presented. Excellent agreement is obtained by comparing experimental X-ray diffraction and the simulated static X-ray structure factor. © 1993.

The Equilibrium Structure of Lithium Salt Solutions in Ether-Functionalized Ammonium Ionic Liquids

Molecular dynamics simulations have been performed for ionic liquids based on a ternary mixture of lithium and ammonium cations and a common anion, bis(trifluoromethylsulfonyl)imide, [Tf2N](-). We address structural changes resulting from adding Li+ in ionic liquids with increasing length of an ether-functionalized chain in the ammonium cation. The calculation of static structure factors reveals the lithium effect on charge ordering and intermediate range order in comparison with the neat ionic liquids. The charge ordering is modified in the lithium solution because the coordination of [Tf2N](-) toward Li+ is much stronger than ammonium cations. Intermediate range order is observed in neat ionic liquids based on ammonium cations with a long chain, but in the lithium solutions, there is also a nonhomogenous distribution of Li+ cations. The presence of Li+ enhances interactions between the ammonium cations due to correlations between the oxygen atom of the ether chain and the nitrogen atom of another ammonium cation.

Dynamik von Polymerschmelzen und Quellverhalten ungeordneter Netzwerke

Deutsch:Mit Hilfe eines parallelen Molekulardynamik-Programmswurden einfach Modelle von Homopolymerschmelzen simuliert.Langkettige Schmelzen zeigten eine sehr gute Übereinstimmungmit den Vorhersagendes Reptationsmodells. Die intermediären Reptationsbereichemit den vorhergesagtenExponenten konnten wesentlich klarer als bisher verifiziertwerden. Es stellteVerschiedene, gebräuchliche Analyse-Methoden führten jedochzuunterschiedlichen Aussagen für die Verhakungslänge. Fürkurze Kettenbzw. kurze Unterketten wurden in dichten Schmelzendie Abweichungen vom Rouse-Modell aufgezeigt. DieseAbweichungenkönnen als Korrelationslocheffekt interpretiert werden undsteht in teilweiser Übereinstimmung zu den Vorhersagenrenormierter Rouse-Modelle. Aus den Schmelzen wurden Netzwerke in einem speziellenZufallsvernetzungsprozeßhergestellt, der die Bildung von Defektstrukturenunterbindet. Diesewurden bzgl. ihres Quell- und Deformationsverhaltenuntersucht. Der maximaleQuellgrad eines Netzwerkes war bereits bei verhaltnismäßigkurzen Kettenlängenverhakungslimitiert. Die Struktur der Ketten in einem biszum osmotischenGleichgewicht gequollenen Netzwerk unterhalb derMaschengröße ist die überstreckter,selbstvermeidender Ketten mit einer Fraktaldimension von D =1.4,jenseits der Maschengröße (bzw. Verhakungslänge) nehmen siedie Struktur einesIrrfluges an (D = 2). Gequollene Netzwerke zeigten, wie auch in Experimenten, einestarkeZunahme von Dichtefluktuationen, welche unter Verstreckunganisotropwurde und in der Streufunktion zu sogenanntenButterfly-Mustern führt.Diese Fluktuationen sind statischer Natur als Folge desEinfrierenseines ungeordneten Zustandes während der Vernetzung.

Dynamics, rheology and critical properties of colloidal fluid mixtures

Liquids under the influence of external fields exhibit a wide range of intriguing phenomena that can be markedly different from the behaviour of a quiescent system. This work considers two different systems — a glassforming Yukawa system and a colloid-polymer mixture — by Molecular Dynamics (MD) computer simulations coupled to dissipative particle dynamics. The former consists of a 50-50 binary mixture of differently-sized, like-charged colloids interacting via a screened Coulomb (Yukawa) potential. Near the glass transition the influence of an external shear field is studied. In particular, the transition from elastic response to plastic flow is of interest. At first, this model is characterised in equilibrium. Upon decreasing temperature it exhibits the typical dynamics of glassforming liquids, i.e. the structural relaxation time τα grows strongly in a rather small temperature range. This is discussed with respect to the mode-coupling theory of the glass transition (MCT). For the simulation of bulk systems under shear, Lees-Edwards boundary conditions are applied. At constant shear rates γ˙ ≫ 1/τα the relevant time scale is given by 1/γ˙ and the system shows shear thinning behaviour. In order to understand the pronounced differences between a quiescent system and a system under shear, the response to a suddenly commencing or terminating shear flow is studied. After the switch-on of the shear field the shear stress shows an overshoot, marking the transition from elastic to plastic deformation, which is connected to a super-diffusive increase of the mean squared displacement. Since the average static structure only depends on the value of the shear stress, it does not discriminate between those two regimes. The distribution of local stresses, in contrast, becomes broader as soon as the system starts flowing. After a switch-off of the shear field, these additional fluctuations are responsible for the fast decay of stresses, which occurs on a time scale 1/γ˙ . The stress decay after a switch-off in the elastic regime, on the other hand, happens on the much larger time scale of structural relaxation τα. While stresses decrease to zero after a switch-off for temperatures above the glass transition, they decay to a finite value for lower temperatures. The obtained results are important for advancing new theoretical approaches in the framework of mode-coupling theory. Furthermore, they suggest new experimental investigations on colloidal systems. The colloid-polymer mixture is studied in the context of the behaviour near the critical point of phase separation. For the MD simulations a new effective model with soft interaction potentials is introduced and its phase diagram is presented. Here, mainly the equilibrium properties of this model are characterised. While the self-diffusion constants of colloids and polymers do not change strongly when the critical point is approached, critical slowing down of interdiffusion is observed. The order parameter fluctuations can be determined through the long-wavelength limit of static structure factors. For this strongly asymmetric mixture it is shown how the relevant structure factor can be extracted by a diagonalisation of a matrix that contains the partial static structure factors. By presenting first results of this model under shear it is demonstrated that it is suitable for non-equilibrium simulations as well.