Search

General features of the $\alpha-\beta$ transition of quartz are investigated. Molecular dynamics methods are mainly used, an analytic treatment being deferred to a work in preparation. A basic preliminary observation is that the transition involves only a subsystem of four normal modes on which the remaining ones just act as a reservoir. The dynamics of the relevant subsystem turns out to be Hamiltonian, being governed by an effective potential that depends on the specific energy of the total system. The effective potential is actually calculated through time averages. It describes the transition as a pitchfork bifurcation, and also explains the phenomenon of the soft mode, since it exhibits a frequency that vanishes at the transition. The critical exponent too is estimated., Comment: 14 pages, 4 figures

We report the results of a numerical investigation, performed in the frame of dynamical systems' theory, for a realistic model of a ionic crystal for which, due to the presence of long--range Coulomb interactions, the Gibbs distribution is not well defined. Taking initial data with a Maxwell-Boltzmann distribution for the mode-energies $E_k$, we study the dynamical fluctuations, computing the moduli of the the energy-changes $|E_k(t)-E_k(0)|$. The main result is that they follow Tsallis distributions, which relax to distributions close to Maxwell-Boltzmann ones; indications are also given that the system remains correlated. The relaxation time $\tau$ depends on specific energy $\varepsilon$, and for the curve $\tau$ vs, $\varepsilon$ one has two results. First, there exists an energy threshold $\varepsilon_0$, above which the curve has the form $$ \tau \cdot \varepsilon \simeq h\ , $$ where, unexpectedly, Planck's constant $h$ shows up. In terms of the standard deviation $\Delta E$ of a mode-energy (for which one has $\Delta E=\varepsilon$), denoting by $\Delta t$ the relaxation time $\tau$, the relation reads $\Delta t \cdot \Delta E \simeq h$, which reminds of the Heisenberg uncertainty relation. Moreover, the threshold corresponds to zero-point energy. Indeed, the quantum value of the latter is $h\nu/2$ ( where $\nu$ is the characterisic infrared frequency of the system), while we find $\varepsilon \simeq h\nu/4$, so that one only has a discrepancy of a factor 2. So it seems that lack of full chaoticity manifests itself, in Statistical Thermodynamics, through quantum-like phenomena.

According to the correspondence principle, classical mechanics and quantum mechanics agree in the semiclassical limit, although presently it has become more and more clear how intriguing would be to try to fix a boundary between them. Here we give a significant example in which the agreement concerns Newtonian trajectories of an electron with initial data corresponding to a quantum ground state. The example is the simplest case in which a chemical bond occurs, i.e. the $H_2^+$ ion. By molecular dynamics simulations for the full system (two protons and one electron) we show that there exist initial data producing an ``effective potential'' among the protons, which superposes in a surprisingly good way the quantum one computed in the Born-Oppenheimer approximation (Fig~1). Preliminarily, following the perturbation procedure first exhibited by Born and Heisenberg in the year 1924, we recall why an effective potential should exist in a classical frame, and also describe the numerical procedure employed in computing it.

In Statistical Mechanics, Tsallis distributions were apparently conceived in connection with systems presenting long--range interactions. In fact, they were observed in numerical computations for models of such a type, as occurring in the approach to equilibrium, i.e., to a Maxwell--Boltzmann distribution. Here we exhibit two apparently new results. The first one is that Tsallis distributions occur also in an ionic--crystal model with long--range Coulomb forces, which is so realistic as to reproduce in an impressively good way the experimental infrared spectra. Thus such distributions may be expected to be actual physical features of crystals. The second result is that Tsallis distributions occur in the standard short--range FPU model too, so that the presence of long--range interactions is not a necessary condition for Tsallis distributions to occur. In fact, this is in agreement with a previous result of the first author in connection with the statistics of return times for the classical FPU model. We thus confirm the thesis advanced by Tsallis himself, that the relevant property for a dynamical system to present Tsallis distributions is that its dynamics should be not fully chaotic, a property which is known to actually pertain to long--range systems., Comment: Proceedings of the conference: XVI course of the Erice School on Complexity

