Network simulation using the simulation language for alternate modeling (SLAM 2)

The simulation language for alternate modeling (SLAM 2) is a general purpose language that combines network, discrete event, and continuous modeling capabilities in a single language system. The efficacy of the system's network modeling is examined and discussed. Examples are given of the symbolism that is used, and an example problem and model are derived. The results are discussed in terms of the ease of programming, special features, and system limitations. The system offers many features which allow rapid model development and provides an informative standardized output. The system also has limitations which may cause undetected errors and misleading reports unless the user is aware of these programming characteristics

Equation of State of nuclear matter in a Virial expansion of nucleons and nuclei

We study the equation of state (EOS) of nuclear matter at subnuclear density in a Virial expansion for a nonideal gas. The gas consists of neutrons, protons, alpha particles, and 8980 species of nuclei with $A \ge 12$ and masses from the finite-range droplet model (FRDM). At very low density, the Virial expansion reduces to nuclear statistical equilibrium. At higher density, the Virial results match smoothly to the relativistic mean field results discussed in our previous paper. We tabulate the resulting EOS at over 73,000 grid points in the temperature range $T$ = 0.158 to 15.8 MeV, the density range $n_B$ = 10$^{-8}$ to 0.1 fm$^{-3}$, and the proton fraction range $Y_P$ = 0.05 to 0.56. In the future we plan to match these low density results to our earlier high density mean field results, and generate a full EOS table for use in supernova and neutron star merger simulations. This Virial EOS is exact in the low density limit.Comment: 15 pages, minor changes, Physical Review C in pres

N-fold way simulated tempering for pairwise interaction point processes

Pairwise interaction point processes with strong interaction are usually difficult to sample. We discuss how Besag lattice processes can be used in a simulated tempering MCMC scheme to help with the simulation of such processes. We show how the N-fold way algorithm can be used to sample the lattice processes efficiently and introduce the N-fold way algorithm into our simulated tempering scheme. To calibrate the simulated tempering scheme we use the Wang-Landau algorithm

Spin and orbital valence bond solids in a one-dimensional spin-orbital system: Schwinger boson mean field theory

A generalized one-dimensional $SU(2)\times SU(2)$ spin-orbital model is studied by Schwinger boson mean-field theory (SBMFT). We explore mainly the dimer phases and clarify how to capture properly the low temperature properties of such a system by SBMFT. The phase diagrams are exemplified. The three dimer phases, orbital valence bond solid (OVB) state, spin valence bond solid (SVB) state and spin-orbital valence bond solid (SOVB) state, are found to be favored in respectively proper parameter regions, and they can be characterized by the static spin and pseudospin susceptibilities calculated in SBMFT scheme. The result reveals that the spin-orbit coupling of $SU(2)\times SU(2)$ type serves as both the spin-Peierls and orbital-Peierles mechanisms that responsible for the spin-singlet and orbital-singlet formations respectively.Comment: 6 pages, 3 figure