Repository: PlasmaPy/PlasmaPy
Pull request #867 (2020-07-13 22:50:04)

Improve readability of sub-modules on the Formulary documentation page
I felt the current documentation page for the formulary module was difficult to ready.  The sub-module list is rather busy and not alphabetized (see **Current Layout** below), making it hard to look at and find what you want.  So, in the **New Layout** I made the following changes...

1. Simplified the `toctree` names so the sidebar just shows the label for the sub-module and not the full module syntax.
1. Hide the `toctree` directive and created a table with the simplified module name and the module syntax.  (I think this greatly improve the readability)
1. Put the `nbgallery` into an **Examples** sub-section, so the galleries appear in an **Examples** drop-down in the sidebar.

----

### Current Layout

![image](https://user-images.githubusercontent.com/29869348/87360669-7bce5800-c51f-11ea-90a9-d049c090a677.png)

----

### New Layout

![image](https://user-images.githubusercontent.com/29869348/87360722-93a5dc00-c51f-11ea-89bd-d1a35147da08.png)


codecov[bot]: # [Codecov](https://codecov.io/gh/PlasmaPy/PlasmaPy/pull/867?src=pr&el=h1) Report
> Merging [#867](https://codecov.io/gh/PlasmaPy/PlasmaPy/pull/867?src=pr&el=desc) into [master](https://codecov.io/gh/PlasmaPy/PlasmaPy/commit/72443f3649f517f50d63ca6dbf260cbbe8905920&el=desc) will **not change** coverage.
> The diff coverage is `n/a`.

[![Impacted file tree graph](https://codecov.io/gh/PlasmaPy/PlasmaPy/pull/867/graphs/tree.svg?width=650&height=150&src=pr&token=LzZ5mnXVad)](https://codecov.io/gh/PlasmaPy/PlasmaPy/pull/867?src=pr&el=tree)

```diff
@@           Coverage Diff           @@
##           master     #867   +/-   ##
=======================================
  Coverage   96.07%   96.07%           
=======================================
  Files          59       59           
  Lines        5225     5225           
=======================================
  Hits         5020     5020           
  Misses        205      205           
```



------

[Continue to review full report at Codecov](https://codecov.io/gh/PlasmaPy/PlasmaPy/pull/867?src=pr&el=continue).
> **Legend** - [Click here to learn more](https://docs.codecov.io/docs/codecov-delta)
> `Δ = absolute <relative> (impact)`, `ø = not affected`, `? = missing data`
> Powered by [Codecov](https://codecov.io/gh/PlasmaPy/PlasmaPy/pull/867?src=pr&el=footer). Last update [72443f3...843e6cc](https://codecov.io/gh/PlasmaPy/PlasmaPy/pull/867?src=pr&el=lastupdated). Read the [comment docs](https://docs.codecov.io/docs/pull-request-comments).


Patch:
diff --git a/changelog/865.feature.rst b/changelog/865.feature.rst
new file mode 100644
index 00000000..c2fa45f1
--- /dev/null
+++ b/changelog/865.feature.rst
@@ -0,0 +1,27 @@
+Added aliases (with trailing underscores) for parameters in the formulary:
+
+    * `plasmapy.formulary.dimensionless.Reynolds_number` -> `~plasmapy.formulary.dimensionless.Re_`
+    * `plasmapy.formulary.dimensionless.Mag_Reynolds` -> `~plasmapy.formulary.dimensionless.Rm_`
+    * `plasmapy.formulary.drifts.ExB_drift` -> `~plasmapy.formulary.drifts.veb_`
+    * `plasmapy.formulary.drifts.force_drift` -> `~plasmapy.formulary.drifts.vfd_`
+    * `plasmapy.formulary.parameters.mass_density` -> `~plasmapy.formulary.parameters.rho_`
+    * `plasmapy.formulary.parameters.Afven_speed` -> `~plasmapy.formulary.parameters.va_`
+    * `plasmapy.formulary.parameters.ion_sound_speed` -> `~plasmapy.formulary.parameters.cs_`
+    * `plasmapy.formulary.parameters.thermal_speed` -> `~plasmapy.formulary.parameters.vth_`
+    * `plasmapy.formulary.parameters.thermal_pressure` -> `~plasmapy.formulary.parameters.pth_`
+    * `plasmapy.formulary.parameters.kappa_thermal_speed` -> `~plasmapy.formulary.parameters.vth_kappa_`
+    * `plasmapy.formulary.parameters.inertial_length` -> `~plasmapy.formulary.parameters.cwp_`
+    * `plasmapy.formulary.parameters.Hall_parameter` -> `~plasmapy.formulary.parameters.betaH_`
+    * `plasmapy.formulary.parameters.gyrofrequency` -> `~plasmapy.formulary.parameters.oc_`, `~plasmapy.formulary.parameters.wc_`
+    * `plasmapy.formulary.parameters.gyroradius` -> `~plasmapy.formulary.parameters.rc_`, `~plasmapy.formulary.parameters.rhoc_`
+    * `plasmapy.formulary.parameters.plasma_frequency` -> `~plasmapy.formulary.parameters.wp_`
+    * `plasmapy.formulary.parameters.Debye_length` -> `~plasmapy.formulary.parameters.lambdaD_`
+    * `plasmapy.formulary.parameters.Debye_number` -> `~plasmapy.formulary.parameters.nD_`
+    * `plasmapy.formulary.parameters.magnetic_pressure` -> `~plasmapy.formulary.parameters.pmag_`
+    * `plasmapy.formulary.parameters.magnetic_energy_density` -> `~plasmapy.formulary.parameters.ub_`
+    * `plasmapy.formulary.parameters.upper_hybrid_frequency` -> `~plasmapy.formulary.parameters.wuh_`
+    * `plasmapy.formulary.parameters.lower_hybrid_frequency` -> `~plasmapy.formulary.parameters.wlh_`
+    * `plasmapy.formulary.parameters.Bohm_diffusion` -> `~plasmapy.formulary.parameters.DB_`
+    * `plasmapy.formulary.quantum.deBroglie_wavelength` -> `~plasmapy.formulary.quantum.lambdaDB_`
+    * `plasmapy.formulary.quantum.thermal_deBroglie_wavelength` -> `~plasmapy.formulary.quantum.lambdaDB_th_`
+    * `plasmapy.formulary.quantum.Fermi_energy` -> `~plasmapy.formulary.quantum.Ef_`
diff --git a/changelog/867.doc.rst b/changelog/867.doc.rst
new file mode 100644
index 00000000..789e9aff
--- /dev/null
+++ b/changelog/867.doc.rst
@@ -0,0 +1,2 @@
+Improve readability of the `plasmapy.formulary` page by replacing the `toctree`
+list with a cleaner reST table.
diff --git a/docs/formulary/index.rst b/docs/formulary/index.rst
index bfdb7580..42cf1252 100644
--- a/docs/formulary/index.rst
+++ b/docs/formulary/index.rst
@@ -11,20 +11,52 @@ physical quantities helpful for plasma physics.
 
 .. toctree::
    :maxdepth: 1
+   :hidden:
 
-   mathematics
-   dielectric
-   dimensionless
-   dispersionfunction
-   distribution
-   drifts
-   magnetostatics
-   parameters
-   quantum
-   relativity
-   braginskii
-   collisions
-   ionization
+   Classical Transport <braginskii>
+   Collisions <collisions>
+   Dielectrics <dielectric>
+   Dimensionless <dimensionless>
+   Dispersion Relations <dispersionfunction>
+   Distribution Functions <distribution>
+   Drifts <drifts>
+   Ionization <ionization>
+   Magnetostatics <magnetostatics>
+   Mathematics <mathematics>
+   Plasma Parameters <parameters>
+   Quantum Relations <quantum>
+   Relativistic Relations <relativity>
+
+.. table::
+   :widths: 5 16
+
+   +----------------------------------------------------+-----------------------------------------+
+   | :doc:`Classical Transport <./braginskii>`          | `plasmapy.formulary.braginskii`         |
+   +----------------------------------------------------+-----------------------------------------+
+   | :doc:`Collisions <./collisions>`                   | `plasmapy.formulary.collisions`         |
+   +----------------------------------------------------+-----------------------------------------+
+   | :doc:`Dielectrics <./dielectric>`                  | `plasmapy.formulary.dielectric`         |
+   +----------------------------------------------------+-----------------------------------------+
+   | :doc:`Dimensionless <./dimensionless>`             | `plasmapy.formulary.dimensionless`      |
+   +----------------------------------------------------+-----------------------------------------+
+   | :doc:`Dispersion Relations <./dispersionfunction>` | `plasmapy.formulary.dispersionfunction` |
+   +----------------------------------------------------+-----------------------------------------+
+   | :doc:`Distribution Functions <./distribution>`     | `plasmapy.formulary.distribution`       |
+   +----------------------------------------------------+-----------------------------------------+
+   | :doc:`Drifts <./drifts>`                           | `plasmapy.formulary.drifts`             |
+   +----------------------------------------------------+-----------------------------------------+
+   | :doc:`Ionization <./ionization>`                   | `plasmapy.formulary.ionization`         |
+   +----------------------------------------------------+-----------------------------------------+
+   | :doc:`Magnetostatics <./magnetostatics>`           | `plasmapy.formulary.magnetostatics`     |
+   +----------------------------------------------------+-----------------------------------------+
+   | :doc:`Mathematics <./mathematics>`                 | `plasmapy.formulary.mathematics`        |
+   +----------------------------------------------------+-----------------------------------------+
+   | :doc:`Plasma Parameters <./parameters>`            | `plasmapy.formulary.parameters`         |
+   +----------------------------------------------------+-----------------------------------------+
+   | :doc:`Quantum Relations <./quantum>`               | `plasmapy.formulary.quantum`            |
+   +----------------------------------------------------+-----------------------------------------+
+   | :doc:`Relativistic Relations <./relativity>`       | `plasmapy.formulary.relativity`         |
+   +----------------------------------------------------+-----------------------------------------+
 
 The subpackage makes heavy use of `astropy.units.Quantity` for handling
 conversions between different unit systems. This is especially important
@@ -43,13 +75,16 @@ For a general overview of how unit-based input works, take a look at the
 following example:
 
 
+Examples
+========
+
 .. nbgallery::
-    :caption: Examples
 
     /notebooks/physics
 
 Notes for developers
 ====================
+
 Values should be returned as an Astropy Quantity in SI units.
 
 If a quantity has several names, then the function name should be the
diff --git a/plasmapy/formulary/dimensionless.py b/plasmapy/formulary/dimensionless.py
index 4c3f5a5e..f4d297a9 100644
--- a/plasmapy/formulary/dimensionless.py
+++ b/plasmapy/formulary/dimensionless.py
@@ -8,7 +8,14 @@
 highly turbulent, while turbulence is negligible at low Reynolds
 numbers.
 """
-__all__ = ["beta", "quantum_theta", "Reynolds_number", "Mag_Reynolds"]
+__all__ = [
+    "beta",
+    "Mag_Reynolds",
+    "quantum_theta",
+    "Re_",
+    "Reynolds_number",
+    "Rm_",
+]
 
 from astropy import constants
 from astropy import units as u
@@ -109,6 +116,8 @@ def Reynolds_number(
 
         Re = \frac{\rho U L }{\mu}
 
+    **Aliases:** `Re_`
+
     Parameters
     ----------
     rho : `~astropy.units.Quantity`
@@ -163,6 +172,10 @@ def Reynolds_number(
     return Re
 
 
+Re_ = Reynolds_number
+""" Alias to :func:`Reynolds_number`. """
+
+
 @validate_quantities(U={"can_be_negative": True})
 def Mag_Reynolds(U: u.m / u.s, L: u.m, sigma: u.S / u.m) -> u.dimensionless_unscaled:
     r"""
@@ -177,6 +190,8 @@ def Mag_Reynolds(U: u.m / u.s, L: u.m, sigma: u.S / u.m) -> u.dimensionless_unsc
     where :math:`\eta = \frac{1}{\mu_0 \sigma}`
     and :math:`\mu_0` is the permeability of free space.
 
+    **Aliases:** `Rm_`
+
     Parameters
     ----------
     U : `~astropy.units.Quantity`
@@ -225,3 +240,7 @@ def Mag_Reynolds(U: u.m / u.s, L: u.m, sigma: u.S / u.m) -> u.dimensionless_unsc
     eta = 1 / (mu0 * sigma)
     Rm = abs(U * L / eta)
     return Rm
+
+
+Rm_ = Mag_Reynolds
+""" Alias to :func:`Mag_Reynolds`. """
diff --git a/plasmapy/formulary/drifts.py b/plasmapy/formulary/drifts.py
index a518d1bc..f40dd094 100644
--- a/plasmapy/formulary/drifts.py
+++ b/plasmapy/formulary/drifts.py
@@ -1,7 +1,12 @@
 """
 Formulas for calculating particle drifts.
 """
-__all__ = ["ExB_drift", "force_drift"]
+__all__ = [
+    "ExB_drift",
+    "force_drift",
+    "veb_",
+    "vfd_",
+]
 
 import astropy.units as u
 import numpy as np
@@ -14,6 +19,8 @@ def ExB_drift(E: u.V / u.m, B: u.T) -> u.m / u.s:
     r"""
     Calculate the "electric cross magnetic" particle drift.
 
+    **Aliases:** `veb_`
+
     Parameters
     ----------
     E : ~astropy.units.Quantity
@@ -61,11 +68,17 @@ def ExB_drift(E: u.V / u.m, B: u.T) -> u.m / u.s:
     return cross / (B * B).sum(-1)
 
 
+veb_ = ExB_drift
+""" Alias to :func:`ExB_drift`. """
+
+
 @validate_quantities
 def force_drift(F: u.N, B: u.T, q: u.C) -> u.m / u.s:
     r"""
     Calculate the general force drift for a particle in a magnetic field.
 
+    **Aliases:** `vfd_`
+
     Parameters
     ----------
     F : ~astropy.units.Quantity
@@ -110,3 +123,7 @@ def force_drift(F: u.N, B: u.T, q: u.C) -> u.m / u.s:
     """
     cross = np.cross(F.si.value, B.si.value) * F.unit * B.unit
     return cross / (q * (B * B).sum(-1))
+
+
+vfd_ = force_drift
+""" Alias to :func:`force_drift`. """
diff --git a/plasmapy/formulary/parameters.py b/plasmapy/formulary/parameters.py
index e1db4115..d8fc296a 100644
--- a/plasmapy/formulary/parameters.py
+++ b/plasmapy/formulary/parameters.py
@@ -3,24 +3,44 @@
 plasma frequency or Debye length.
 """
 __all__ = [
-    "mass_density",
     "Alfven_speed",
-    "ion_sound_speed",
-    "thermal_speed",
-    "thermal_pressure",
-    "kappa_thermal_speed",
-    "Hall_parameter",
-    "gyrofrequency",
-    "gyroradius",
-    "plasma_frequency",
+    "Bohm_diffusion",
+    "betaH_",
+    "cs_",
+    "cwp_",
+    "DB_",
     "Debye_length",
     "Debye_number",
+    "gyrofrequency",
+    "gyroradius",
+    "Hall_parameter",
     "inertial_length",
-    "magnetic_pressure",
+    "ion_sound_speed",
+    "kappa_thermal_speed",
+    "lambdaD_",
+    "lower_hybrid_frequency",
     "magnetic_energy_density",
+    "magnetic_pressure",
+    "mass_density",
+    "nD_",
+    "oc_",
+    "plasma_frequency",
+    "pmag_",
+    "pth_",
+    "rc_",
+    "rho_",
+    "rhoc_",
+    "thermal_pressure",
+    "thermal_speed",
+    "ub_",
     "upper_hybrid_frequency",
-    "lower_hybrid_frequency",
-    "Bohm_diffusion",
+    "va_",
+    "vth_",
+    "vth_kappa_",
+    "wc_",
+    "wp_",
+    "wlh_",
+    "wuh_",
 ]
 
 import numbers
@@ -80,6 +100,8 @@ def mass_density(
 ) -> u.kg / u.m ** 3:
     """Utility function to merge two possible inputs for particle charge.
 
+    **Aliases:** `rho_`
+
     Parameters
     ----------
     density : ~astropy.units.Quantity
@@ -132,6 +154,10 @@ def mass_density(
     return rho
 
 
+rho_ = mass_density
+""" Alias to :func:`mass_density`. """
+
+
 @check_relativistic
 @validate_quantities(density={"can_be_negative": False})
 def Alfven_speed(
@@ -140,6 +166,8 @@ def Alfven_speed(
     r"""
     Return the Alfvén speed.
 
+    **Aliases:** `va_`
+
     Parameters
     ----------
     B : ~astropy.units.Quantity
@@ -229,6 +257,10 @@ def Alfven_speed(
     return V_A
 
 
+va_ = Alfven_speed
+""" Alias to :func:`Alfven_speed`. """
+
+
 @check_relativistic
 @validate_quantities(
     T_i={"can_be_negative": False, "equivalencies": u.temperature_energy()},
@@ -249,6 +281,8 @@ def ion_sound_speed(
     r"""
     Return the ion sound speed for an electron-ion plasma.
 
+    **Aliases:** `cs_`
+
     Parameters
     ----------
     T_e : ~astropy.units.Quantity
@@ -412,6 +446,10 @@ def ion_sound_speed(
     return V_S
 
 
+cs_ = ion_sound_speed
+""" Alias to :func:`ion_sound_speed`. """
+
+
 # This dictionary defines coefficients for thermal speeds
 # calculated for different methods and values of ndim.
 # Created here to avoid re-instantiating on each call
@@ -439,6 +477,8 @@ def thermal_speed(
     Return the most probable speed for a particle within a Maxwellian
     distribution.
 
+    **Aliases:** `vth_`
+
     Parameters
     ----------
     T : ~astropy.units.Quantity
@@ -560,6 +600,10 @@ def thermal_speed(
     return np.sqrt(coef * k_B * T / m)
 
 
+vth_ = thermal_speed
+""" Alias to :func:`thermal_speed`. """
+
+
 @validate_quantities(
     T={"can_be_negative": False, "equivalencies": u.temperature_energy()},
     n={"can_be_negative": False},
@@ -568,6 +612,8 @@ def thermal_pressure(T: u.K, n: u.m ** -3) -> u.Pa:
     r"""
     Return the thermal pressure for a Maxwellian distribution.
 
+    **Aliases:** `pth_`
+
     Parameters
     ----------
     T : ~astropy.units.Quantity
@@ -609,6 +655,10 @@ def thermal_pressure(T: u.K, n: u.m ** -3) -> u.Pa:
     return n * k_B * T
 
 
+pth_ = thermal_pressure
+""" Alias to :func:`thermal_pressure`. """
+
+
 @check_relativistic
 @validate_quantities(
     T={"can_be_negative": False, "equivalencies": u.temperature_energy()}
@@ -619,6 +669,8 @@ def kappa_thermal_speed(
     r"""Return the most probable speed for a particle within a Kappa
     distribution.
 
+    **Aliases:** `vth_kappa_`
+
     Parameters
     ----------
     T : ~astropy.units.Quantity
@@ -716,6 +768,10 @@ def kappa_thermal_speed(
     return vTh * coeff
 
 
+vth_kappa_ = kappa_thermal_speed
+""" Alias to :func:`kappa_thermal_speed`. """
+
+
 @validate_quantities(
     n={"can_be_negative": False},
     T={"can_be_negative": False, "equivalencies": u.temperature_energy()},
@@ -735,6 +791,8 @@ def Hall_parameter(
 
     All parameters apply to `particle`.
 
+    **Aliases:** `betaH_`
+
     Parameters
     ----------
     n : ~astropy.units.quantity.Quantity
@@ -789,6 +847,10 @@ def Hall_parameter(
     return gyro_frequency / coll_rate
 
 
+betaH_ = Hall_parameter
+""" Alias to :func:`Hall_parameter`. """
+
+
 @validate_quantities(
     validations_on_return={
         "units": [u.rad / u.s, u.Hz],
@@ -799,6 +861,8 @@ def Hall_parameter(
 def gyrofrequency(B: u.T, particle="e-", signed=False, Z=None) -> u.rad / u.s:
     r"""Calculate the particle gyrofrequency in units of radians per second.
 
+    **Aliases:** `oc_`, `wc_`
+
     Parameters
     ----------
     B : ~astropy.units.Quantity
@@ -897,6 +961,13 @@ def gyrofrequency(B: u.T, particle="e-", signed=False, Z=None) -> u.rad / u.s:
     return omega_ci
 
 
+oc_ = gyrofrequency
+""" Alias to :func:`gyrofrequency`. """
+
+wc_ = gyrofrequency
+""" Alias to :func:`gyrofrequency`. """
+
+
 @validate_quantities(
     Vperp={"can_be_nan": True},
     T_i={"can_be_nan": True, "equivalencies": u.temperature_energy()},
@@ -911,6 +982,8 @@ def gyroradius(
 ) -> u.m:
     r"""Return the particle gyroradius.
 
+    **Aliases:** `rc_`, `rhoc_`
+
     Parameters
     ----------
     B : ~astropy.units.Quantity
@@ -1051,6 +1124,13 @@ def gyroradius(
     return r_Li
 
 
+rc_ = gyroradius
+""" Alias to :func:`gyroradius`. """
+
+rhoc_ = gyroradius
+""" Alias to :func:`gyroradius`. """
+
+
 @validate_quantities(
     n={"can_be_negative": False},
     validations_on_return={
@@ -1062,6 +1142,8 @@ def gyroradius(
 def plasma_frequency(n: u.m ** -3, particle="e-", z_mean=None) -> u.rad / u.s:
     r"""Calculate the particle plasma frequency.
 
+    **Aliases:** `wp_`
+
     Parameters
     ----------
     n : ~astropy.units.Quantity
@@ -1155,6 +1237,10 @@ def plasma_frequency(n: u.m ** -3, particle="e-", z_mean=None) -> u.rad / u.s:
     return omega_p.si
 
 
+wp_ = plasma_frequency
+""" Alias to :func:`plasma_frequency`. """
+
+
 @validate_quantities(
     T_e={"can_be_negative": False, "equivalencies": u.temperature_energy()},
     n_e={"can_be_negative": False},
@@ -1163,6 +1249,8 @@ def Debye_length(T_e: u.K, n_e: u.m ** -3) -> u.m:
     r"""Calculate the characteristic decay length for electric fields,
      due to charge screening.
 
+    **Aliases:** `lambdaD_`
+
     Parameters
     ----------
     T_e: ~astropy.units.Quantity
@@ -1223,6 +1311,10 @@ def Debye_length(T_e: u.K, n_e: u.m ** -3) -> u.m:
     return lambda_D
 
 
+lambdaD_ = Debye_length
+""" Alias to :func:`Debye_length`. """
+
+
 @validate_quantities(
     T_e={"can_be_negative": False, "equivalencies": u.temperature_energy()},
     n_e={"can_be_negative": False},
@@ -1231,6 +1323,8 @@ def Debye_number(T_e: u.K, n_e: u.m ** -3) -> u.dimensionless_unscaled:
     r"""Return the number of electrons within a sphere with a radius
     of the Debye length.
 
+    **Aliases:** `nD_`
+
     Parameters
     ----------
     T_e : ~astropy.units.Quantity
@@ -1290,6 +1384,10 @@ def Debye_number(T_e: u.K, n_e: u.m ** -3) -> u.dimensionless_unscaled:
     return N_D
 
 
+nD_ = Debye_number
+""" Alias to :func:`Debye_number`. """
+
+
 @validate_quantities(
     n={"can_be_negative": False},
     validations_on_return={"equivalencies": u.dimensionless_angles()},
@@ -1299,6 +1397,8 @@ def inertial_length(n: u.m ** -3, particle: particles.Particle) -> u.m:
     r"""
     Calculate a charged particle's inertial length.
 
+    **Aliases:** `cwp_`
+
     Parameters
     ----------
     n : ~astropy.units.Quantity
@@ -1356,11 +1456,21 @@ def inertial_length(n: u.m ** -3, particle: particles.Particle) -> u.m:
     return c / omega_p
 
 
+cwp_ = inertial_length
+"""
+Alias to :func:`inertial_length`.
+
+* Name is shorthand for :math:`c / \\omega_p`.
+"""
+
+
 @validate_quantities
 def magnetic_pressure(B: u.T) -> u.Pa:
     r"""
     Calculate the magnetic pressure.
 
+    **Aliases:** `pmag_`
+
     Parameters
     ----------
     B : ~astropy.units.Quantity
@@ -1415,11 +1525,17 @@ def magnetic_pressure(B: u.T) -> u.Pa:
     return (B ** 2) / (2 * mu0)
 
 
+pmag_ = magnetic_pressure
+""" Alias to :func:`magnetic_pressure`. """
+
+
 @validate_quantities
 def magnetic_energy_density(B: u.T) -> u.J / u.m ** 3:
     r"""
     Calculate the magnetic energy density.
 
+    **Aliases:** `ub_`
+
     Parameters
     ----------
     B : ~astropy.units.Quantity
@@ -1474,6 +1590,10 @@ def magnetic_energy_density(B: u.T) -> u.J / u.m ** 3:
     return magnetic_pressure(B)
 
 
+ub_ = magnetic_energy_density
+""" Alias to :func:`magnetic_energy_density`. """
+
+
 @validate_quantities(
     n_e={"can_be_negative": False},
     validations_on_return={
@@ -1486,6 +1606,8 @@ def upper_hybrid_frequency(B: u.T, n_e: u.m ** -3) -> u.rad / u.s:
     r"""
     Return the upper hybrid frequency.
 
+    **Aliases:** `wuh_`
+
     Parameters
     ----------
     B : ~astropy.units.Quantity
@@ -1542,6 +1664,10 @@ def upper_hybrid_frequency(B: u.T, n_e: u.m ** -3) -> u.rad / u.s:
     return omega_uh
 
 
+wuh_ = upper_hybrid_frequency
+""" Alias to :func:`upper_hybrid_frequency`. """
+
+
 @validate_quantities(
     n_i={"can_be_negative": False},
     validations_on_return={
@@ -1554,6 +1680,8 @@ def lower_hybrid_frequency(B: u.T, n_i: u.m ** -3, ion="p+") -> u.rad / u.s:
     r"""
     Return the lower hybrid frequency.
 
+    **Aliases:** `wlh_`
+
     Parameters
     ----------
     B : ~astropy.units.Quantity
@@ -1631,6 +1759,10 @@ def lower_hybrid_frequency(B: u.T, n_i: u.m ** -3, ion="p+") -> u.rad / u.s:
     return omega_lh
 
 
+wlh_ = lower_hybrid_frequency
+""" Alias to :func:`lower_hybrid_frequency`. """
+
+
 @validate_quantities(
     T_e={"can_be_negative": False, "equivalencies": u.temperature_energy()},
     B={"can_be_negative": False},
@@ -1651,6 +1783,8 @@ def Bohm_diffusion(T_e: u.K, B: u.T) -> u.m ** 2 / u.s:
     where :math:`k_B` is the Boltzmann constant
     and :math:`e` is the fundamental charge.
 
+    **Aliases:** `DB_`
+
     Parameters
     ----------
     T_e: `~astropy.units.Quantity`
@@ -1692,3 +1826,7 @@ def Bohm_diffusion(T_e: u.K, B: u.T) -> u.m ** 2 / u.s:
     """
     D_B = k_B * T_e / (16 * e * B)
     return D_B
+
+
+DB_ = Bohm_diffusion
+""" Alias to :func:`Bohm_diffusion`. """
diff --git a/plasmapy/formulary/quantum.py b/plasmapy/formulary/quantum.py
index bca445d2..b36bc7ac 100644
--- a/plasmapy/formulary/quantum.py
+++ b/plasmapy/formulary/quantum.py
@@ -4,12 +4,15 @@
 
 """
 __all__ = [
+    "chemical_potential",
     "deBroglie_wavelength",
-    "thermal_deBroglie_wavelength",
+    "Ef_",
     "Fermi_energy",
+    "lambdaDB_",
+    "lambdaDB_th_",
     "Thomas_Fermi_length",
+    "thermal_deBroglie_wavelength",
     "Wigner_Seitz_radius",
-    "chemical_potential",
 ]
 
 import numpy as np
@@ -31,6 +34,8 @@ def deBroglie_wavelength(V: u.m / u.s, particle) -> u.m:
     r"""
     Calculates the de Broglie wavelength.
 
+    **Aliases:** `lambdaDB_`
+
     Parameters
     ----------
     V : ~astropy.units.Quantity
@@ -128,6 +133,10 @@ def deBroglie_wavelength(V: u.m / u.s, particle) -> u.m:
     return lambda_dBr
 
 
+lambdaDB_ = deBroglie_wavelength
+""" Alias to :func:`deBroglie_wavelength`. """
+
+
 @validate_quantities(
     T_e={"can_be_negative": False, "equivalencies": u.temperature_energy()},
     validations_on_return={"can_be_negative": False},
@@ -136,6 +145,8 @@ def thermal_deBroglie_wavelength(T_e: u.K) -> u.m:
     r"""
     Calculate the thermal deBroglie wavelength for electrons.
 
+    **Aliases:** `lambdaDB_th_`
+
     Parameters
     ----------
     T_e: ~astropy.units.Quantity
@@ -181,6 +192,10 @@ def thermal_deBroglie_wavelength(T_e: u.K) -> u.m:
     return lambda_dbTh
 
 
+lambdaDB_th_ = thermal_deBroglie_wavelength
+""" Alias to :func:`thermal_deBroglie_wavelength`. """
+
+
 @validate_quantities(
     n_e={"can_be_negative": False}, validations_on_return={"can_be_negative": False}
 )
@@ -188,6 +203,8 @@ def Fermi_energy(n_e: u.m ** -3) -> u.J:
     r"""
     Calculate the kinetic energy in a degenerate electron gas.
 
+    **Aliases:** `Ef_`
+
     Parameters
     ----------
     n_e : ~astropy.units.Quantity
@@ -242,6 +259,10 @@ def Fermi_energy(n_e: u.m ** -3) -> u.J:
     return energy_F
 
 
+Ef_ = Fermi_energy
+""" Alias to :func:`Fermi_energy`. """
+
+
 @validate_quantities(
     n_e={"can_be_negative": False}, validations_on_return={"can_be_negative": False}
 )


Test patch:
diff --git a/plasmapy/formulary/tests/test_dimensionless.py b/plasmapy/formulary/tests/test_dimensionless.py
index 34567ae6..910b8f72 100644
--- a/plasmapy/formulary/tests/test_dimensionless.py
+++ b/plasmapy/formulary/tests/test_dimensionless.py
@@ -8,6 +8,8 @@
     Reynolds_number,
     beta,
     quantum_theta,
+    Re_,
+    Rm_,
 )
 from plasmapy.utils import RelativityError
 
@@ -70,3 +72,10 @@ def test_Mag_Reynolds():
 
     with pytest.raises(u.UnitTypeError):
         Mag_Reynolds(2.2 * u.kg, L, sigma)
+
+
+def test_dimensionless_aliases():
+    r"""Test all aliases defined in dimensionless.py"""
+
+    assert Re_ is Reynolds_number
+    assert Rm_ is Mag_Reynolds
diff --git a/plasmapy/formulary/tests/test_drifts.py b/plasmapy/formulary/tests/test_drifts.py
index 48c8231e..11517501 100644
--- a/plasmapy/formulary/tests/test_drifts.py
+++ b/plasmapy/formulary/tests/test_drifts.py
@@ -38,6 +38,9 @@ def test_ExB_3d_array(self):
             result, (2 / 3) * u.Quantity([[[0, 0, 1]]], unit=u.m / u.s)
         )
 
+    def test_alias(self):
+        assert drifts.veb_ is drifts.ExB_drift
+
 
 class Test_force_drift:
     def test_force_x_B_1d_arrays(self):
@@ -75,3 +78,6 @@ def test_nonsensical_units(self):
 
         with pytest.raises(u.UnitTypeError):
             drifts.force_drift(F, B, q)
+
+    def test_alias(self):
+        assert drifts.vfd_ is drifts.force_drift
diff --git a/plasmapy/formulary/tests/test_parameters.py b/plasmapy/formulary/tests/test_parameters.py
index bc5435f8..8d2af7d2 100644
--- a/plasmapy/formulary/tests/test_parameters.py
+++ b/plasmapy/formulary/tests/test_parameters.py
@@ -9,21 +9,42 @@
 from plasmapy.formulary.parameters import (
     Alfven_speed,
     Bohm_diffusion,
+    betaH_,
+    cs_,
+    cwp_,
+    DB_,
     Debye_length,
     Debye_number,
     gyrofrequency,
     gyroradius,
+    Hall_parameter,
     inertial_length,
     ion_sound_speed,
     kappa_thermal_speed,
+    lambdaD_,
     lower_hybrid_frequency,
     magnetic_energy_density,
     magnetic_pressure,
     mass_density,
+    nD_,
+    oc_,
     plasma_frequency,
+    pmag_,
+    pth_,
+    rc_,
+    rho_,
+    rhoc_,
     thermal_pressure,
     thermal_speed,
+    ub_,
     upper_hybrid_frequency,
+    va_,
+    vth_,
+    vth_kappa_,
+    wc_,
+    wp_,
+    wlh_,
+    wuh_,
 )
 from plasmapy.particles.exceptions import InvalidParticleError
 from plasmapy.utils.exceptions import (
@@ -1102,3 +1123,28 @@ def test_Bohm_diffusion():
 
     with pytest.raises(u.UnitTypeError):
         Bohm_diffusion(2.2 * u.kg, B)
+
+
+def test_parameters_aliases():
+    r"""Test all aliases defined in parameters.py"""
+
+    assert rho_ is mass_density
+    assert va_ is Alfven_speed
+    assert cs_ is ion_sound_speed
+    assert vth_ is thermal_speed
+    assert pth_ is thermal_pressure
+    assert vth_kappa_ is kappa_thermal_speed
+    assert betaH_ is Hall_parameter
+    assert oc_ is gyrofrequency
+    assert wc_ is gyrofrequency
+    assert rc_ is gyroradius
+    assert rhoc_ is gyroradius
+    assert wp_ is plasma_frequency
+    assert lambdaD_ is Debye_length
+    assert nD_ is Debye_number
+    assert cwp_ is inertial_length
+    assert pmag_ is magnetic_pressure
+    assert ub_ is magnetic_energy_density
+    assert wuh_ is upper_hybrid_frequency
+    assert wlh_ is lower_hybrid_frequency
+    assert DB_ is Bohm_diffusion
diff --git a/plasmapy/formulary/tests/test_quantum.py b/plasmapy/formulary/tests/test_quantum.py
index ea3e38d0..a968b75a 100644
--- a/plasmapy/formulary/tests/test_quantum.py
+++ b/plasmapy/formulary/tests/test_quantum.py
@@ -13,6 +13,9 @@
     chemical_potential,
     deBroglie_wavelength,
     thermal_deBroglie_wavelength,
+    Ef_,
+    lambdaDB_,
+    lambdaDB_th_,
 )
 
 
@@ -226,3 +229,11 @@ def test_fail1(self):
             f"and should not be equal to {fail1}."
         )
         assert testTrue, errStr
+
+
+def test_quantum_aliases():
+    r"""Test all aliases defined in quantum.py"""
+
+    assert Ef_ is Fermi_energy
+    assert lambdaDB_ is deBroglie_wavelength
+    assert lambdaDB_th_ is thermal_deBroglie_wavelength