Magnetic ordering in strong anisotropic ferromagnetic bilayers with dipolar and antiferromagnetic interlayer exchange interactions: Monte Carlo approach.

In this work, multilayer films consisting of two strong-anisotropic ferromagnetic layers antiferromagnetically coupled by a nonmagnetic spacer are studied by Monte Carlo simulations. The system is modeled by an Ising-based Hamiltonian that depends on both the intralayer exchange and dipolar constants and on the interlayer exchange constant (IEC). The ground state of the monolayers (null IEC) corresponds to alternate stripe domains with width h defined by the ratio between the exchange and dipolar constants (δ). The results show that IEC alters the energy balance that controls the stripe domain formation, leading to a ground state characterized by in-plane stripes out-plane antiferromagnetically coupled. When temperature increases two regimes are identified: an IEC-dominated regime where the orientational and positional orders are simultaneously lost in both layers, driving the system to the tetragonal liquid (TL) phase, and a dipolar-dominated one where signs of layers decoupling and the onset of positional disorder are observed. The last could be related with an intermediate nematic phase (NM). From the study of the nonequilibrium dynamics, the phase transitions to TL phase are characterized as continuous and those to the NM one as Kosterlitz-Thouless type. Also, for both layers the critical temperatures are the same and increase with IEC magnitude. Furthermore, the obtained critical exponents depend on the IEC values, which is indicative of a weak universality. For the dipolar-dominated regime, the decoupling between layers is also evidenced by the difference between their critical exponents.

Evidence of Kosterlitz-Thouless phase transitions in the Ising model with dipolar interactions.

The ferromagnetic bidimensional Ising model with dipolar interactions has been proposed to model ultrathin films with strong out-of-plane anisotropy. The phase diagram presents a rich phenomenology that includes low-temperature phases characterized by stripes of width n (h_{n}) and a high-temperature phase with domains of stripes with mutually perpendicular orientations, named tetragonal liquid (TL). The latter phase can be reached by two possible ways. One of them is the direct transition h_{n} to TL, and the other one is through an intermediate phase with orientational order but short-range positional disorder, named nematic phase (NM). The regions of the phase diagram where these transitions occur, as well as their character, remain an open question and are the object of the present work. In order to clarify this topic, intensive Monte Carlo simulations were performed by employing short-time dynamics as the main tool for studying the phase transition behavior. The dynamic evolution of the orientational order parameter and its moments are measured for selected values of the ratio between the ferromagnetic exchange and dipolar constants, called δ. The obtained results indicate that the intermediate NM phase is present for δ≥2 in narrow ranges of temperatures. Also, the results suggest that both transitions, i.e., h_{n}-NM and NM-TL, have a Kosterlitz-Thouless character. This type of topological transition is observed in continuous bidimensional models and have been proposed for discrete ones, as in the case of the present work.

Nanoclusters of crystallographically aligned nanoparticles for magnetic thermotherapy: aqueous ferrofluid, agarose phantoms and ex vivo melanoma tumour assessment.

Magnetic hyperthermia is an oncological therapy where magnetic nanostructures, under a radiofrequency field, act as heat transducers increasing tumour temperature and killing cancerous cells. Nanostructure heating efficiency depends both on the field conditions and on the nanostructure properties and mobility inside the tumour. Such nanostructures are often incorrectly bench-marketed in the colloidal state and using field settings far off from the recommended therapeutic values. Here, we prepared nanoclusters composed of iron oxide magnetite nanoparticles crystallographically aligned and their specific absorption rate (SAR) values were calorimetrically determined in physiological fluids, agarose-gel-phantoms and ex vivo tumours extracted from mice challenged with B16-F0 melanoma cells. A portable, multipurpose applicator using medical field settings; 100 kHz and 9.3 kA m-1, was developed and the results were fully analysed in terms of nanoclusters' structural and magnetic properties. A careful evaluation of the nanoclusters' heating capacity in the three milieus clearly indicates that the SAR values of fluid suspensions or agarose-gel-phantoms are not adequate to predict the real tissue temperature increase or the dosage needed to heat a tumour. Our results show that besides nanostructure mobility, perfusion and local thermoregulation, the nanostructure distribution inside the tumour plays a key role in effective heating. A suppression of the magnetic material effective heating efficiency appears in tumour tissue. In fact, dosage had to be increased considerably, from the SAR values predicted from fluid or agarose, to achieve the desired temperature increase. These results represent an important contribution towards the design of more efficient nanostructures and towards the clinical translation of hyperthermia.

Phase transitions and multicritical behavior in the Ising model with dipolar interactions.

In this work, the phase diagram of the ferromagnetic Ising model with dipole interactions is revisited with the aim of determining the nature of the phase transition between stripe-ordered phases with width n (h_{n}) and tetragonal liquid (TL) phases. Extensive Monte Carlo simulations are performed in order to study the short-time dynamic behavior of the observables for selected values of the ratio between the ferromagnetic exchange and dipolar constants δ. The obtained results indicate that the h_{1}-TL phase transition line is continuous up to δ=1.2585, while for the h_{2}-TL line a weak first-order character is found in the interval 1.2585≤δ≤1.36 and becomes continuous for 1.37≤δ≤1.9. This result suggests the existence of a tricritical point close to δ=1.37. When it is appropriate, the complete set of critical exponents is obtained, and in all the studied cases they depend on δ but do not belong to the Ising universality class. Furthermore, short-time dynamic studies reveal that at the point where the mentioned lines meet the h_{1}-h_{2} line, i.e., at δ=1.2585, the critical phase corresponding to the h_{1}-TL transition coexists with the h_{2} phase.

