Skip to content

Instantly share code, notes, and snippets.

View JesusCastroFernandez's full-sized avatar
💭
Happy

JesusCastroFernandez

💭
Happy
View GitHub Profile
from qiskit import QuantumCircuit, Aer, execute
from qiskit.visualization import plot_histogram
# Definición del circuito cuántico
num_qubits = 1 # Número total de qubits
num_shots = 1000 # Número de ejecuciones del circuito
circuit = QuantumCircuit(num_qubits, num_qubits)
# Parámetros del estado cuántico |ψ⟩ = α|0⟩ + β|1⟩
from qiskit import QuantumCircuit, Aer, execute
from qiskit.visualization import plot_histogram
# Definición del circuito cuántico
num_qubits = 3 # Número total de qubits
num_shots = 1000 # Número de ejecuciones del circuito
circuit = QuantumCircuit(num_qubits, num_qubits)
# Preparación del estado inicial
import numpy as np
import scipy.constants as constants
# Función para calcular la probabilidad de encontrar al electrón en un intervalo dado
def calcular_probabilidad(intervalo, psi_funcion):
x = np.linspace(intervalo[0], intervalo[1], 1000)
probabilidad = np.abs(psi_funcion(x))**2
probabilidad_total = np.trapz(probabilidad, x)
return probabilidad_total
# Datos
me = mp = 9.11e-31 # Masa del electrón y del positrón en kg
e = 1.602e-19 # Carga elemental en culombios
h = 6.62e-34 # Constante de Planck en J·s
# Cálculo de la energía total producida
masa_total = me + mp # Masa total de los dos partículas en kg
energia_total = masa_total * c**2 # Energía total producida en julios
# Cálculo de la frecuencia de la radiación producida
# Datos
lambda1 = 200e-9 # Longitud de onda lambda1 en metros
potencial1 = 2 # Potencial de frenado para lambda1 en voltios
lambda2 = 240e-9 # Longitud de onda lambda2 en metros
potencial2 = 1 # Potencial de frenado para lambda2 en voltios
e = 1.6e-19 # Carga elemental en culombios
c = 3e8 # Velocidad de la luz en m/s
# Cálculo del trabajo de extracción del metal
trabajo_extraccion = potencial1 * e # Trabajo de extracción en julios
# a. Electrón con velocidad medida con una incertidumbre del 0,003%
velocidad_electron = 7000e3 # Velocidad del electrón en m/s
incertidumbre_velocidad_electron = 0.003 / 100 # Incertidumbre de la velocidad como fracción decimal
# Cálculo de la incertidumbre en la determinación de la posición del electrón
incertidumbre_posicion_electron = h / (2 * m_electron * velocidad_electron) * incertidumbre_velocidad_electron
# b. Partícula de 50 g que se desplaza a una velocidad de 300 m/s, medida con la misma incertidumbre
masa_particula = 50e-3 # Masa de la partícula en kg
velocidad_particula = 300 # Velocidad de la partícula en m/s
# Constantes
h = 6.62607015e-34 # Constante de Planck en J*s
m_electron = 9.10938356e-31 # Masa del electrón en kg
c = 299792458 # Velocidad de la luz en m/s
# a. Electrón acelerado por una diferencia de potencial (ΔV)
delta_V = 100 # Diferencia de potencial en voltios
energia_electron = delta_V * 1.602176634e-19 # Energía del electrón en Joules
lambda_electron = h / (m_electron * c) * (1 / energia_electron) # Longitud de onda asociada en metros
# Datos
longitud_onda = 1.5e-6 # Longitud de onda de la radiación incidente en metros
velocidad_fotoelectrones = 105 # Velocidad de los fotoelectrones emitidos en m/s
masa_electron = 9.1e-31 # Masa del electrón en kg
# Cálculo del trabajo de extracción
energia_cinetica = 0.5 * masa_electron * velocidad_fotoelectrones**2 # Energía cinética de los fotoelectrones en Joules
trabajo_extraccion = energia_cinetica # Suponiendo que no hay energía adicional (como potencial) a tener en cuenta
# Cálculo de la frecuencia umbral de fotoemisión
from qiskit import QuantumCircuit, ClassicalRegister, QuantumRegister, execute, Aer
# Función para calcular la energía máxima de los fotoelectrones
def calcular_energia_maxima(frecuencia, trabajo_extraccion):
h = 6.62607015e-34 # Constante de Planck en J*s
energia_foton = h * frecuencia # Energía de un fotón en Joules
energia_maxima = energia_foton - trabajo_extraccion # Energía máxima de los fotoelectrones en Joules
# Convertir la energía máxima a electronvoltios (eV)
energia_maxima_ev = energia_maxima / 1.602176634e-19
from qiskit import QuantumCircuit, ClassicalRegister, QuantumRegister, execute, Aer
# Función para calcular la velocidad de las OEM
def calcular_velocidad_ondas(l):
c = 299792458 # Velocidad de la luz en m/s
v = c / l # Calcular la velocidad de las ondas electromagnéticas
return v
# Función para calcular la energía de los fotones