Molecular dynamics simulations of the dipolar-induced formation of magnetic nanochains and nanorings

Iron, cobalt and nickel nanoparticles, grown in the gas phase, are known to arrange in chains and bracelet-like rings due to the long-range dipolar interaction between the ferromagnetic (or super-paramagnetic) particles. We investigate the dynamics and thermodynamics of such magnetic dipolar nanoparticles for low densities using molecular dynamics simulations and analyze the influence of temperature and external magnetic fields on two- and three-dimensional systems. The obtained phase diagrams can be understood by using simple energetic arguments.Comment: 6 pages, 6 figure

Free Energy Approach to the Formation of an Icosahedral Structure during the Freezing of Gold Nanoclusters

The freezing of metal nanoclusters such as gold, silver, and copper exhibits a novel structural evolution. The formation of the icosahedral (Ih) structure is dominant despite its energetic metastability. This important phenomenon, hitherto not understood, is studied by calculating free energies of gold nanoclusters. The structural transition barriers have been determined by using the umbrella sampling technique combined with molecular dynamics simulations. Our calculations show that the formation of Ih gold nanoclusters is attributed to the lower free energy barrier from the liquid to the Ih phases compared to the barrier from the liquid to the face-centered-cubic crystal phases

Particle-size dependence of orbital order-disorder transition in LaMnO3

The latent heat (L) associated with the orbital order-disorder transition at T_JT is found to depend significantly on the average particle size (d) of LaMnO3. It rises slowly with the decrease in d down to ~100 nm and then jumps by more than an order of magnitude in between d ~ 100 nm and ~30 nm. Finally, L falls sharply to zero at a critical particle size d_c ~ 19 nm. The transition temperature T_JT also exhibits an almost similar trend of variation with the particle size, near d ~ 30 nm and below, even though the extent of variation is relatively small. The zero-field-cooled (ZFC) and field-cooled (FC) magnetization versus temperature study over a temperature range 10-300 K reveals that the antiferromagnetic transition temperature decreases with d while the temperature range, over which the ZFC and FC data diverge, increases with the drop in d. The FC magnetization also is found to increase sharply with the drop in particle size. A conjecture of nonmonotonic variation in orbital domain structure with decrease in particle size - from smaller domains with large number of boundaries to larger domains with small number of boundaries due to lesser lattice defects and, finally, down to even finer domain structures with higher degree of metastability - along with increase in surface area in core-shell structure, could possibly rationalize the observed L versus d and T_JT versus d patterns. Transmission electron microscopy data provide evidence for presence of core-shell structure as well as for increase in lattice defects in finer particles.Comment: 26 pages including 5 figures; pdf only; accepted for publication in Phys. Rev.

Premelting of Thin Wires

Recent work has raised considerable interest on the nature of thin metallic wires. We have investigated the melting behavior of thin cylindrical Pb wires with the axis along a (110) direction, using molecular dynamics and a well-tested many-body potential. We find that---in analogy with cluster melting---the melting temperature $T_m (R)$ of a wire with radius $R$ is lower than that of a bulk solid, $T_m^b$, by $T_m (R) = T_m^b -c/R$. Surface melting effects, with formation of a thin skin of highly diffusive atoms at the wire surface, is observed. The diffusivity is lower where the wire surface has a flat, local (111) orientation, and higher at (110) and (100) rounded areas. The possible relevance to recent results on non-rupturing thin necks between an STM tip and a warm surface is addressed.Comment: 10 pages, 4 postscript figures are appended, RevTeX, SISSA Ref. 131/94/CM/S

Why do gallium clusters have a higher melting point than the bulk?

Density functional molecular dynamical simulations have been performed on Ga$_{17}$ and Ga$_{13}$ clusters to understand the recently observed higher-than-bulk melting temperatures in small gallium clusters [Breaux {\em et al.}, Phys. Rev. Lett. {\bf 91}, 215508 (2003)]. The specific-heat curve, calculated with the multiple-histogram technique, shows the melting temperature to be well above the bulk melting point of 303 K, viz. around 650 K and 1400 K for Ga$_{17}$ and Ga$_{13}$, respectively. The higher-than-bulk melting temperatures are attributed mainly to the covalent bonding in these clusters, in contrast with the covalent-metallic bonding in the bulk.Comment: 4 pages, including 6 figures. accepted for publication in Phys. Rev. Let

Finite size melting of spherical solid-liquid aluminium interfaces

We have investigated the melting of nano-sized cone shaped aluminium needles coated with amorphous carbon using transmission electron microscopy. The interface between solid and liquid aluminium was found to have spherical topology. For needles with fixed apex angle, the depressed melting temperature of this spherical interface, with radius $R$, was found to scale linearly with the inverse radius $1/R$. However, by varying the apex angle of the needles we show that the proportionality constant between the depressed melting temperature and the inverse radius changes significantly. This lead us to the conclusion that the depressed melting temperature is not controlled solely by the inverse radius $1/R$. Instead we found a direct relation between the depressed melting temperature and the ratio between the solid-liquid interface area and the molten volume.Comment: to appear in Philosophical Magazine (2009

Melting behavior of ultrathin titanium nanowires

The thermal stability and melting behavior of ultrathin titanium nanowires with multi-shell cylindrical structures are studied using molecular dynamic simulation. The melting temperatures of titanium nanowires show remarkable dependence on wire sizes and structures. For the nanowire thinner than 1.2 nm, there is no clear characteristic of first-order phase transition during the melting, implying a coexistence of solid and liquid phases due to finite size effect. An interesting structural transformation from helical multi-shell cylindrical to bulk-like rectangular is observed in the melting process of a thicker hexagonal nanowire with 1.7 nm diameter.Comment: 4 pages, 4 figure

Thermodynamics of tin clusters

We report the results of detailed thermodynamic investigations of the Sn$_{20}$ cluster using density-functional molecular dynamics. These simulations have been performed over a temperature range of 150 to 3000 K, with a total simulation time of order 1 ns. The prolate ground state and low-lying isomers consist of two tricapped trigonal prism (TTP) units stacked end to end. The ionic specific heat, calculated via a multihistogram fit, shows a small peak around 500 K and a shoulder around 850 K. The main peak occurs around 1200 K, about 700 K higher than the bulk melting temperature, but significantly lower than that for Sn$_{10}$. The main peak is accompanied by a sharp change in the prolate shape of the cluster due to the fusion of the two TTP units to form a compact, near spherical structure with a diffusive liquidlike ionic motion. The small peak at 500 K is associated with rearrangement processes within the TTP units, while the shoulder at 850 K corresponds to distortion of at least one TTP unit, preserving the overall prolate shape of the cluster. At all temperatures observed, the bonding remains covalent.Comment: Latex File and EPS Figures. 18 pages,11 Figures. Submitted to Phys. Rev.

Simulation of the thermally induced austenitic phase transition in NiTi nanoparticles

The reverse martensitic ("austenitic") transformation upon heating of equiatomic nickel-titanium nanoparticles with diameters between 4 and 17 nm is analyzed by means of molecular-dynamics simulations with a semi-empirical model potential. After constructing an appropriate order parameter to distinguish locally between the monoclinic B19' at low and the cubic B2 structure at high temperatures, the process of the phase transition is visualized. This shows a heterogeneous nucleation of austenite at the surface of the particles, which propagates to the interior by plane sliding, explaining a difference in austenite start and end temperatures. Their absolute values and dependence on particle diameter are obtained and related to calculations of the surface induced size dependence of the difference in free energy between austenite and martensite.Comment: 6 pages, 4 figures, accepted for publication in "The European Physical Journal B