Monte Carlo Study of the Growth of L12 ordered domains in fcc A3B binary alloys

A Monte Carlo study of the late time growth of L12-ordered domains in a fcc A3B binary alloy is presented. The energy of the alloy has been modeled by a nearest-neighbor interaction Ising Hamiltonian. The system exhibits a fourfold degenerated ground state and two kinds of interfaces separating ordered domains: flat and curved antiphase boundaries. Two different dynamics are used in the simulations: the standard atom-atom exchange mechanism and the more realistic vacancy-atom exchange mechanism. The results obtained by both methods are compared. In particular we study the time evolution of the excess energy, the structure factor and the mean distance between walls. In the case of atom-atom exchange mechanism anisotropic growth has been found: two characteristic lengths are needed in order to describe the evolution. Contrarily, with the vacancyatom exchange mechanism scaling with a single length holds. Results are contrasted with existing experiments in Cu3Au and theories for anisotropic growth.

Modeling premartensitic effects in Ni2MnGa: A mean-field and Monte Carlo simulation study

ic first-order transition line ending in a critical point. This critical point is responsible for the existence of large premartensitic fluctuations which manifest as broad peaks in the specific heat, not always associated with a true phase transition. The main conclusion is that premartensitic effects result from the interplay between the softness of the anomalous phonon driving the modulation and the magnetoelastic coupling. In particular, the premartensitic transition occurs when such coupling is strong enough to freeze the involved mode phonon. The implication of the results in relation to the available experimental data is discussed.

Comparison between molecular dynamics abd Monte Carlo simulations of an ordering process in a binary alloy

Ordering in a binary alloy is studied by means of a molecular-dynamics (MD) algorithm which allows to reach the domain growth regime. Results are compared with Monte Carlo simulations using a realistic vacancy-atom (MC-VA) mechanism. At low temperatures fast growth with a dynamical exponent x>1/2 is found for MD and MC-VA. The study of a nonequilibrium ordering process with the two methods shows the importance of the nonhomogeneity of the excitations in the system for determining its macroscopic kinetics.

Modeling changes in organic carbon stocks for distinct soils in southeastern brazil after four eucalyptus rotations using the century model

Soil organic matter (SOM) plays an important role in carbon (C) cycle and soil quality. Considering the complexity of factors that control SOM cycling and the long time it usually takes to observe changes in SOM stocks, modeling constitutes a very important tool to understand SOM cycling in forest soils. The following hypotheses were tested: (i) soil organic carbon (SOC) stocks would be higher after several rotations of eucalyptus than in low-productivity pastures; (ii) SOC values simulated by the Century model would describe the data better than the mean of observations. So, the aims of the current study were: (i) to evaluate the SOM dynamics using the Century model to simulate the changes of C stocks for two eucalyptus chronosequences in the Rio Doce Valley, Minas Gerais State, Brazil; and (ii) to compare the C stocks simulated by Century with the C stocks measured in soils of different Orders and regions of the Rio Doce Valley growing eucalyptus. In Belo Oriente (BO), short-rotation eucalyptus plantations had been cultivated for 4.0; 13.0, 22.0, 32.0 and 34.0 years, at a lower elevation and in a warmer climate, while in Virginópolis (VG), these time periods were 8.0, 19.0 and 33.0 years, at a higher elevation and in a milder climate. Soil samples were collected from the 0-20 cm layer to estimate C stocks. Results indicate that the C stocks simulated by the Century model decreased after 37 years of poorly managed pastures in areas previously covered by native forest in the regions of BO and VG. The substitution of poorly managed pastures by eucalyptus in the early 1970´s led to an average increase of C of 0.28 and 0.42 t ha-1 year-1 in BO and VG, respectively. The measured C stocks under eucalyptus in distinct soil Orders and independent regions with variable edapho-climate conditions were not far from the values estimated by the Century model (root mean square error - RMSE = 20.9; model efficiency - EF = 0.29) despite the opposite result obtained with the statistical procedure to test the identity of analytical methods. Only for lower soil C stocks, the model over-estimated the C stock in the 0-20 cm layer. Thus, the Century model is highly promising to detect changes in C stocks in distinct soil orders under eucalyptus, as well as to indicate the impact of harvest residue management on SOM in future rotations.

