Repository: PlasmaPy/PlasmaPy

Issue:
Create aliases for common formulary functions.
### Update (2020-12-17)

Most of the common `plasmapy.formulary` functionality has been given aliases.  This issue is being left open since aliases may be desired on future functionality.  If there's functionality that does not have an alias and you'd like one, then please leave a comment.

----

Explicit function names are great for readability and clarity, and the primary functions should remain like that.  However, explicit names can be long and become burdensome for users.  Users either work with the long names or have to individually import each function and rename it to a more manageable name.  Thus, I propose we make aliases for, at least, the most common formulary functionality.  Implementation would looks something like...

```python
__all__ = ["Alfven_speed", "va_"]

...

def Alfven_speed(...):
    ...
    return speed

va_ = Alfven_speed
"""Alias to :func:`Alfven_speed`."""

# or

#: Alias for :func:`Alfven_speed`.
va_ = Alfven_speed
```
### Development points

* Aliases will have a trailing underscore to signify they are aliases and allow users to still use the shortest abbreviations for their code. For example,

  ```python
  from plasmapy.formulary import va_

  va = va_(...)
  ```
* Aliases shout be documented in the docstring of the parent function after the main description but before the parameters list.

  ```python
  def Alfven_speed(...):
      """
      Main function description

      **Aliases:**  `va_`

      Parameters
      ----------
      ...
      """
  ```

* Aliases should be tested in the the parent function `pytest`s.  This can be done with the following simple `assert` statement.

  ```python
  assert va_ is Alfven_speed
  ```

### List of functions to make aliases for...

(By no means is this a complete list and will be updated as needed.)

- [x] `Alfven_speec` -> `va_` (PR #865)
- [x] `Bohm_diffusion` -> `DB_`  (PR #865)
- [x] `deBroglie_wavelength` -> `lambdaDB_` (PR #865)
- [x] `Debye_length` -> `lambdaD_` (PR #865)
- [x] `Debye_number` -> `nD_` (PR #865)
- [x] `ExB_drift` -> `veb_` (PR #865)
- [x] `Fermi_energy` -> `Ef_` (PR #865)
- [x] `force_drift` -> `vfd_` (PR #865)
- [x] `gyrofrequency` -> `oc_` `wc_` (PR #865)
- [x] `gyroradius` -> `rc_` `rhoc_` (PR #865)
- [x] `inertial_length` -> `cwp_` (PR #865)
- [x] `ion_sound_speed` -> `cs_` (PR #865)
- [x] `ionization_balance` -> `Z_bal_` (PR #860, #851)
- [x] `kappa_thermal_speed` -> `vth_kappa_` (PR #865)
- [x] `lower_hybrid_frequency` -> `wlh_` (PR #865)
- [x] `Mag_Reynolds` -> `Rm_` (PR #865)
- [x] `magnetic_energy_density` -> `ub_` (PR #865)
- [x] `magnetic_pressure` -> `pmag_` (PR #865)
- [x] `mass_density` -> `rho_` (PR #865)
- [x] `plasma_frequency` = `wp_` (PR #865)
- [x] `Reynolds_number` -> `Re_` (PR #865)
- [x] `thermal_deBroglie_wavelength` -> `lambdaDB_th_` (PR #865)
- [x] `thermal_speed` -> `vth_` (PR #865)
- [x] `thermal_pressure` -> `pth_` (PR #865)
- [x] `upper_hybrid_frequency` -> `wuh_` (PR #865)
- [ ] ...

Discussion:
codecov[bot]: # [Codecov](https://codecov.io/gh/PlasmaPy/PlasmaPy/pull/865?src=pr&el=h1) Report
> Merging [#865](https://codecov.io/gh/PlasmaPy/PlasmaPy/pull/865?src=pr&el=desc) into [master](https://codecov.io/gh/PlasmaPy/PlasmaPy/commit/511c25bce06fda6e393e03fc1d422901e8779591&el=desc) will **increase** coverage by `0.02%`.
> The diff coverage is `100.00%`.

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

```diff
@@            Coverage Diff             @@
##           master     #865      +/-   ##
==========================================
+ Coverage   96.07%   96.09%   +0.02%     
==========================================
  Files          59       59              
  Lines        5217     5245      +28     
==========================================
+ Hits         5012     5040      +28     
  Misses        205      205              
```


| [Impacted Files](https://codecov.io/gh/PlasmaPy/PlasmaPy/pull/865?src=pr&el=tree) | Coverage Δ | |
|---|---|---|
| [plasmapy/formulary/dimensionless.py](https://codecov.io/gh/PlasmaPy/PlasmaPy/pull/865/diff?src=pr&el=tree#diff-cGxhc21hcHkvZm9ybXVsYXJ5L2RpbWVuc2lvbmxlc3MucHk=) | `100.00% <100.00%> (ø)` | |
| [plasmapy/formulary/drifts.py](https://codecov.io/gh/PlasmaPy/PlasmaPy/pull/865/diff?src=pr&el=tree#diff-cGxhc21hcHkvZm9ybXVsYXJ5L2RyaWZ0cy5weQ==) | `100.00% <100.00%> (ø)` | |
| [plasmapy/formulary/parameters.py](https://codecov.io/gh/PlasmaPy/PlasmaPy/pull/865/diff?src=pr&el=tree#diff-cGxhc21hcHkvZm9ybXVsYXJ5L3BhcmFtZXRlcnMucHk=) | `99.00% <100.00%> (+0.10%)` | :arrow_up: |
| [plasmapy/formulary/quantum.py](https://codecov.io/gh/PlasmaPy/PlasmaPy/pull/865/diff?src=pr&el=tree#diff-cGxhc21hcHkvZm9ybXVsYXJ5L3F1YW50dW0ucHk=) | `67.05% <100.00%> (+1.20%)` | :arrow_up: |

------

[Continue to review full report at Codecov](https://codecov.io/gh/PlasmaPy/PlasmaPy/pull/865?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/865?src=pr&el=footer). Last update [511c25b...8227650](https://codecov.io/gh/PlasmaPy/PlasmaPy/pull/865?src=pr&el=lastupdated). Read the [comment docs](https://docs.codecov.io/docs/pull-request-comments).

rocco8773: @anthony-vo As per our discussion on riot, we've decided that all aliases should be marked with the trailing underscore.  This will allow users to still defined the short named variable for their code use.  Basically, this is  how the user would interact with an alias...

```python
from plasmapy.formulary import va_

va = va_(...)
```

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/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