diff --git a/Project.toml b/Project.toml
index b0527d9..85ed6db 100644
--- a/Project.toml
+++ b/Project.toml
@@ -10,6 +10,7 @@ Dates = "ade2ca70-3891-5945-98fb-dc099432e06a"
 PyFortran90Namelists = "e44308e6-bd5b-11e9-2850-49daf8f1ec40"
 QuantumESPRESSOBase = "51b62caa-b28f-11e9-38c2-1f67cb498e05"
 ReadableRegex = "cbbcb084-453d-4c4c-b292-e315607ba6a4"
+StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 VersionParsing = "81def892-9a0e-5fdd-b105-ffc91e053289"
 
 [weakdeps]
diff --git a/src/PWscf/output/each.jl b/src/PWscf/output/each.jl
index 87f483f..fc108ab 100644
--- a/src/PWscf/output/each.jl
+++ b/src/PWscf/output/each.jl
@@ -1,3 +1,5 @@
+using StaticArrays: SVector
+
 export eachstep,
     eachiteration,
     eachiterationhead,
@@ -6,6 +8,9 @@ export eachstep,
     eachunconvergedenergy,
     eachconvergedenergy,
     each_energy_by_step,
+    eachatomicforceblock,
+    eachatomicforce,
+    eachtotalforce,
     eachcellparameterscard,
     eachatomicpositionscard,
     eachtimeditem
@@ -283,6 +288,91 @@ eachconvergedenergy(str::AbstractString) =
 
 function each_energy_by_step end
 
+const FORCES_ACTING_ON_ATOMS_BLOCK = Regex(
+    raw"Forces acting on atoms (cartesian axes, Ry/au):" *
+    capture(lazy_zero_or_more(ANY)) *
+    rs"Total force =[ \t]*" *
+    capture(rs"([-+]?[0-9]*\.[0-9]+|[0-9]+\.?[0-9]*)") *
+    rs"[ \t]+Total SCF correction =[ \t]*" *
+    capture(rs"([-+]?[0-9]*\.[0-9]+|[0-9]+\.?[0-9]*)"),
+)
+
+struct EachAtomicForceBlock <: Each
+    iterator::Base.RegexMatchIterator
+end
+
+function Base.iterate(iter::EachAtomicForceBlock)
+    iterated = iterate(iter.iterator)
+    if isnothing(iterated)
+        return nothing
+    else
+        matched, state = iterated
+        return matched.match, state
+    end
+end
+function Base.iterate(iter::EachAtomicForceBlock, state)
+    iterated = iterate(iter.iterator, state)
+    if isnothing(iterated)
+        return nothing
+    else
+        matched, state = iterated
+        return matched.match, state
+    end
+end
+
+Base.eltype(::Type{EachAtomicForceBlock}) = String
+
+Base.IteratorSize(::Type{EachAtomicForceBlock}) = Base.SizeUnknown()
+
+eachatomicforceblock(str::AbstractString) =
+    EachAtomicForceBlock(eachmatch(FORCES_ACTING_ON_ATOMS_BLOCK, str))
+
+const FORCE_ACTING_ON_ATOM = Regex(
+    rs"atom\s+(\d+)\s+type\s+(\d+)\s+force\s+=\s+(-?\d+\.\d+)\s+(-?\d+\.\d+)\s+(-?\d+\.\d+)"
+)
+
+struct AtomicForce <: PWOutputItem
+    atom::Int64
+    type::Int64
+    force::SVector{3,Float64}
+end
+
+function Base.parse(::Type{AtomicForce}, str::AbstractString)
+    obj = tryparse(AtomicForce, str)
+    isnothing(obj) ? throw(ParseError("no matched string found!")) : return obj
+end
+function Base.tryparse(::Type{AtomicForce}, str::AbstractString)
+    matched = match(FORCE_ACTING_ON_ATOM, str)
+    if isnothing(matched)
+        return nothing
+    else
+        atom, type = map(Base.Fix1(parse, Int64), matched.captures[1:2])
+        force = map(Base.Fix1(parse, Float64), matched.captures[3:5])
+        return AtomicForce(atom, type, force)
+    end
+end
+
+eachatomicforce(str::AbstractString) = EachParsed{AtomicForce}(FORCE_ACTING_ON_ATOM, str)
+
+const TOTAL_FOCE = Regex(
+    rs"Total force =[ \t]*" * capture(rs"([-+]?[0-9]*\.[0-9]+|[0-9]+\.?[0-9]*)")
+)
+
+struct TotalForce <: PWOutputItem
+    force::Float64
+end
+
+function Base.parse(::Type{TotalForce}, str::AbstractString)
+    obj = tryparse(TotalForce, str)
+    isnothing(obj) ? throw(ParseError("no matched string found!")) : return obj
+end
+function Base.tryparse(::Type{TotalForce}, str::AbstractString)
+    matched = match(TOTAL_FOCE, str)
+    return isnothing(matched) ? nothing : TotalForce(parse(Float64, matched[1]))
+end
+
+eachtotalforce(str::AbstractString) = EachParsed{TotalForce}(TOTAL_FOCE, str)
+
 eachcellparameterscard(str::AbstractString) =
     EachParsed{CellParametersCard}(CELL_PARAMETERS_BLOCK, str)