Vertical melting of a stack of membranes

A stack of tensionless membranes with nonlinear curvature energy and vertical harmonic interaction is studied. At low temperatures, the system forms a lamellar phase. At a critical temperature, the stack disorders vertically in a melting-like transition. Received 26 April 2000 and Received in final form 16 August 2000     

A path integral approach to effective non-linear medium

In this article, we propose a new method to compute the effective properties of non-linear disordered media. We use the fact that the effective constants can be defined through the minimum of an energy functional. We express this minimum in terms of a path integral allowing us to use many-body techniques. We obtain the perturbation expansion of the effective constants to second order in disorder, for any kind of non-linearity. We apply our method to the case of strong non-linearities (i.e. , where is fluctuating from point to point), and to the case of weak non-linearity (i.e. where and fluctuate from point to point). Our results are in agreement with previous ones, and could be easily extended to other types of non-linear problems in disordered systems. Received: 13 May 1998 / Accepted: 27 July 1998     

A Raman study of polymerised C60

60 are reported. Although prepared according to different routes the Raman spectra of the two polymeric phases of C₆₀ show a quantitative agreement with respect to mode positions and intensity. We conclude from this that both materials have the same structure at least in the short range order, i.e. the same type of bonding and co-ordination between neighbouring C₆₀ molecules. An investigation of the time dependence of the thermal decomposition of high pressure polymerised C₆₀ is also presented. The rate of decomposition of the polymeric phase is found to be multi-exponential at all temperatures investigated. From an Arrhenius-type analysis of the short time data and the long time data, respectively, the activation energy for thermal dissociation of polymeric bonds was found to increase with time. Received: 20 September 1996/Accepted: 11 November 1996     

EPR investigation of a-Si:H aerosol particles formed under silane thermal decomposition

The dangling bond defects were investigated in a-Si:H particles formed under silane thermal decomposition in flow reactor. EPR together with hydrogen evolution method were used. The experimental results allowed us to conclude that there are two kinds of dangling bond defects in a-Si:H aerosol particles. The defects of the first kind are localized on the surface of interconnected microvoids and microchannels (surface dangling bonds) and those of the second kind are embedded in amorphous silicon network (volume dangling bonds). The thermal equilibration of dangling bonds and temperature dependence of equilibrium dangling bond concentration were investigated. It was found that at temperatures > 400 K the dangling bond concentrationNApplied Magnetic Resonance _s reversibly depends on sample temperature. The volume dangling bond concentration increases with temperature increasing (the effective activation energy of dangling bond formationU > 0), and the surface dangling bond concentration decreases with temperature increasing (U < 0). It has been found that EPR line is considerably asymmetric for samples with high hydrogen content and for low hydrogen content the EPR line is weakly asymmetric. A conclusion was drawn that the asymmetry degree depends on amorphous silicon lattice distortions. This conclusion has been confirmed by EPR spectra simulations.     

Energetics and electronic properties of BN nanocones with pentagonal rings at their apexes

The geometric structures, energetics and electronic properties of the recently discovered BN nanocones are investigated using first-principles calculations based on the density-functional theory. We have proposed one particular structure for BN nanocones associated with the 240^° disclination, derived by the extraction of four 60^° segments, presenting as characteristic four pentagons at the apex and termination in two atoms. The cones are simulated by three clusters containing 58 B plus N atoms and additional 12 H atoms to saturate the dangling bonds at the edge. The most stable configuration is obtained when the two terminating atoms are one B and one N. For the cases where the two terminating atoms are of the same kind, the tip with B atoms is determined to have lower binding energy than with N atoms. The local densities of states of these BN nanocones are investigated and sharp states are found in the regions close (below and above) to the Fermi energy. Received 14 October 2002 / Received in final form 6 December 2002 Published online 11 February 2003 RID="a" ID="a"e-mail: ppiquini@smail.ufsm.br     

A computer simulation study of the ground-state configurations of Fe and Fe-Al clusters

