Looking for Guidance on Field-Cooling

[Pravin Sharma]

Hello VAMPIRE Users,

 I’m working on the field-cooling calculations, and I’m a bit unsure about the exact quantities we should extract and plot for proper analysis. Could someone please point me in the right direction or share the key variables typically used for field-cooling plots?

I can save the .chk files, which I will use later for hysteresis calculations. However, I am confused about what to plot for field cooling? 

I am sharing the input files, which I am using to calculate the field-cooling simulations.

#------------------------------------------

# Creation attributes

#------------------------------------------

create:crystal-structure = bcc

create:periodic-boundaries-x

create:periodic-boundaries-y

#create:periodic-boundaries-z

#------------------------------------------

# System Dimensions

#----------------------------------------

dimensions:unit-cell-size = 2.82 !A

dimensions:system-size-x = 20.0 !nm

dimensions:system-size-y = 20.0 !nm

dimensions:system-size-z = 10.0 !nm

#------------------------------------------

# Material Files:

#------------------------------------------

material:file = bilayer.mat

#------------------------------------------

# Simulation attributes:

#------------------------------------------

#sim:temperature = 0.0

#sim:equilibration-temperature = 50.0

sim:minimum-temperature = 50.0

sim:maximum-temperature = 300.0

sim:temperature-increment = 5.0

sim:equilibration-time-steps = 50000

sim:loop-time-steps = 100000

sim:time-steps-increment = 10

sim:time-step = 1e-15

#extraline

sim:applied-field-strength = 1.0 !T

sim:applied-field-unit-vector = 0.999,0.001,0  

sim:cooling-time = 50 !ps

sim:cooling-function = gaussian

sim:integrator-random-seed = 1000

sim:save-checkpoint = end

#------------------------------------------

# Program and integrator details

#------------------------------------------

sim:program = field-cool

sim:integrator = monte-carlo

#------------------------------------------

# Data output

#------------------------------------------

output:column-headers

output:fixed-width = 0.1

output:temperature

output:magnetisation

output:mean-magnetisation-length

output:magnetisation-length

config:atoms

config:atoms-output-rate = 1000

Best regards,
Pravin

VAMPIRE User