Phase transitions and critical phenomena in the two-dimensional Ising model with dipole interactions: A short-time dynamics study.

The ferromagnetic Ising model with antiferromagnetic dipole interactions is investigated by means of Monte Carlo simulations, focusing on the characterization of the phase transitions between the tetragonal liquid and stripe of width h phases. The dynamic evolution of the physical observables is analyzed within the short-time regime for 0.5≤δ≤1.3, where δ is the ratio between the short-range exchange and the long-range dipole interaction constants. The obtained results for the interval 0.5≤δ≤1.2 indicate that the phase transition line between the h=1 stripe and tetragonal liquid phases is continuous. This finding contributes to clarifying the controversy about the order of this transition. This controversy arises from the difficulties introduced in the simulations due to the presence of long-range dipole interactions, such as an important increase in the simulation times that limits the system size used, strong finite size effects, as well as to the existence of multiple metastable states at low temperatures. The study of the short-time dynamics of the model allows us to avoid these hindrances. Moreover, due to the fact that the finite-size effects do not significantly affect the power-law behavior exhibited in the observables within the short-time regime, the results could be attributed to those corresponding to the thermodynamic limit. As a consequence of this, a careful characterization of the critical behavior for the whole transition line is performed by giving the complete set of critical exponents.

Shearing self-propelled suspensions: arrest of coarsening and suppression of giant density fluctuations.

We study the effect of a linear shear flow on a collection of interacting active, self-propelled particles modeled via the Vicsek model. The imposed flow has a dramatic effect on the behavior of the model. We find that in the presence of shear there is no order-disorder transition, and that coarsening of the domains is arrested. Shear also suppresses the so-called giant density fluctuations that are observed in the quiescent limit.

Monte Carlo study of the phase transition in the critical behavior of the Ising model with shear.

The critical behavior of the Ising model with nonconserved dynamics and an external driving field mimicking a shear profile is analyzed by studying its dynamical evolution in the short-time regime. Starting from high-temperature disordered configurations (fully disordered configurations, FDC), the critical temperature Tc is determined when the order parameter, defined as the absolute value of the transversal spin profile, exhibits a power-law behavior with an exponent that is a combination of some of the critical exponents of the transition. For each value of the shear field magnitude, labeled as gamma, Tc has been estimated and two stages have been found: (1) a growing stage at low values of gamma, where Tc approximately gammapsi and psi=0.52(3), and (2) a saturation regime at large gamma. The same values of Tc(gamma) were found studying the dynamical evolution from the ground-state configuration with all spins pointing in the same direction. By combining the exponents of the corresponding power laws obtained from each initial configuration, the set of critical exponents was calculated. These values, at large external field magnitude, define a critical behavior different from that of the Ising model and of other driven lattice gases.

Monte Carlo results for the Ising model with shear.

We study the kinetics of domain growth in the Ising model with nonconserved dynamics under the action of a stochastic driving field that mimics the action of a shear flow. At late times, we found multistriped configurations with constant transversal size and linear growth in the direction of the flow. In cases with weak shear, a regime characterized by the decreasing of the transversal size is found that could correspond to previous theoretical investigations. This behavior is confirmed by the analysis of the structure factor patterns.

Influence of pinning effects on the electrochemical formation of silver patterns in agarose-containing sols and gels.

The formation of silver patterns via electrolysis from aqueous silver sulfate + x% w/v agarose sol and gel media, with and without supporting electrolyte, in a quasi-two-dimensional (2D) cylindrical cell at room temperature, is utilized as a reference system to investigate the complexity of pinning effects. From pattern morphology and electrochemical data, both delocalized and localized pinning in the bulk dominate the drift of the growth front, depending on the concentration of agarose in the heterogeneous media. Delocalized pinning results from mobile, small agarose aggregates at the growth front and from their accumulation by the front drift. For gels, localized pinning comes from their own percolated structure. A depinning/pinning transition is observed in going from sols to gels. The relative contribution of diffusion and advection in mass-transport-controlled silver electrodeposition depends on the plating bath composition. On the other hand, silver ion attachment to the cathode appears to be interfered with by some screening caused by weakly adsorbed, mobile agarose aggregates at the metal surface without slowing down the rate of the electron-transfer step at the cathode. Their relative contribution of a delocalized, localized pinning and screening effect to a great extent determines the morphology and transition in the growth mode of silver patterns in both media. The analysis of charge and current transients and the corresponding silver pattern morphologies for open and dense radial patterns is made. Results are qualitatively simulated with a novel, rather simple cellular automaton algorithm.

Irreversible phase transitions driven by an oscillatory parameter in a far-from-equilibrium system.

The dynamic response of a forest-fire model to the harmonic variation of an external parameter is studied by means of numerical simulations. Second-order irreversible phase transitions driven by the harmonic input are reported. The location of such transitions depends on both the amplitude and period of the input signal. By means of epidemic studies the relevant critical exponents can be determined, which allow us to place the reported transitions in the universality class of directed percolation. This conclusion is also supported by a field theoretical calculation.