Using the noncentral embedded atom model potential recently proposed by Besson and Morillo for bulk alloys (), we performed computer simulations to predict the ground-state configurations of and clusters (). The computed structures of clusters are in general agreement with such theoretical results as have been obtained by density functional calculations (i.e. for ). The results for Fe-rich clusters show surface segregation of Al, which is in keeping with the findings of a previous study of clusters. Received 29 April 1999     

Connecting cluster dynamics and protein folding

The relaxation dynamics of clusters can be interpreted in terms of the topographies of their potential surfaces. Systems with short-range potentials have sawtooth-like potential surfaces with small drops in energy from one local minimum to the next and few-body motions as the clusters move from one minimum to another; such systems readily take on amorphous structures. These are called “glass-formers". Systems with long-range forces have potentials whose topographies are like rough staircases, with some large drops in energy from one minimum to the next; their well-to-well passages involve very collective motions and such systems are excellent structure-seekers. They find their way to well-ordered, highly selective structures under almost all circumstances. These characteristics generalize to describe the potential surfaces and folding behavior of polypeptides and proteins. The forces are effective long-range forces due to the polymer chain. Staircase topographies emerge both from direct sampling of potential surfaces and from the inversion of the kinetics generated by a much more aaabstract topological model, from which folding pathways can be inferred. Received 4 December 2000     

Assessing the validity of theoretical results

We present recent theoretical results for the V₃ and Au₄ clusters. Calculations of the V₃ doublet system indicate that the 6-311+G(d) basis set is sufficiently flexible to provide reliable minimum energy structures and vibrational frequencies, that these structures and frequencies are insensitive to spin contamination of the wave function when the BPW91 functional is used, and that changing to the B3LYP functional may result in very different structures and frequencies. A computationally less expensive scalar relativistic treatment of Au₄ clusters gives structural properties that are in good agreement with those obtained using a four-component method. Received 30 November 2000     

Martensitic transformations: first-principles calculations combined with molecular-dynamics simulations

Results are presented of first-principles total-energy calculations and molecular-dynamics simulations of structural transformations in magnetic transition metal alloys like Fe_1-xNi_x. While first-principles calculations allow to identify those structures having the lower total energy, molecular-dynamics simulations can be used to trace out the dependence of the transformation on temperature, composition, concentration of defects etc. We have used the method of the semi-empiric embedded-atom potential in the molecular-dynamics simulations which yields remarkable good results for the structural changes. Received: 3 February 1998 / Revised: 4 May 1998 / Accepted: 24 June 1998     

Ground-state geometries and stability of NaMg clusters using ab initio molecular dynamics method

The ground-state geometries, energetics and the stability of ( n =1-12) clusters are studied using ab initio molecular dynamics method. Our results indicate that the ground-state geometries of large clusters () are different from those of clusters where a trivalent impurity Al is added to the same monovalent host Na. Other features observed are an early appearance of 3-dimensional structure and a pentagonal growth path from n =6 up to n =11. As expected, the ground-state geometry of is not an icosahedron but can be viewed as a distorted form of one of the low lying geometries of cluster. In the energetically favored structures impurity atom Mg is never located at the center of the cluster. The stability analysis based on the energetics shows (8 valence electrons) to be the most stable. In addition there is a remarkable even-odd pattern observed in the dissociation energy and the second difference in energy which is absent in earlier studies of and clusters. Received: 16 September 1998 / Received in final form: 15 February 1999     

Continuum light emission from sputtered species of graphite during ion beam irradiation

