FazBrowse GitHub Viewer
|
Trending
|
URL:
|
Home
Tools:
[Download Repo ZIP]
[View Raw Code]
[Original HTTPS Page]
Python-2/physics/hubble_parameter.py at master · https-github-com-nzysoft/Python-2 · GitHub
Uh oh!
There was an error while loading.
Please reload this page
.
https-github-com-nzysoft
/
Python-2
Public
forked from
TheAlgorithms/Python
Notifications
You must be signed in to change notification settings
Fork
1
Star
0
Code
Pull requests
0
Actions
Projects
Security and quality
0
Insights
Additional navigation options
Code
Pull requests
Actions
Projects
Security and quality
Insights
Expand file tree
Breadcrumbs
Python-2
/
physics
/
hubble_parameter.py
Copy path
More file actions
More file actions
Latest commit
History
History
History
109 lines (84 loc) · 3.24 KB
Breadcrumbs
Python-2
/
physics
/
hubble_parameter.py
Copy path
File metadata and controls
109 lines (84 loc) · 3.24 KB
Raw
Copy raw file
Download raw file
Open symbols panel
Edit and raw actions
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
"""
Title : Calculating the Hubble Parameter
Description : The Hubble parameter H is the Universe expansion rate
in any time. In cosmology is customary to use the redshift redshift
in place of time, becausethe redshift is directily mensure
in the light of galaxies moving away from us.
So, the general relation that we obtain is
H = hubble_constant*(radiation_density*(redshift+1)**4
+ matter_density*(redshift+1)**3
+ curvature*(redshift+1)**2 + dark_energy)**(1/2)
where radiation_density, matter_density, dark_energy are the relativity
(the percentage) energy densities that exist
in the Universe today. Here, matter_density is the
sum of the barion density and the
dark matter. Curvature is the curvature parameter and can be written in term
of the densities by the completeness
curvature = 1 - (matter_density + radiation_density + dark_energy)
Source :
https://www.sciencedirect.com/topics/mathematics/hubble-parameter
"""
def
hubble_parameter
(
hubble_constant
:
float
,
radiation_density
:
float
,
matter_density
:
float
,
dark_energy
:
float
,
redshift
:
float
,
)
->
float
:
"""
Input Parameters
----------------
hubble_constant: Hubble constante is the expansion rate today usually
given in km/(s*Mpc)
radiation_density: relative radiation density today
matter_density: relative mass density today
dark_energy: relative dark energy density today
redshift: the light redshift
Returns
-------
result : Hubble parameter in and the unit km/s/Mpc (the unit can be
changed if you want, just need to change the unit of the Hubble constant)
>>> hubble_parameter(hubble_constant=68.3, radiation_density=1e-4,
... matter_density=-0.3, dark_energy=0.7, redshift=1)
Traceback (most recent call last):
...
ValueError: All input parameters must be positive
>>> hubble_parameter(hubble_constant=68.3, radiation_density=1e-4,
... matter_density= 1.2, dark_energy=0.7, redshift=1)
Traceback (most recent call last):
...
ValueError: Relative densities cannot be greater than one
>>> hubble_parameter(hubble_constant=68.3, radiation_density=1e-4,
... matter_density= 0.3, dark_energy=0.7, redshift=0)
68.3
"""
parameters
=
[
redshift
,
radiation_density
,
matter_density
,
dark_energy
]
if
any
(
0
>
p
for
p
in
parameters
):
raise
ValueError
(
"All input parameters must be positive"
)
if
any
(
1
<
p
for
p
in
parameters
[
1
:
4
]):
raise
ValueError
(
"Relative densities cannot be greater than one"
)
else
:
curvature
=
1
-
(
matter_density
+
radiation_density
+
dark_energy
)
e_2
=
(
radiation_density
*
(
redshift
+
1
)
**
4
+
matter_density
*
(
redshift
+
1
)
**
3
+
curvature
*
(
redshift
+
1
)
**
2
+
dark_energy
)
hubble
=
hubble_constant
*
e_2
**
(
1
/
2
)
return
hubble
if
__name__
==
"__main__"
:
import
doctest
# run doctest
doctest
.
testmod
()
# demo LCDM approximation
matter_density
=
0.3
print
(
hubble_parameter
(
hubble_constant
=
68.3
,
radiation_density
=
1e-4
,
matter_density
=
matter_density
,
dark_energy
=
1
-
matter_density
,
redshift
=
0
,
)
)
Back
|
FazBrowse Home
|
New Git URL