FleCSPH is a smoothed particle hydrodynamics simulation tool, based on the compile-time configurable framework FleCSI. The asynchronous distributed tree topology combined with a fast multipole method allows FleCSPH to efficiently compute hydrodynamics and long range particle-particle interactions. FleCSPH provides initial data generators, particle relaxation techniques, and standard evolution drivers, which can be easily modified and extended to user-specific setups. Data input/output uses the H5part format, compatible with modern visualization software.

We present new results for ground-state candidate energies of Mg-rich olivine (MRO) clusters and use the binding energies of these clusters to determine their nucleation rates in stellar outflows, with particular interest in the en- vironments of core-collapse supernovae (CCSNe). Low-lying structures of clusters (Mg2SiO4)n 2 \le n \le 13 are determined from a modified minima hopping algorithm using an empirical silicate potential in the Buckingham form. These configurations are further refined and optimized using the den- sity functional theory code Quantum Espresso. Utilizing atomistic nucleation theory, we determine the critical size and nucleation rates of these clusters. We find that configurations and binding energies in this regime are very dis- similar from those of the bulk lattice. Clusters grow with SiO_4-MgO layering and exhibit only global, rather than local, symmetries. When compared to classical nucleation theory we find suppressed nucleation rates at most temperatures and pressures, with enhanced nucleation rates at very large pressures. This implies a slower progression of silicate dust formation in stellar environments than previously assumed.

The classical theory of grain nucleation suffers from both theoretical and predictive deficiencies. We strive to alleviate these deficiencies in our understanding of dust formation and growth by utilizing an atomistic model of nucleation. Carbon cluster geometries are determined with a set of global minimization algorithms. Using density functional theory, the binding energies of carbon clusters from $n=2$ to $n=99$ are then calculated. These energies are used to calculate the critical size and nucleation rate of carbon clusters. We find that the critical cluster size is largely determined by the changes in geometry of the clusters. Clusters with size $n=27$ and $n=8$, roughly corresponding to the transition from ring-to-fullerene geometry and chain-to-ring geometry respectively, are the critical sizes across the range of temperature and saturation where nucleation is significant. In contrast to the classical theory, nucleation is enhanced at low-temperatures, and suppressed at high temperatures. These results will be applied to a modified chemical evolution code using results from supernova simulations.

The binding energies of n < 100 carbon clusters are calculated using the ab-initio density functional theory code Quantum Espresso. Carbon cluster geometries are determined using several levels of classical techniques and further refined using density functional theory. The resulting energies are used to compute the work of cluster formation and the nucleation rate in a saturated, hydrogen-poor carbon gas. Compared to classical calculations that adopt the capillary approximation, we find that nucleation of carbon clusters is enhanced at low temperatures and depressed at high temperatures. This difference is ascribed to the different behavior of the critical cluster size. We find that the critical cluster size is at n = 27 or n = 8 for a broad range of temperatures and saturations, instead of being a smooth function of such parameters. The results of our calculations can be used to follow carbonaceous cluster/grain formation, stability, and growth in hydrogen poor environments, such as the inner layers of core-collapse supernova and supernova remnants.

Ascertaining the core-collapse supernova mechanism is a complex, and yet unsolved, problem dependent on the interaction of general relativity, hydrodynamics, neutrino transport, neutrino-matter interactions, and nuclear equations of state and reaction kinetics. Ab initio modeling of core-collapse supernovae and their nucleosynthetic outcomes requires care in the coupling and approximations of the physical components. We have built our multi-physics CHIMERA code for supernova modeling in 1-, 2-, and 3-D, using ray-by-ray neutrino transport, approximate general relativity, and detailed neutrino and nuclear physics. We discuss some early results from our current series of exploding 2D simulations and our work to perform computationally tractable simulations in 3D using the "Yin-Yang" grid.