Unifying Dark Energy & Dark Matter: Investigating the Chaplygin Gas

Click For Summary
SUMMARY

This discussion focuses on the unification of dark energy and dark matter through the Chaplygin gas model, specifically analyzing a second-order differential equation presented in a research paper. The equation in question is 36 (3/2*ab′′ + (1 − w)ab′ − (1 + w)(1 − 3αw)b = 0, which requires two boundary conditions for a complete solution. The user seeks clarification on the necessity of an additional boundary condition beyond b′(aeq) = 0, and discusses the implications of plotting the function b/b(a_{eq}) using Mathematica.

PREREQUISITES
  • Understanding of differential equations, particularly second-order differential equations.
  • Familiarity with the Chaplygin gas model in cosmology.
  • Knowledge of boundary conditions in mathematical physics.
  • Experience with Mathematica for plotting functions and solving equations.
NEXT STEPS
  • Research the implications of boundary conditions in solving second-order differential equations.
  • Explore the Chaplygin gas model and its applications in cosmology.
  • Learn how to use Mathematica for solving differential equations and plotting results.
  • Investigate the relationship between dark energy and dark matter in modern cosmological theories.
USEFUL FOR

Researchers in cosmology, physicists studying dark energy and dark matter, and students interested in advanced mathematical modeling and computational tools like Mathematica.

kptsilva
Messages
33
Reaction score
0
Hello every one!

How many of you have learned about the unification of dark energy and dark matter?
I'm doing my research regarding the Chaplygin gas and I've come across this research paper.
arxiv.org/pdf/grqc/0202064.pdf

At the end, they solve the equation 36 (3/2*ab′′ + (1 − w)ab′ − (1 + w)(1 − 3αw)b = 0
, and only use one boundary which is b′(aeq) = 0. Don't we need another boundary condition since this is a second order DE?


Any knowledge will be much appreciated!
 
Astronomy news on Phys.org
The other initial condition must be the value of [itex]b(a_{eq})[/itex]. But since what is being plotted is [itex]b / b(a_{eq})[/itex], this is naturally one on the graph.
 
thank you for the idea!
I'd like to plot this using mathematica, so first i need solve this for b. So it requires two boundary conditions. I can't use the condition b(aeq)=1 since the plot b/b(aeq) vs. a will not be different from b vs. a, right?
 

Similar threads

  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 24 ·
Replies
24
Views
5K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 7 ·
Replies
7
Views
3K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 13 ·
Replies
13
Views
4K
  • · Replies 11 ·
Replies
11
Views
4K
  • · Replies 46 ·
2
Replies
46
Views
10K
  • · Replies 11 ·
Replies
11
Views
3K