Introduction to Physics Applied to Medicine and Biology

Medical physics is the application of physics concepts and methods to the diagnosis and treatment of human diseases. Biophysics is an interdisciplinary field that applies approaches and methods traditionally used in physics to study biological systems.

Physics in Cancer Diagnosis and Treatment

Radiotherapy

The use of ionizing radiation to kill cancer cells and shrink tumors:

• Planning: Dose calculations and distribution
• Administration: Precise radiation control
• Safety: Protection of healthy tissues
• Efficacy: Treatment optimization

Medical Imaging

Techniques based on physics principles for diagnosis:

• X-rays: Tissue density images
• Computed Tomography (CT): Detailed 3D images
• Magnetic Resonance Imaging (MRI): Soft tissue images
• Positron Emission Tomography (PET): Functional images

Nanotechnology

Use of nanoparticles for targeted drug delivery:

• Targeted Delivery: Specific cancer cells
• Damage Minimization: Preservation of healthy cells
• Controlled Release: Controlled drug release
• Image Guidance: Real-time visualization

Cancer Biophysics

Cancer biophysics studies the physical properties of cancer cells and tissues:

• Cell Stiffness: Cancer cells are softer
• Metastasis: Facilitates spread to other organs
• Mechanical Properties: Changes in elasticity
• New Approaches: Physics-based diagnosis and treatment

Computational Applications

Radiotherapy Simulation

import numpy as np
import matplotlib.pyplot as plt
from scipy.ndimage import gaussian_filter

def simulate_radiation_dose(tumor_center, tumor_size, beam_center, beam_width):
    """Simulates radiation dose distribution."""
    
    # Create 3D grid
    x = np.linspace(0, 20, 100)
    y = np.linspace(0, 20, 100)
    z = np.linspace(0, 20, 50)
    
    X, Y, Z = np.meshgrid(x, y, z)
    
    # Define tumor
    tumor_mask = np.sqrt((X - tumor_center[0])**2 + 
                         (Y - tumor_center[1])**2 + 
                         (Z - tumor_center[2])**2) < tumor_size
    
    # Define radiation beam
    beam_mask = np.exp(-((X - beam_center[0])**2 + 
                         (Y - beam_center[1])**2) / (2 * beam_width**2))
    
    # Calculate dose
    dose = beam_mask * 100  # Maximum dose of 100 Gy
    
    # Apply depth attenuation
    for i in range(len(z)):
        dose[:, :, i] *= np.exp(-0.1 * z[i])
    
    return X, Y, Z, dose, tumor_mask

# Simulation parameters
tumor_center = (10, 10, 10)
tumor_size = 2
beam_center = (10, 10, 0)
beam_width = 1

# Run simulation
X, Y, Z, dose, tumor_mask = simulate_radiation_dose(
    tumor_center, tumor_size, beam_center, beam_width
)

# Visualize dose in central plane
plt.figure(figsize=(12, 4))

plt.subplot(1, 3, 1)
plt.contourf(X[:, :, 25], Y[:, :, 25], dose[:, :, 25], levels=20)
plt.colorbar(label='Dose (Gy)')
plt.title('Dose in Central Plane')
plt.xlabel('X (cm)')
plt.ylabel('Y (cm)')

plt.subplot(1, 3, 2)
plt.contourf(X[:, 50, :], Z[:, 50, :], dose[:, 50, :], levels=20)
plt.colorbar(label='Dose (Gy)')
plt.title('Dose in Sagittal Plane')
plt.xlabel('X (cm)')
plt.ylabel('Z (cm)')

plt.subplot(1, 3, 3)
plt.contourf(Y[50, :, :], Z[50, :, :], dose[50, :, :], levels=20)
plt.colorbar(label='Dose (Gy)')
plt.title('Dose in Coronal Plane')
plt.xlabel('Y (cm)')
plt.ylabel('Z (cm)')

plt.tight_layout()
plt.show()

Cell Mechanics Analysis

def analyze_cell_mechanics(elastic_modulus, cell_size, applied_force):
    """Analyzes cell mechanical properties."""
    
    # Hooke's Law: F = k * x
    # k = E * A / L (spring constant)
    
    # Cell area (assuming spherical shape)
    cell_area = 4 * np.pi * (cell_size/2)**2
    
    # Spring constant
    spring_constant = elastic_modulus * cell_area / cell_size
    
    # Deformation
    deformation = applied_force / spring_constant
    
    # Elastic energy
    elastic_energy = 0.5 * spring_constant * deformation**2
    
    return {
        'spring_constant': spring_constant,
        'deformation': deformation,
        'elastic_energy': elastic_energy
    }

# Compare normal vs cancer cells
normal_cell = analyze_cell_mechanics(1000, 10, 1)  # kPa, μm, nN
cancer_cell = analyze_cell_mechanics(500, 10, 1)   # Softer

print(f"Normal Cell - Constant: {normal_cell['spring_constant']:.1f} nN/μm")
print(f"Cancer Cell - Constant: {cancer_cell['spring_constant']:.1f} nN/μm")
print(f"Stiffness Ratio: {normal_cell['spring_constant']/cancer_cell['spring_constant']:.1f}")

Imaging Techniques

• Digital Radiography: High-resolution images
• Spiral CT: Fast volume acquisition
• Functional MRI: Brain activity images
• PET-CT: Anatomy and function combination

Learning Resources

• Books/Journals: Physics in Medicine & Biology
• Courses: MIT OpenCourseWare: radiation and nuclear science courses, Coursera medical physics search
• Tools: Geant4, MCNP, MATLAB

---

Understanding these physical properties can lead to new approaches for cancer diagnosis and treatment.