DZone
Thanks for visiting DZone today,
Edit Profile
  • Manage Email Subscriptions
  • How to Post to DZone
  • Article Submission Guidelines
Sign Out View Profile
  • Post an Article
  • Manage My Drafts
Over 2 million developers have joined DZone.
Log In / Join
Refcards Trend Reports
Events Video Library
Refcards
Trend Reports

Events

View Events Video Library

Related

  • Python in 2026: uv vs Poetry vs pip: The Definitive Comparison
  • Real-Time Face Recognition Using OpenCV, Dlib, and Python
  • From Polling to PubSub: Building an Asynchronous OPC UA Stack in Python
  • Introducing RAI Audit Kit: Evidence-Grade Responsible AI Audits in Python

Trending

  • Agents, Tools, and MCP: A Mental Model That Actually Helps
  • Retesting Best Practices for Agile Teams: A Quick Guide to Bug Fix Verification
  • Agentic Testing: Moving Quality From Checkpoint to Control Layer
  • Differential Flamegraphs in Java in Jeffrey Microscope
  1. DZone
  2. Coding
  3. Languages
  4. Physical Constants in Python

Physical Constants in Python

Python and SciPy library are pretty powerful scientific tools. In this post, we put this to the test using a famous problem from the world of physics.

By 
John Cook user avatar
John Cook
·
Oct. 05, 18 · Tutorial
Likes (2)
Comment
Save
Tweet
Share
9.1K Views

Join the DZone community and get the full member experience.

Join For Free

You can find a large collection of physical constants in scipy.constants. The most frequently used constants are available directly, and hundreds more are in a dictionary physical_constants.

The fine structure constant α is defined as a function of other physical constants:

The following code shows that the fine structure constant and the other constants that go into it are available in scipy.constants.

    import scipy.constants as sc

    a = sc.elementary_charge**2
    b = 4 * sc.pi * sc.epsilon_0 * sc.hbar * sc.c
    assert( abs(a/b - sc.fine_structure) < 1e-13 )

Eddington's Constant

In the 1930s, Arthur Eddington believed that the number of photons in the observable universe was exactly the Eddington number:

Since at the time the fine structure constant was thought to be 1/136, this made the number of photons a nice even 136 × 2 256. Later he revised his number when it looked like the fine structure constant was 1/137. According to the Python code above, the current estimate is more like 1/137.036.

Eddington was a very accomplished scientist, though he had some ideas that seem odd today. His number is a not a bad estimate, though nobody believes it could be exact.

Related Posts

The constants in scipy.constants have come up in a couple previous blog posts.

The post on Koide's coincidence shows how to use the physical_constants dictionary, which includes not just the physical constant values but also their units and uncertainty.

The post on Benford's law shows that the leading digits of the constants in scipy.constants follow the logarithmic distribution observed by Frank Benford (and earlier by Simon Newcomb).

Python (language)

Published at DZone with permission of John Cook. See the original article here.

Opinions expressed by DZone contributors are their own.

Related

  • Python in 2026: uv vs Poetry vs pip: The Definitive Comparison
  • Real-Time Face Recognition Using OpenCV, Dlib, and Python
  • From Polling to PubSub: Building an Asynchronous OPC UA Stack in Python
  • Introducing RAI Audit Kit: Evidence-Grade Responsible AI Audits in Python

Partner Resources

×

Comments

The likes didn't load as expected. Please refresh the page and try again.

  • RSS
  • X
  • Facebook

ABOUT US

  • About DZone
  • Support and feedback
  • Community research

ADVERTISE

  • Advertise with DZone

CONTRIBUTE ON DZONE

  • Article Submission Guidelines
  • Become a Contributor
  • Core Program
  • Visit the Writers' Zone

LEGAL

  • Terms of Service
  • Privacy Policy

CONTACT US

  • 3343 Perimeter Hill Drive
  • Suite 215
  • Nashville, TN 37211
  • [email protected]

Let's be friends:

  • RSS
  • X
  • Facebook