Light emission during sputtering of graphite targets with 1–10 keV Ne⁺, Kr⁺ and Xe⁺ beams has been investigated in the 180-600 nm wavelength range. Beside the characteristic lines of sputtered C₁ and C₁ ⁺, a continuum superimposed with a number of broad structures was observed in the 250-520 nm range, and having a maximum at 386 nm. Mass analysis of the sputtered flux confirmed the presence of negative carbon clusters C _m ^- , C₂ ^- being the dominant one. Ion beam parameters i.e. ion mass, energy, current density and ion dose were varied to identify the origin of the continuum emission. On the basis of the experimental results, it is suggested that the continuum is predominantly due to the overlapping of various band systems of sputtered C₂ with a small contribution from the heavier sputtered carbon clusters C_m (m>2). Received 24 September 1999 and Received in final form 11 February 2000     

Vibrational spectroscopy of ice interfaces

Vibrational spectroscopy via infrared-visible sum-frequency generation was used to probe interfaces of hexagonal ice(0001) with vapor, a hydrophobic solid, and a hydrophilic solid. The OH stretch vibrational modes for hydrogen-bonded and free dangling OH bonds were measured at various temperatures below bulk melting. At the vapor/ice interface, the dangling OH mode provided information about surface melting. Surface structural disordering appeared to set in at ∼200 K and increased dramatically with temperature. Received: 15 January 2002 / Published online: 11 June 2002     

Overlap properties and adsorption transition of two Hamiltonian paths

We consider a model of two (fully) compact polymer chains, coupled through an attractive interaction. These compact chains are represented by Hamiltonian paths (HP), and the coupling favors the existence of common bonds between the chains. We use a (n=0 component) spin representation for these paths, and we evaluate the resulting partition function within a homogeneous saddle point approximation. For strong coupling (i.e. at low temperature), one finds a phase transition towards a “frozen” phase where one chain is completely adsorbed onto the other. By performing a Legendre transform, we obtain the probability distribution of overlaps. The fraction of common bonds between two HP, i.e. their overlap q, has both lower () and upper () bounds. This means in particular that two HP with overlap greater than coincide. These results may be of interest in (bio)polymers and in optimization problems. Received 4 December 1998 and Received in final form 10 March 1999     

Vitrification in a 2D Ising model with mobile bonds

A bond-disordered two-dimensional Ising model is used to simulate Kauzmann's mechanism of vitrification in liquids, by a Glauber Monte Carlo simulation. The rearrangement of configurations is achieved by allowing impurity bonds to hop to nearest neighbors at the same rate as the spins flip. For slow cooling, the theoretical minimum energy configuration is approached, characterized by an amorphous distribution of locally optimally arranged impurity bonds. Rapid cooling to low temperatures regularly finds bond configurations of higher energy, which are both a priori rare and severely restrictive to spin movement, providing a simple realization of kinetic vitrification. A supercooled liquid regime is also found, and characterized by a change in sign of the field derivative of the spin-glass susceptibility at a finite temperature. Received 3 August 2000 and Received in final form 9 March 2001     

Density functional calculations of dangling bonds for CN and SiN films

Hyperfine interaction constants (HFICs) of dangling bonds for CN and SiN films were calculated by density-functional theory. The averaged ¹⁴N isotropic HFICs for C dangling bonds are almost equal to those for Si dangling bonds. The anisotropic ¹⁴N HFICs calculated for C dangling bonds are larger than those for Si dangling bonds by a factor 2. The calculated results were compared with experimental results obtained by electron spin resonance and electron nuclear double resonance. It was indicated that the carbon dangling bonds are located such that they avoid N atoms in CN films.     

Calculation of the structure of carbon clusters based on fullerene-like C<Subscript>24</Subscript> and C<Subscript>48</Subscript> molecules

Equilibrium structures obtained by linking with valence bonds the carbon carcasses of two fullerene-like molecules have been studied by molecular dynamics simulation. In free fullerene, carbon atoms form sp² hybridized bonds, but at places of links between fullerenes, sp³ hybridized bonds are formed, which determines the changes in the properties of such structures. In the literature, the topology of diamond-like phases is described, but equilibrium clusters based on fullerene-like molecules are underexplored. The right angles between the C–C bonds are energetically unfavorable, and the reduction in the energy of clusters in the process of relaxation is connected with the optimization of valence angles, which leads to a reduction in the symmetry of clusters and, in a number of cases, even to disruption of some valence bonds. It is shown that different fashions of linking two fullerenes result in the formation of clusters with different structures and energies. Different initial conditions can lead to different configurations of clusters with the same topology. Among the analyzed clusters, a structure with the minimum potential energy per atom was found. The results of this work contribute to the study of the real structure of carbon clusters.     