Calibration of the century, apsim and ndicea models of decomposition and n mineralization of plant residues in the humid tropics

The aim of this study was to calibrate the CENTURY, APSIM and NDICEA simulation models for estimating decomposition and N mineralization rates of plant organic materials (Arachis pintoi, Calopogonium mucunoides, Stizolobium aterrimum, Stylosanthes guyanensis) for 360 days in the Atlantic rainforest bioma of Brazil. The models´ default settings overestimated the decomposition and N-mineralization of plant residues, underlining the fact that the models must be calibrated for use under tropical conditions. For example, the APSIM model simulated the decomposition of the Stizolobium aterrimum and Calopogonium mucunoides residues with an error rate of 37.62 and 48.23 %, respectively, by comparison with the observed data, and was the least accurate model in the absence of calibration. At the default settings, the NDICEA model produced an error rate of 10.46 and 14.46 % and the CENTURY model, 21.42 and 31.84 %, respectively, for Stizolobium aterrimum and Calopogonium mucunoides residue decomposition. After calibration, the models showed a high level of accuracy in estimating decomposition and N- mineralization, with an error rate of less than 20 %. The calibrated NDICEA model showed the highest level of accuracy, followed by the APSIM and CENTURY. All models performed poorly in the first few months of decomposition and N-mineralization, indicating the need of an additional parameter for initial microorganism growth on the residues that would take the effect of leaching due to rainfall into account.

Potencial de uso do modelo century e sig para avaliar o impacto da agricultura sobre estoques regionais de carbono orgânico do solo

O uso e as mudanças no uso e manejo dos solos figuram entre os principais fatores determinantes dos estoques de carbono orgânico do solo (COS). A modelagem dinâmica espacialmente explícita é uma técnica que vem sendo usada com sucesso em avaliações das alterações nos estoques regionais de COS. Assim, este trabalho objetivou utilizar o modelo Century associado a técnicas de geoprocessamento e sensoriamento remoto para avaliar os efeitos das alterações do uso agrícola e manejo do solo na dinâmica do carbono orgânico do solo (COS) em propriedades rurais do Distrito Santana, Ijuí - RS. Para isso, foi criada uma base de dados geoespaciais em ambiente ArcGis9.x com os planos temáticos: solos, elementos da paisagem e início do uso agrícola. Após reconstituídos os cenários históricos de manejo, o modelo Century 4.0 foi inicializado com dados edafoclimáticos locais, calibrado em duas etapas e validado, obedecendo-se à seguinte ordem: calibração principal - validação - calibração para generalizações. Na calibração principal, foram ajustados parâmetros internos do modelo, valendo-se da adição de carbono (C) pelo milho, trigo e soja e dos estoques de COS medidos em 2007 na camada de 0 a 20 cm da mata nativa e de uma lavoura adjacente, cuja conversão ocorreu entre 1901 e 1930 (lavoura mais antiga, localizada na área 1), sob Latossolo de topo. Em seguida, o Century, assim ajustado, foi estatisticamente validado com base em dados observados em 10 lavouras iniciadas em diferentes épocas, sob a mesma classe de solo e posição na paisagem, amostradas em 2007 na camada de 0 a 20 cm e distribuídas em quatro áreas homogêneas (áreas 1, 2, 3 e 4) dentro do Distrito Santana. Na calibração para generalizações, foram criadas quatro opções de "floresta subtropical" para representar a mata nativa de cada classe de solo em estudo (Latossolo, Chernossolo, Neossolo Regolítico e Neossolo Flúvico), a partir da opção de floresta tropical gerada na etapa de calibração principal, mantendo-se os ajustes já realizados naquela etapa e alterando apenas o parâmetro referente à máxima produção bruta mensal de biomassa pela floresta para que a variável de saída SOMTC (COS simulado) se estabilizasse, após o modelo ser executado por 3.000 anos (execução de equilíbrio), coincidindo com os estoques de COS medidos no solo sob vegetação nativa respectiva de cada classe de solo; em seguida, procedeu-se às simulações correspondentes a 54 unidades de simulação (US), por meio da interface i-Century, nas condições de solos cultivados com os cenários de manejo históricos (de 1901 até 2007). O Century, depois de calibrado, estimou adequadamente a dinâmica do C, reproduzindo a evolução dos estoques de COS ocorrida no Distrito Santana. Segundo essas estimativas, a adoção de práticas conservacionistas de manejo garante a redução das perdas de C causadas pelo processo erosivo e decomposição microbiana, havendo com isso aumento no estoque de COS.

