Purchasing new hardware? Read our latest product comparisons

Simulated universe spawns the whole array of galaxies


August 22, 2012

Left shows galaxies from AREPO simulation, right shows actual galaxies from Hubble image (Image: Harvard-Smithsonian Center for Astrophysics and NASA)

Left shows galaxies from AREPO simulation, right shows actual galaxies from Hubble image (Image: Harvard-Smithsonian Center for Astrophysics and NASA)

Image Gallery (9 images)

A new approach for simulating the birth and evolution of galaxies and cosmic filaments within the Universe has been developed by researchers at the Harvard-Smithsonian Center for Astrophysics together with their colleagues at the Heidelberg Institute for Theoretical Studies. It is called AREPO, and has been used to simulate the evolution of our Universe from only 380,000 years after the Big Bang to the present. The full variety of spiral, elliptical, peculiar, and dwarf galaxies appear in the simulated Universe.

Scientists have been using computers to simulate the evolution of the Universe for about two decades. The cosmological model used for most simulations is the Lambda Cold Dark Matter (LCDM) model, which is considered the standard model for Big Bang cosmology. The model includes general relativity with a cosmological constant, and a predominant portion of dark matter and energy. It successfully describes the existence and structure of the cosmic microwave background, the large scale structure in the distribution of galaxies, the abundances of the light elements, and the accelerating expansion of the Universe. However, hitching the LCDM model to a 14 billion year simulation leading to the present of our Universe is not a trivial task.

The quality of such simulations has increased rapidly, owing as much to refinement of the simulation techniques as to the titanic increase in computational power over the past 20 years. High-accuracy simulations of the behavior of dark matter and energy in the LCDM model have been carried out. However, these are relatively simple simulations in that dark matter and energy are only thought to interact through their gravity. Including normal matter and energy is considerably more complex, because they do interact in ways that are important influences on the evolution of the Universe.

To model the formation of galactic structure requires treatment of many physical phenomena beyond that of gravity, including the effects of magnetism, hydrodynamics, radiation, multiple types of matter and energy, and shock waves among others. In most cases our knowledge of the physics underlying these phenomena is sufficient, but how to best simulate their mutual interactions on a computer is still being improved.

The striking spiral structure of galaxy M81 in the Big Dipper (Photo: NASA)

A glaring example is offered by galactic structure. For example, our cosmic neighborhood is full of large spiral galaxies like our Milky Way, and the Andromeda galaxy (M31). Such spirals are common, but cosmological evolution simulations do not generally produce universes containing large spiral galaxies. Rather they produce clumps of matter making up roughly spherical amorphous galaxies without anything like the broad disks and extended arms of a typical spiral galaxy. This is shown clearly in the following video where the gas distribution is shown for a galaxy from a cosmological simulation using the codes GADGET (bottom) and AREPO (top). The nested views show similar large scale structure (left panel), but reveal very different gas distributions around galaxies (right panel).

Two basic classes of simulation methods have been used to simulate the evolution of the Universe. These are Lagrangian Smoothed Particle Hydrodynamics (SPH) and Eulerian mesh-based hydrodynamics with Adaptive Mesh Refinement (AMR). Both of these methods have serious limitations for handling such complex physics. SPH treats shock waves and discontinuities poorly (a great deal of structure evolution takes place through collision of smaller chunks of matter), and suppresses the hydrodynamic instabilities which are largely responsible for growth of rotation in interacting bodies. The GADGET code referred to above is an SPH code.

AMR also has significant limitations, despite a long history in simulating hydrodynamics in many real world problems. They violate Galilean-invariance, which simply means that the results change if a massive body is moving across spatial coordinates instead of being at a constant position. This is not how physics works. The physical result should not depend on how an observer moves. For the same reason, a highly structured moving body will tend to lose its structure simply as a result of its motion in the simulation. Particularly unwelcome behavior is associated with AMR simulation of matter flows where the kinetic energy is large compared to the thermal energy: the usual situation in galactic evolution. This sort of non-physical prediction hampers one's trust in AMR simulation.

AREPO was designed to combine the best parts of SPH and AMR simulation methods while leaving behind the troublesome parts. The basic difference is in the spatial model for the normal and dark matter and energy in the Universe. The initial condition is split up into an unstructured set of smaller volumes by a mesh. AREPO uses a mesh that flexes and moves in space to match the motions of the underlying gas, stars, dark matter, and dark energy. The AREPO mesh can continuously and consistently adjust to local conditions to insure proper treatment of physics associated with normal matter and energy. AREPO retains the benefits of the SPH approach while avoiding most of the inaccuracies. It also retains the high-accuracy treatment of shocks, shear waves, and fluid instabilities characteristic of AMR. The devil is in the details, but the simulation approach used in AREPO is a great advance over previous methods.

The AREPO simulation of the Universe ran on Harvard's Odyssey high-performance supercomputer, using in total 1024 processor cores. The only information used as input to the LCDM model was the distribution of the 3K background radiation. Fourteen billion years was compressed into only a few months, leading to the video below, which gives a tour through the the last ten billion years of the evolution of the Universe ... or at least of a universe.

Source: Harvard-Smithsonian Center for Astrophysics

About the Author
Brian Dodson From an early age Brian wanted to become a scientist. He did, earning a Ph.D. in physics and embarking on an R&D career which has recently broken the 40th anniversary. What he didn't expect was that along the way he would become a patent agent, a rocket scientist, a gourmet cook, a biotech entrepreneur, an opera tenor and a science writer. All articles by Brian Dodson

Why do people keep repeating this story about the "big bang" theory when it was shown back in the 1960s or 70s, that after 20 years of comparing the spectral lines of many thousands of stars, no two of them were the same?? They were as different as snowflakes according to the guy who was doing the research. I saw it in a movie in high school about Alvin Toffler's book Future Shock. If all of the stars we could see at that time were made out of one homogenous glob of matter that suddenly and inexplicably decided to expand, then why aren't all of the stars giving off the same exact spectral lines? I think that you could at least find a couple of them that were pretty close in that regard, right? This is like saying that milk from the bottom of the carton is different from the milk at the top of it. And even MORE important than that, from WHERE did that initial glob of mass come from that suddenly had a hissy fit and blew itself to pieces?!?! And taking it out to the Ne Plus Ultra a bit more, how far out does space go and what's on the other side of the boundary?? Surely there must be answers to these questions, why aren't they being asked by anyone? Or, is space really just a viewpoint (consideration) of dimension, and what is real is what we AGREE upon to be real??

Your ever humble logician servant, Randy

Expanded Viewpoint

they have forgotten some vital data. the way i see it a bang everything spreads out. then as the things slow they get caught in gravity wells of larger items. the first item to be caught in a gravity well will help shape the system bcs its gravity well will shape the obits of those that get caught later. and the only degrees i have or for cooking and for drafting.

Hahn Jackson
Post a Comment

Login with your Gizmag account:

Related Articles
Looking for something? Search our articles