Stability and fragmentation of multiply ionized clusters

The stability of neutral, singly and multiply ionized silicon clusters, (N = 2-7, M = 0, , , ), has been investigated using an ab initio density functional method. We show that the fragmentation effect significantly affects the structure of mass-spectra of multiply ionized silicon clusters. For clusters, the clusters with a large fragmentation energy are found to correspond to the high peaks at N = 4 and 6 in mass-spectra. For clusters, a peak at N = 5 in mass-spectra has been predicted to be especially high. Received: 9 June 1997 / Revised: 8 January 1998 / Accepted: 25 February 1998     

Energy characteristics of carbon clusters with passivated bonds

A modified phenomenological model is proposed for calculating the formation energy of carbon nanoclusters which makes it possible to analyze the regions of existence of clusters of various forms. A new parameter of the model, which corresponds to passivation of broken carbon bonds, affects the shape of the equilibrium optimum clusters, i.e., those having a minimum energy for a fixed number of atoms. Analytic dependences of equilibrium-configuration states determining the existence of spheroidal closed clusters, nanopipes, and fragments of a graphite plane, on the broken-bond energy parameter obtained in this model are presented. Fiz. Tverd. Tela (St. Petersburg) 41, 809–812 (May 1999)     

Broad histogram relation is exact

The Broad Histogram is a method designed to calculate the energy degeneracy g(E) from microcanonical averages of certain macroscopic quantities N _up and N _dn. These particular quantities are defined within the method, and their averages must be measured at constant energy values, i.e. within the microcanonical ensemble. Monte-Carlo simulational methods are used in order to perform these measurements. Here, the mathematical relation allowing one to determine g(E) from these averages is shown to be exact for any statistical model, i.e. any energy spectrum, under completely general conditions. We also comment about some troubles concerning the measurement of the quoted microcanonical averages, when one uses a particular approach, namely the energy random walk dynamics. These troubles appear when movements corresponding to different energy jumps are performed using the same probability, and also when the correlations between successive averaging states are not adequately treated: they have nothing to do with the method itself. Received: 29 May 1998 / Received in final form: 12 June 1998 / Accepted: 13 June 1998     

Melting behavior of large disordered sodium clusters

The melting-like transition in disordered sodium clusters Na₉₂ and Na₁₄₂ is studied by performing density functional constant-energy molecular dynamics simulations. The orbital-free version of the density functional formalism is used. In Na₁₄₂ the atoms are distributed in two distinct shells (surface and inner shells) and this cluster melts in two steps: the first one, at ≈130 K, is characterized by the development of a high intrashell atomic mobility, and the second, homogeneous melting at ≈270 K, involves diffusive motion of all the atoms across the whole cluster volume. On the contrary, the melting of Na₉₂ proceeds smoothly over a very broad temperature interval, without any abrupt change in the thermal or structural indicators. The occurrence of two steps in the melting transition is suggested to be related to the existence of a grouping of the atoms in radial shells, even if those shells present a large degree of internal disorder. It then appears that a cluster can be considered fully amorphous (totally disordered) only when there are no radial regions of low atomic density separating shells. The isomer of Na₉₂ studied here fulfills this criterion and its thermal behavior can be considered as representative of that expected for fully amorphous clusters. Disordered Na₁₄₂, on the other hand, that has a discernible structure of an inner and a surface shell, should be considered as not fully disordered. The thermal behavior of these two clusters is also compared to that of icosahedral (totally ordered) sodium clusters of the same sizes. Received 5 February 2001 and Received in final form 21 May 2001