Monte Carlo simulation of kilovolt electron transport in solids

A Monte Carlo procedure to simulate the penetration and energy loss of low¿energy electron beams through solids is presented. Elastic collisions are described by using the method of partial waves for the screened Coulomb field of the nucleus. The atomic charge density is approximated by an analytical expression with parameters determined from the Dirac¿Hartree¿Fock¿Slater self¿consistent density obtained under Wigner¿Seitz boundary conditions in order to account for solid¿state effects; exchange effects are also accounted for by an energy¿dependent local correction. Elastic differential cross sections are then easily computed by combining the WKB and Born approximations to evaluate the phase shifts. Inelastic collisions are treated on the basis of a generalized oscillator strength model which gives inelastic mean free paths and stopping powers in good agreement with experimental data. This scattering model is accurate in the energy range from a few hundred eV up to about 50 keV. The reliability of the simulation method is analyzed by comparing simulation results and experimental data from backscattering and transmission measurements.

Elastic constants of the beta-Cu-Zn alloy system: A Monte Carlo study

Monte Carlo calculations of the isothermal elastic constants of the beta-CuxZn1-x alloy system as a function of the composition have been carried out. We assume the atoms interact via a two-body Morse potential function and use numerical values for the potential parameters evaluated taking into account experimental data. We find a quite good agreement between our results and the expected experimental behavior.

Closure of the Monte Carlo dynamical equation in the spherical Sherrington-Kirkpatrick model

We study the analytical solution of the Monte Carlo dynamics in the spherical Sherrington-Kirkpatrick model using the technique of the generating function. Explicit solutions for one-time observables (like the energy) and two-time observables (like the correlation and response function) are obtained. We show that the crucial quantity which governs the dynamics is the acceptance rate. At zero temperature, an adiabatic approximation reveals that the relaxational behavior of the model corresponds to that of a single harmonic oscillator with an effective renormalized mass.

Magnetic microstructures from magnetic force microscopy and Monte Carlo simulation in CoFe-Ag-Cu granular films

CoFe-Ag-Cu granular films, prepared by rf sputtering, displayed magnetic domain microstructures for ferromagnetic concentrations above about 32% at, and below the percolation threshold. All samples have a fcc structure with an (111) texture perpendicular to the film plane. Magnetic force microscopy (MFM) showed a variety of magnetic domain microstructures, extremely sensitive to the magnetic history of the sample, which arise from the balance of the ferromagnetic exchange, the dipolar interactions and perpendicular magnetocrystalline anisotropy, MFM images indicate that in virgin samples, magnetic bubble domains with an out-of-plane component of the magnetization are surrounded by a quasicontinuous background of opposite magnetization domains. The application of a magnetic field in different geometries drastically modifies the microstructure of the system in the remanent state: i) for an in-plane field, the MFM images show that most of the magnetic moments are aligned along the film plane, ii) for an out-of-plane field, the MFM signal increases about one order of magnitude, and out-of-plane striped domains with alternating up and down magnetization are stabilized. Numerical simulations show that a variety of metastable domain structures (similar to those observed experimentally) can be reached, depending on magnetic history, in systems with competing perpendicular anisotropy, exchange and dipolar interactions.