It is well known that Gibbs' statistical mechanics is not justified for systems presenting long-range interactions, such as plasmas or galaxies. In a previous work we considered a realistic FPU-like model of an ionic crystal (and thus with long-range interactions), and showed that it reproduces the experimental infrared spectra from 1000 K down to 7 K, provided one abandons the Gibbs identification of temperature in terms of specific kinetic energy, at low temperatures. Here we investigate such a model in connection with its ergodicity properties. The conclusion we reach is that at low temperatures ergodicity does not occur, and thus the Gibbs prescriptions are not dynamically justified, up to geological time scales. We finally give a preliminary result indicating how the so-called `nonclassical' q-statistics show up in the realistic ionic-crystal model. How to formulate a consistent statistical mechanics, with the corresponding suitable identification of temperature in such nonergodicity conditions, remains an open problem, which apparently constitutes the modern form of the FPU problem.

It was recently shown that the experimental infrared spectra of ionic crystals at room temperature are very well reproduced by classical realistic models, and here new results are reported on the temperature dependence of the spectra, for the LiF crystal. The principal aim of the present work is however to highlight the deep analogy existing between the problem of spectra in ionic crystal models on the one hand, and that of energy equipartition in the Fermi--Pasta--Ulam model, on the other. Indeed at low temperatures the computations of the spectra show that the dynamics of the considered system is not completely chaotic, so that the use of the Boltzmann--Gibbs statistics is put in question, as in the Fermi--Pasta--Ulam case. Here, however, at variance with the equipartition problem, a first positive indication is given on the modifications that should be introduced in a classical statistical treatment: the new results at low temperatures show that it is indeed the Clausius identification of temperature that has to be modified. In fact, at very low temperatures the theoretical spectra fail to reproduce the experimental ones, if the temperature is taken as proportional to mean kinetic energy, but agreement is recovered through the only expedient of introducing a suitable temperature rescaling. Analogous results are also found in connection with thermal expansion., Comment: Sent for publication in Physica A

This paper is concerned with theoretical estimates of the refractive--index curves for quartz, obtained by the Kubo formul\ae\ in the classical approximation, through MD simulations for the motions of the ions. Two objectives are considered. The first one is to understand the role of nonlinearities in situations where they are very large, as at the $\alpha$--$\beta$ structural phase transition. We show that on the one hand they don't play an essential role in connection with the form of the spectra in the infrared. On the other hand they play an essential role in introducing a chaoticity which involves a definite normal mode. This might explain why that mode is Raman active in the $\alpha$ phase, but not in the $\beta$ phase. The second objective concerns whether it is possible in a microscopic model to obtain normal mode frequencies, or peak frequencies in the optical spectra, that are in good agreement with the experimental data for quartz. Notwithstanding a lot of effort, we were unable to find results agreeing better than about 6%, as apparently also occurs in the whole available literature. We interpret this fact as indicating that some essential qualitative feature is lacking in all models which consider, as the present one, only short--range repulsive potentials and unretarded long--range electric forces., Comment: 16 pages, 3 tables and 8 figures

The theoretical dispersion curves $n(\omega)$ (refractive index versus frequency) of ionic crystals in the infrared domain are expressed, within the Green--Kubo theory, in terms of a time correlation function involving the motion of the ions only. The aim of this paper is to investigate how well the experimental data are reproduced by a classical approximation of the theory, in which the time correlation functions are expressed in terms of the ions orbits. We report the results of molecular dynamics (MD) simulations for the ions motions of a LiF lattice of 4096 ions at room temperature. The theoretical curves thus obtained are in surprisingly good agreement with the experimental data, essentially over the whole infrared domain. This shows that at room temperature the motion of the ions develops essentially in a classical regime.

The dimension-four genuine anomalous quartic couplings are studied in processes of six-fermion production via e+e- collisions. Complete tree-level electroweak calculations are performed including initial-state-radiation and beamstrahlung. The analysis of final-state distributions can be used to find kinematical cuts to enhance the effects of anomalous couplings. For the parameters of the custodial-symmetry-conserving anomalous couplings a sensitivity in the range between 10^{-3} and 10^{-2} can be expected at 1 TeV., Comment: LaTeX (using elsart.sty). 15 pages, 6 figures. Submitted to Nuclear Physics B